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Abstract 



We present theoretical and observational studies of non-Gaussian fluctuations in the cosmic mi- 
crowave background (CMB) radiation anisotropy. We use the angular bispectrum and trispectrum, 
the harmonic transform of the angular three- and four-point correlation functions. If the primordial 
fluctuations are non-Gaussian, then this non-Gaussianity will be apparent in the CMB sky. 

Non-linearity in inflation produces the primordial non-Gaussianity. We predict the primary 
angular bispectrum from inflation down to arcminutes scales, and forecast how well we can measure 
the primordial non-Gaussian signal. In addition to that, secondary anisotropy sources in the low- 
redshift universe also produce non-Gaussianity, so do foreground emissions from extragalactic or 
interstellar microwave sources. We study how well we can measure these non-Gaussian signals, 
including the primordial signal, separately. We find that when we can compute the predicted form 
of the bispectrum, it becomes a "matched filter" for finding non-Gaussianity in the data, being 
very powerful tool of measuring weak non-Gaussian signals and of discriminating between different 
non-Gaussian components. We find that slow-roll inflation produces too small bispectrum to be 
detected by any experiments; thus, any detection strongly constrains this class of models. We also 
find that the secondary bispectrum from coupling between the Sunyaev-Zel'dovich effect and the 
weak lensing effect, and the foreground bispectrum from extragalactic point sources, give detectable 
non-Gaussian signals on small angular scales. 

We test Gaussianity of the COBE DMR sky maps, by measuring all the modes of the angular 
bispectrum down to the DMR beam size. We compare the data with the simulated Gaussian 
realizations, finding no significant signal of the bispectrum on the mode-by-mode basis. We also find 
that the previously reported detection of the bispectrum is consistent with a statistical fluctuation. 
By fitting the theoretical prediction to the data for the primary bispectrum, we put a constraint 
on non-linearity in inflation. Simultaneously fitting the foreground bispectra, which are estimated 
from interstellar dust and synchrotron template maps, shows that neither dust nor synchrotron 
emissions contribute significantly to the bispectrum at high Galactic latitude. We thus conclude 
that the angular bispectrum finds no significant non-Gaussian signals in the DMR data. 

We present the first measurement of the angular trispectrum on the DMR sky maps, further 
testing Gaussianity of the DMR data. By applying the same method as used for the bispectrum 
to the DMR data, we find no significant non-Gaussian signals in the trispectrum. Therefore, the 
angular bispectrum and trispectrum show that the DMR sky map is comfortably consistent with 
Gaussianity. 



The methods that we have developed in this thesis can readily be applied to the MAP data, 
and will enable us to pursue non-Gaussian CMB fluctuations with the unprecedented sensitivity. 
We show that high-sensitivity measurement of the CMB bispectrum and trispectrum will probe the 
physics of the early universe as well as the astrophysics in the low-redshift universe, independently 
of the CMB power spectrum. 
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Chapter 1 

Introduction 



1.1 Why Pursue Non-Gaussianity in CMB? 

Modern understanding to emergence of inhomogeneity in the universe through the cosmic history 
is outstanding; the inhomogeneity is quantum in origin. Then, it becomes classical through its evo- 
lution, producing fluctuations in the cosmic microwave background (CMB) radiation, and creating 
complex structures, such as galaxies, seen in the present universe. 

A theory of early universe, the cosmic inflation ( puth, 1981 ; Sato, 1981 ; Albrecht and Stein 



hardt, 1982; Linde, 1982), has predicted the emergence of the quantum fluctuations in early uni- 
verse. Inflation not only resolves several serious issues in the old Big-Bang cosmology, but also gives 
a mechanism to produce inhomogeneity in the universe, and makes specific testable predictions for 



the global structure and the inhomogeneity in the universe ( puth and Pi, 1982 ; Hawking, 1982 



|Starobinsky, 198"2 ; Bardeen et al., 1983| ). The predictions may concisely be summarized as follows: 



(a) The observable universe is spatially flat. 

(b) The observable universe is homogeneous and isotropic on large angular scales, apart from tiny 
fluctuations. 

(c) The primordial inhomogeneity has a specific spatial pattern, so-called the scale-invariant 
fluctuation. 

(d) Statistics of the primordial inhomogeneity obey Gaussian statistics. 

All these predictions but (d) have passed challenging observations. Among those observations, 
the most firm evidence supporting the predictions has come from measurement of the angular 
distribution of CMB. Being the oldest observable object in the universe, CMB is the best clue to 
early universe. 

Observed isotropy of CMB with 0.001% accuracy^ supports the prediction (b) ( Mather et al. 



1990| ; pmoot et al., 1992[) . The measured angular distribution of the CMB anisotropy, 0.001% 



*Apart from anisotropics due to the local motion of the Earth, which is of order 0.1% (Smoot et al., 1991) 
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inhomogeneity, supports the primordial fluctuation distribution being consistent with the prediction 
(c) in the flat universe, the prediction (a) (Bennett et al., 1996; Miller et al., 1999; |de Bernadis et 
al., 20001 ; [Hanany et al., 200C| ). 

Using various statistical techniques, many authors have attempted to test the prediction (d), 
the Gaussianity of the primordial inhomogeneity, using the CMB anisotropy on large angular scales 
(~ 7°) QKogut et al., 1996b| ; [Ferreira et al., 1998| ; pando et al., 1998| ; |Bromley and Tegmark, 1999| ; 
Banday et al., 200C ; Mukherjee et al., 200C ; Magueijo, 2000 ; Sandvik and Magueijo, lOOOl ; | arreirq 
|et al., 20"00 ), on intermediate scales (~ 1°) ( Park et al., 2001] ), and on small scales (~ 10') ( |Wu et 



al., 20011) . 

In contrast to the predictions (a)-(c) for which the observations suggest no controversy, previous 
work on searching for non-Gaussian CMB anisotropics has come to different conclusions from one 
another: some do claim detection ( |Ferreira et al., 199S ; Pando et al., 1998; Magueijo, "2000D , and 
the others do not ( Kogut et al., 1996b| ; [Sandvik and Magueijo, "2000t |Barreiro et al., 200q ; Park et| 
al., 2001| ; |Wu et al., 200l| ). Furthermore, [Bromley and Tegmark (1999| ) and Panday et al. (200q ) 
claim the non-Gaussian signal of Ferreira et al. (1998| ) to be non-cosmological in origin. Mukherjee] 
et al. (200C| ) revise the conclusion of pando et al. (1998|) by addressing ambiguity in their method, 
from which the conclusion crucially depends upon an orientation of the data on the sky, and show 
no evidence for the non-Gaussianity. The existence of non- Gaussianity in CMB is controversial. 

Since non-Gaussianity has infinite degrees of freedom, testing a Gaussian hypothesis is difficult; 
one statistical method showing CMB consistent with Gaussian does not mean CMB being really 
Gaussian. In this sense, different statistical methods can come to different conclusions. In addition 
to the difficulty, cosmological non-Gaussianity is hard to measure. Instrumental and environmen- 
tal effects in observations easily produce spurious non-Gaussian signals. Astronomical microwave 
sources such as interstellar dust emissions also produce strong non-Gaussian fluctuations on the 
sky. Hence, we must be as careful as possible when searching for cosmological non-Gaussianity. 

Pursuit of cosmological non-Gaussianity in CMB is a challenging test of inflation, the origin 
of inhomogeneity in the universe. In this dissertation, we present theoretical and observational 
studies of the CMB non-Gaussianity. Our primary goal is to test inflation with non-Gaussian CMB 
fluctuations. 

As a statistical tool of searching for non-Gaussianity, we use the angular three- and four-point 
correlation functions in harmonic space, the angular bispectrum and trispectrum, which are sensitive 
to weakly non-Gaussian fluctuations. An advantage of these angular n-point harmonic spectra over 
other statistics is that they are predictable from not only inflation, but also secondary sources in 
the low-redshift universe (the Sunyaev-ZeFdovich effect, weak gravitational lensing effect, and so 
on), extragalactic foreground emissions, global topology of the universe, and so on. When we fit 
those predictions to the measured spectra, the predictions become "matched filters" for detecting 
weak non-Gaussianity in the data, and much more powerful than just a null test of Gaussianity. 

This dissertation is organized as follows. The rest of this chapter will overview what the CMB 
sky looks like, and the previous study of CMB non-Gaussianity, followed by a heuristic description 
of non-Gaussian fluctuation production in inflation. 
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In chapter ^J, we go through primordial fluctuation generation mechanism in inflation, and in- 
vestigate how non-linear curvature perturbations are generated and propagated through the CMB 
anisotropy. We present a possible quantum-to-classical transition mechanism of the quantum fluc- 
tuations on super horizon scales. 

In chapter [3|, we study statistical properties of the angular n-point harmonic spectrum for n = 2 
(power spectrum), 3 (bispectrum), and 4 (trispectrum). We present practical methods of measuring 
the angular bispectrum and trispectrum from observational data with many pixels. 

In chapter ||, we make theoretical predictions for the CMB angular bispectrum, which include 
primary contribution from inflation, secondary contribution from the Sunyaev-Zel'dovich effect 
and the weak-lensing effect, and foreground contribution from extragalactic radio and infrared 
astronomical sources. We estimate signal-to- noise ratios of detecting each contribution with CMB 
experiments, further discussing how well we can measure each contribution separately. 

In chapter [|, we measure the angular bispectrum on the COBE DMR four-year data, testing 
Gaussianity of the data. By fitting the theoretical prediction for the primary bispectrum to the data, 
we constrain non-linearity in inflation; also fitting foreground bispectra from interstellar microwave 
emissions takes into account the effect of non-Gaussian interstellar dust and synchrotron emissions 
at high Galactic latitude. 

In chapter [| we measure the angular trispectrum on the DMR data, further testing Gaussianity 
of the data. 



1.2 Inhomogeneity in Microwave Sky 

The cosmic microwave background radiation (CMB) is the isotropic microwave radiation filling 
the sky. The temperature is precisely measured to be 2.73 K ( [Mather et al., 199C| ), and the peak 
intensity at v = 160 GHz (A = 1.9 mm) is 370 MJy str" 1 . What if subtracting this mean radiation 
from the sky, what are we left with, completely dark, literally nothing? 



1.2.1 COBE DMR sky maps 



The Differential Microwave Radiometer (DMR) mounted on the Cosmic Background Explorer 
(COBE), the satellite for full sky measurement of CMB launched by NASA in 1989, has revealed 



that what we are left with is the tiny inhomogeneity; the r.m.s. amplitude is about 30 /iK ( Smoot 



et al., 1992 ), 0.001% of the mean temperature. 

Various statistical analyses on the DMR sky maps have shown the angular distribution of 
the CMB anisotropy remarkably consistent with the scale-invariant fluctuation, the prediction of 
inflation flBennett et al., 1996| |G6rski et al., 19961 ; [Hinshaw et al., 1996| ; Wright et al., 1996[ ). 
The scale-invariant fluctuation implies that the r.m.s. amplitude of the CMB anisotropy is nearly 
independent of angular scales. DMR has measured 35 fiK r.m.s. fluctuations on 7° scale, 29 jiK 
on 10° scale ( [Banday et al., 19971 ). 

DMR comprises 3 dual-horn antennas working at 31.5, 53, and 90 GHz. The combination of 
53 and 90 GHz maps gives the most sensitive sky map to CMB, while 31.5 GHz map is twice as 
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Figure 1.1: COBE DMR Sky Map 
COBE DMR four-year 53 + 90 GHz sky map in Galactic projection, using the HEALPix 
pixelization (Gorski et al., 199S) with 1?83 pixel size, leaving 12,288 pixels. The left is the 
raw map, while the right map has been smoothed with a 7° FWHM Gaussian. 



noisier as the other channels. Figure 1.1 shows the combined DMR full sky map. The left shows 
the raw map in which the instrumental noise dominates appearance of the map; the right shows the 
smoothed map with the DMR beam in which the instrumental noise is filtered out, giving better 
appearance of CMB signals. The mean signal-to-noise ratio of hot and cold spots in the smoothed 
map is 2, while a few prominent spots have 3 to 4. Hence, we cannot say much about structures of 
the CMB anisotropy relying on the map basis; however, we can do say on the statistical basis. 

Statistically, structures in the DMR map are inconsistent with pure instrumental noise; on the 
contrary, the structure has a distinct angular correlation pattern represented by the scale-invariant 
fluctuation. To quantify this, it is useful to calculate the angular power spectrum, Ci, the harmonic 
transform of the angular two-point correlation function, which measures how much fluctuation 
power exists on a given angular scale, 9 ~ tt/Z. Figure L2 plots the measured C\ on the DMR 
map. What is actually plotted is 1(1 + 1)Ci/2tt, roughly mean squares of fluctuations at I. The 
scale-invariant fluctuation implies C\ oc [/(/ + 1)] 1 ( Peebles, 1982 ), and hence 1(1 + remains 
constant (solid line), so do the data points in the figure; the data points fit the inflation's prediction 
well. The dashed line shows a more accurate prediction, taking into account the effects of general 
relativistic photon-baryon fluid dynamics before the decoupling as well as of time evolution of 
gravitational potential field after the decoupling. The agreement with the data becomes better, 
further confirming that the DMR angular power spectrum is consistent with inflation. 



Gaussianity of the DMR data has been tested with various statistical methods (Kogut et al. 



1996b; Ferreira et al., 1998; Pando et al., 1998]; Bromley and Tegmark, 1999 


; Banday et al., 2000fc 


Mukherjee et al., 2000; 


Magueijo, 2000; 


Sandvik and Magueijo, 2000; 


Barreiro et al., 2000). 


Ferreira 



|et al. (1998 ) and Magueijo (2000| ) claim positive detection of non-Gaussian signals using the angular 
bispectrum, and Pando et al. (1998) claim detection using the wavelet analysis. The latter non- 
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Figure 1.2: COBE DMR Angular Power Spectrum 

The CMB angular power spectrum, Ci, measured on the COBE DMR four-year map. 
The plotted quantity, 1(1 + l)C;/(27r), represents the r.m.s. amplitude of fluctuations 
at a given angular scale, I. The data points (filled circles) are uncorrelated with each 
other (Tegmark and Hamilton, 1997). The solid line shows the scale-invariant power 
spectrum, 1(1 + 1)C/ = constant, while the dashed line shows a ACDM spectrum. 
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Gaussian signal has appeared to be less significant than they claim, as it disappears when the 
DMR map is rotated by 180° (Mukherjee et al., 200C; Barreiro et al., 2000p . For the former two 



bispectrum analyses, Bromley and Tegmark (1999| ) and Banday et al. (20"0(i| ) claim that the Ferreiraj 



eFli37Jl998|) 's signal is non-cosmological, but their claims do not account for the Magueijo "(2000|) , s 



signal. In this thesis, we will argue that the reported non-Gaussian signals are not a matter of 
origin, but statistical fluctuations. 

1.2.2 Post- COBE era 

After the discovery of COBE, pursuit of the CMB anisotropy has been oriented toward measurement 
of the angular power spectrum, Ci, on smaller angular scales, i.e., larger I. Particularly, many efforts 
have been made to measure Ci at I ~ 200 (9 ~ 1°), where inflation predicts a prominent peak in 



/(/ + 1)C/ as a consequence of flatness of the universe, the prediction of inflation (|Kamionkowski et 
|aL, 1994) . 



By early 2000, there has been strong evidence for the peak ( Miller et al., 1999| ); in the end 
of 2000, BOOMERanG and MAXIMA, balloon-borne CMB experiments, have detected the peak 
( de Bernadis et al., 2000 ; Hanany et al., 2000D , further supporting inflation. Figure L3 compares 



the data from three experiments, which probe different angular scales from one another, with a 
prediction from inflation. The agreement between the data and the prediction is outstanding. 

Measurement of Ci so far, however, has assumed CMB Gaussian. If CMB is not Gaussian, 
then the measured Ci is biased. Moreover, when we fit the measured C; to a theoretical Ci, 
we need to know the covariance matrix of C\. Since the covariance matrix of C\ is the four-point 
harmonic spectrum, the trispectrum, we have to investigate non-Gaussian signals in the trispectrum 
to construct the covariance matrix accurately. 

Even if non-Gaussianity is small, we have to take it into account in analyzing Ci; the next 
generation satellite experiments, MAP and Planck, will measure C\ with 1% or better accuracy, 
and we will use the measured C\ to determine many of cosmological parameters with 10% or better 
accuracy. Unless non-Gaussian effect is much smaller than the observational uncertainty (who 
knows?), we have to take the effect into account, to achieve the accurate measurement of the 
parameters. 

The state-of-the-art balloon-borne experiments provide not only C\, but also high resolution, 
high signal-to- noise ratio CMB maps. The left of figure |1.4| shows the BOOMERanG sky map, 
which covers roughly 3% of the sky, 1,800 square degrees, with 10' (0?17) angular resolution. The 
high signal-to-noise ratios in the map show the observed structures in the map not instrumental 
noise, but CMB. 

Visually, the CMB anisotropy in the map looks very much Gaussian, implying non-Gaussianity 
is weak, if any. As yet the BOOMERanG team has not performed Gaussianity tests on the data; 
however, the MAXIMA team has done on their data. The MAXIMA map covers 124 square degrees 
of the sky with 10' angular resolution, leaving 5,972 pixels. They have found that the one-point 
probability density distribution function and the Minkowski functionals measured on the map are 
consistent with Gaussianity. 
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Figure 1.3: COBE DMR, QMASK, and BOOMERanG Angular Power Spectra 

The CMB angular power spectrum, measured by the three experiments probing different 
angular scales. The circles are the COBE DMR data (full sky coverage with a Galactic cut, 
7° angular resolution), the triangles are the QMASK data (648-square-degree sky coverage, 
0?68 angular resolution), and the squares are the BOOMERanG data (1,800-square-degree 
sky coverage, 0?17 angular resolution). The dashed line shows a prediction of inflation, 
a ACDM spectrum with il m = 0.3, ft A = 0.7, fl h = 0.04, h = 0.7, and n = 0.95. 
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Figure 1.4: Sub-degree CMB Maps: BOOMERanG and QMASK 

The left is the BOOMERanG sky map, which covers 3% of the sky, 1,800 square degrees, 
with 10' (0?17) angular resolution (dc Bcrnadis ct al, 200C). The beam size is indicated 
by the filled circle at the bottom-right corner. The right is the QMASK sky map, which 
covers 648 square degrees of the sky with 0?68 angular resolution (Xu et al, 2001). Both 
maps have been scaled to match each other in size. The QMASK map is one-third of the 
BOOMERanG map. 



The genus statistic has been applied to the QMASK map ( Xu et al., 2001| ), a high signal-to-noise 
ratio CMB map combining the balloon-borne QMAP ( Devlin et al., 1998| ; Herbig et al., 1998 : dc 
Oliveira-Costa et al., 1998] ) and the ground-based Saskatoon experiments (Netterfield et al., 1997), 



showing the map consistent with Gaussian (Park et al., 2001). The QMASK map is shown on the 
right of figure 1.4, where the size of the map has been scaled to match the BOOMERanG map. 
The QMASK map covers 648 square degrees of the sky, one-third of the BOOMERanG map, with 
0?68 angular resolution, leaving 6,495 pixels. 

The number of independent pixels that has been used for Gaussianity tests on the DMR (3,881), 
MAXIMA (5,972), and QMASK (6,495) data is similar to each other; although their angular scales 
are quite different, they have similar power of testing Gaussianity. We will have a significant progress 
when the BOOMERanG data, in which the number of pixels is 57,103, test the Gaussianity. 

The just-launched CMB satellite, the Microwave Anisotropy Probe (MAP), will increase our 
power of pursuing non-Gaussianity substantially. It will survey the full sky with the angular 
resolution better than 14', and have £ 10 6 pixels. The statistical methods developed in this thesis 
can readily be applied to the MAP data, enabling us to test the Gaussianity with unprecedented 
sensitivity. 



1.3 Non-Gaussian Fluctuations in Inflation 
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1.3.1 Adiabatic production of non-Gaussian fluctuations 

While inflation predicts Gaussian CMB fluctuations to very good accuracy, strictly speaking, non- 
linearity in inflation produces weakly non-Gaussian fluctuations, which propagates through CMB. 
Although the exact treatment is complicated, we present a basic idea behind it concisely here 
( Komatsu and Spcrgel, 2001b] ). 



The curvature perturbations, generate the CMB anisotropy, AT/T. The linear perturbation 
theory gives a linear relation between $ and AT/T, 

AT , . 

— ~5T$, (1.1) 

where gx is the radiation transfer function. For temperature fluctuations on super-horizon scales 



at the decoupling epoch, the Sachs- Wolfe effect (Sachs and Wolfe, 1967) dominates, and = —1/3 



for adiabatic fluctuations. On sub-horizon scales, oscillates (acoustic oscillation), and we need 
to solve the Boltzmann photon transport equations coupled with the Einstein equations for g^. 
It follows from the relation, AT oc <3?, that AT is Gaussian, if is Gaussian. As we will see, 
non-linearity in inflation makes <3? weakly non-Gaussian. 

Even if $ is Gaussian, AT/T can be non-Gaussian. According to the general relativistic cos- 
mological perturbation theory, there is a non-linear relation between AT/T and <£: 

AT 

~Y ~ 9T ($ + h$ 2 ) • (1.2) 

Here, the second term with a coefficient of order unity, /$ ~ 0(1), is the higher-order correction 
arising from the second-order perturbation theory ( [Pyne and Carroll, 1996 ). It produces non- 



Gaussian fluctuations; thus, even if $ is Gaussian, AT becomes weakly non-Gaussian. 

Is <& Gaussian? Non-linearity in inflation makes <3? weakly non-Gaussian. By expanding the 
fluctuation dynamics in inflation up to the second order, we obtain a non-linear relation between 
$ and inflaton fluctuations, 5<j>: 

$ ~ m£g 9 (54> + m^fsjStf) . (1.3) 



Salopek and Bond (1990| ) show that this relation is a non-linear solution for curvature perturbations 



on super horizon scales; the solution gives g$ ~ 0(10) and ~ 0(1O _1 ) for a class of slowly-rolling 
single-field inflation models. 

Quantum fluctuations produce Gaussian 5<j). If the dynamics of 5<p is simple enough to keep 
itself Gaussian throughout the evolution, then we can stop our consideration here; however, it is not 
necessarily true. For example, non-trivial interaction terms in the equation of motion for inflaton 
fields ( |Fark et al., 1993| ), or a non-linear coupling between long- wavelength classical fluctuations 
and short-wavelength quantum fluctuations in the context of stochastic inflation ( [Starobinsky, 1986 : 



Gangui et al., 1994| ), can make S(j) weakly non-Gaussian, resulting in a non-linear relation between 



and a Gaussian field, rj, 

W ~ 95^ (v + m^ 1 f v v 2 ) , (1.4) 
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Figure 1.5: Three Steps for Non-Gaussian CMB Anisotropy 
Adiabatic generation of non-Gaussian CMB anisotropics. First, inflation generates Gaussian 
quantum fluctuations, which become non-Gaussian inflaton fluctuations through a non- 
linear coupling between them (Eq.fll.4|)). Then, the inflaton fluctuations become more non- 



Gaussian curvature perturbations through a non- linear relation between them (Eq.(1.3)). 
Finally, the curvature perturbations become more non-Gaussian CMB anisotropics through 
non-linear gravitational effects (Eq.( |l.2| )). 
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where rj represents initially produced quantum fluctuations, and gg^ ~ 1 and f v ~ O(10 



Figure 1.5 summarizes the above three steps in the opposite order. 

Collecting all the above contributions, we obtain a non-linear relationship between AT/T and 

AT 



9T 



$ L + (u + <?i% + g^a^fv) *l] , (1-5) 



where $l = g^g&<t>'m~p\ r l ~ 10m~^rj is an auxiliary Gaussian curvature perturbation. It may be 
useful to define a non-linear coupling parameter, /nl = /$ + g® 1 fs<j> + 9^ 9ss fv ^he ^ rs ^ term 
i n /nLj the second order gravity effect ~ 0(1), is dominant compared with the other two terms 
~ 0(1O~ 2 ), non-linearity in slow-roll inflation. Note that /nl corresponds to —$3/2 in |Gangui et al.| 
|(1994f) and — a<j> in |Verde et al. (2000| ). Using /nl, we rewrite equation ( |l.5|) as AT(x)/T ~ ^t^Kx), 



where 

%2 /„\ /-t,2 



$(x) = c& L (x) + /nl [$£(x) - ($ L (x))j , (1.6) 

the angular bracket denoting the statistical ensemble average. 

To see intuitively what non-Gaussian fluctuations that we have considered here look like, in 



figure 1.6 we plot one-point probability density distribution function (p.d.f) of the CMB anisotropy. 
We compare Gaussian p.d.f with non-Gaussian p.d.f of adiabatic fluctuations produced in inflation. 
The dashed line is Gaussian distribution, i.e., no non-linear perturbations are included (/nl = 0). 
The solid line includes non-linear coupling of order /nl = 1000, while the dotted line of order 
/nl = 5000. It follows from the figure that positive /nl gives negatively skewed p.d.f; negative /nl 
gives positively skewed p.d.f. The larger |/nl| is, the more skewed p.d.f becomes. 

We find that /nl = 1000 gives virtually identical p.d.f to the Gaussian p.d.f. In chapter we 
will show that even with the ideal CMB experiment, we can measure |/nl| no smaller than 60, if 
we use the skewness of the one-point p.d.f. If we use the bispectrum, however, we can measure 
as small |/nl| as 3. The bispectrum is thus much more sensitive to the adiabatic non-Gaussianity 
than the skewness of the one-point p.d.f is. In chapter [|, we will show that the angular bispectrum 
measured on the COBE DMR sky map constrains |/nl| < 1-6 x 10 3 . The next generation satellite 
experiments, MAP and Planck, will improve the constraint substantially. 

Since the minimum |/nl| detectable with CMB experiments is 3, the second order gravity effect, 
/nl ~ C(l)> wn l n °t produce detectable non-Gaussianity in CMB, nor will slow-roll inflation, 

/ NL ~ 0(1O- 2 ). 

Yet, inflation may not be so simple. Any significant deviation from slow-roll, or features in 



an inflaton potential ( Kofman et al., 1991|) , could produce a bigger |/nl|> bigger non-Gaussianity. 



While we have restricted ourselves to adiabatic fluctuations, the deviation from Gaussianity can be 
more significant if inflation produces non-negligible isocurvature fluctuations ( [Linde and Mukhanov 



1997: Bucher and Zhu, 1997; Peebles, 1997). Measurement of non-Gaussian CMB anisotropics thus 



potentially constrains non-linearity, "slow-rollness" , and "adiabaticity" in inflation. In the next 
subsection, we concisely describe a possible mechanism to produce isocurvature fluctuations in 
inflation. 
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Figure 1.6: Adiabatic Non-Gaussianity: One-point p.d.f 
One-point probability density distribution function (p.d.f) of the CMB anisotropy, AT/T, 
comparing a Gaussian p.d.f with non-Gaussian p.d.f 's of non-linear adiabatic fluctuations 
produced in inflation. The dashed line plots Gaussian distribution. The solid line plots non- 
Gaussian distribution for /nl = 1000, the dotted line for /nl = 5000. The larger /nl is, the 
more negatively skewed p.d.f becomes. If /nl < 0, then p.d.f becomes positively skewed. 
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1.3.2 Isocurvature fluctuations 

Isocurvature fluctuations do not perturb spatial curvature at the initial fluctuation-generation 
epoch. In inflation, in addition to a scalar field responsible for adiabatic fluctuations, another scalar 
field, a, may produce isocurvature density fluctuations with amplitude of 5p a /p a ~ H 2 / (da / dt) , 
where H is the Hubble parameter during inflation. This formula assumes that a rolls down on 
its potential very slowly. In some cases, this fluctuation amplitude is about the same as adiabatic 
density fluctuations generated by a scalar field, cj), which drives inflation: 5pa-/pa- ~ Sp^/p^ ~ 
H 2 1 (dcj) I dt) . This happens when both fluctuations are produced in a similar way, through the 
quantum-fluctuation production in inflation. 

Even if dp/p is similar to each other, the energy density, p, can be significantly different. Since 
4> drives inflation, its energy density, p^, dominates the total energy density of the universe during 
inflation: p^ 3> p a \ thus, it gives 5p a ~ bp^ (p^j p^) <C Sp^. Then, the density fluctuations generate 
the curvature perturbations, <£. Since 5p a <C Sp^, a makes negligible contribution to compared 
with (p, i.e., <5/9(j does not generate the curvature perturbations, being an isocurvature mode. In 
this model, 5p a is Gaussian, as the quantum fluctuations have produced it linearly. 

If a moves fast, then the quantum fluctuations produce 5p a non-linearly; we have non-Gaussian 
density fluctuations. [Linde and Mukhanov (1997 ) have proposed a massive- free field oscillating 



about its potential minimum as a possible non-Gaussian isocurvature-fluctuation production mech- 
anism in inflation. The idea is as follows. When a field rolls down on a potential, V(a), very slowly, 
quantum fluctuations of a, 5a, produce the energy density fluctuations of 5p a ~ (dV/da)5a; thus, 
5p a is linear in 5a, being Gaussian. In contrast, when a field oscillates rapidly about a = 0, there is 
no mean field; for example, a massive-free scalar field with a potential V(a) = m 2 a 2 /2 for m ^ H 
produces the density fluctuations of 

5p a ~ m 2 a5a + m 2 (5af = m 2 {5af. (1.7) 

Here, m ^ H ensures that a has rolled down to a = quickly, and oscillates. Hence, 5p a is quadratic 
in 5a, being non-Gaussian. 

After the initial generation of isocurvature fluctuations, a may produce the curvature pertur- 
bations through the evolution. If a does not decay, or decays only very slowly, the energy density 
decreases as a~ 3 . On the other hand, the radiation energy density that is produced during the 
reheating phase by a decaying scalar field 4> that has driven inflation decreases as a -4 ; thus, at 
some point in the cosmic evolution, the cr-field energy density dominates the universe, producing 
the curvature perturbations, 

d>(x)=r ? 2 (x)-(r ? 2 (x)\, (1.8) 



and hence the CMB anisotropics, AT/T ~ gT&- Here, r\ is a Gaussian fluctuation field which is 
related to 5a, and <?t is the isocurvature radiation transfer function. The Sachs-Wolfe effect gives 
9T = -2. 

The CMB experiments show that isocurvature fluctuations do not contribute to the curvature 
perturbations very much; on the contrary, their contribution is negligible compared with adiabatic 
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contribution. Figure 1.3 compares a prediction for the CMB angular power spectrum from adiabatic 
fluctuations with the data. The agreement is very good, and there is no need to invoke isocurvature 
fluctuations. Moreover, the isocurvature fluctuations predict a very different form of the power 
spectrum; thus, the data have excluded possibility of the isocurvature fluctuations dominating the 
observed CMB power spectrum at high significance. 

Yet, there could exist isocurvature fluctuations in inflation. Generally speaking, if there are 
many scalar fields, there must exist isocurvature fluctuations. It is rather unusual to assume only 
one scalar field during inflation, for currently viable theories of the high energy particle physics 
predict existence of many kinds of scalar fields in a very high energy regime. There is, however, little 
hope to detect their signatures in the CMB power spectrum, as they are so weak compared with 
adiabatic fluctuations. Instead, searching for non-Gaussian signals in CMB is a promising strategy 
to look for some of those isocurvature fluctuations which are generally much more non-Gaussian 
than the adiabatic fluctuations. 

Since AT/T is quadratic in a Gaussian variable, one-point p.d.f of AT/T is the x 2 distribution 
with one degree of freedom. Figure |l.7| plots the one-point p.d.f of the isocurvature model (solid 
line) in comparison with Gaussian p.d.f (dashed line). The predicted p.d.f is highly non-Gaussian. 
If we assume the isocurvature CMB fluctuations dominating the universe, then the predicted non- 
Gaussian p.d.f may look too non-Gaussian to be consistent with observations; however, the COBE 
DMR data do not exclude this model on the basis of the non-Gaussianity because of the large 
beam-smoothing effect (Novikov et al., 2000). The bigger the beam is, the closer the smoothed 
X 2 distribution is to Gaussian distribution (Novikov et al., 2000). The CMB experiments probing 
much smaller angular scales than DMR will test the isocurvature non-Gaussian models. 

While we do not explore the isocurvature fluctuations so extensively, we present in appendix |C| 
an analytic prediction for the CMB angular bispectrum generated from the isocurvature fluctuations 
that we have described in this section. This formula may be used to fit the measured bispectrum; 
by doing so, we can constrain the model independently of the angular power spectrum. Although 
the one-point p.d.f of the isocurvature fluctuations is very similar to Gaussian distribution for large- 
beam CMB experiments, the bispectrum may still be powerful enough to detect the non-Gaussian 
signals. 
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Figure 1.7: Isocurvature Non-Gaussianity: One-point p.d.f 

One-point probability density distribution function (p.d.f) of the CMB temperature 
anisotropy, AT/T, comparing a Gaussian p.d.f with a non-Gaussian p.d.f of isocurva- 
ture fluctuations produced in inflation ( Linde and Mukhanov, 1997 ). The dashed line 
plots Gaussian distribution, while the solid line plots the non-Gaussian distribution. Note 
that we have assumed no beam smoothing; the large beam smoothing makes the non- 
Gaussian distribution similar to the Gaussian distribution (Novikov ct al., 2000). 



Chapter 2 

Perturbation Theory in Inflation 



2.1 Inflation — Overview 

During inflation, the universe expands exponentially. It implies the Hubble parameter, H(t) = 
dlna/dt, the expansion rate of the universe, being nearly constant in time, and the expansion scale 
factor, a(t), given by 

a(t) = a(t ) exp H(t')dt^j « a{t ) exp [H(t) (t - t )] . (2.1) 

The exponential expansion drives the observable universe spatially flat, for as the universe expands 
rapidly, a small section on a surface of a three-sphere of the universe approaches flat (we live on 
the section). Thus, inflation predicts flatness of the universe, and recent CMB experiments have 
confirmed the prediction ( Miller et al., 1999| ; de Bernadis et al., 2000 ; Hanany et al., 2000| ). 



What makes the exponential expansion possible? One finds that neither matter nor radiation 
can make it; on the contrary, their energy density, p, and pressure, p, make the universe decelerate. 
Since the universe accelerates only when p + 3p < 0, one needs a negative pressure component 
dominating the universe. How can it be possible? 

A spatially homogeneous scalar field, (ft, with a potential, V((p), provides negative pressure, 
making the exponential expansion possible. The energy density is p<f, = \{d<p/dt) 2 + V((f)), while 
the pressure is p^ = ^(dcfi/dt) 2 — V((j)), giving 

P4> + 3p^ = 2 [{d(j)/dtf - V(4>)} . (2.2) 

Hence, one finds that (d(p/dt) 2 < V(4>) suffices to accelerate the universe. This slowly-rolling scalar 
field is a key ingredient of inflation; by assuming a slowly-rolling scalar field dominant in early 
universe, the universe expands exponentially. 

While what is <fi and how it comes to dominate the universe are still in debate, a simple model 
sketched in figure ^Tl| works well. In the phase (a), <f> rolls down on V{4>) slowly, driving the 
universe to expand exponentially. In the phase (b), eft oscillates rapidly, terminating inflation. 
After inflation ends, interactions of 4> with other particles lead 4> to decay with a decay rate of T^, 
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Figure 2.1: Classical Evolution of Scalar Field 

Classical evolution of a scalar field, (j>, in a potential, V(cf>). (a) cf> rolls down on V(cf>) 
slowly, driving the universe to expand exponentially, (b) (j) oscillates rapidly about = 0, 
terminating inflation; it then decays into particles and radiation because of interactions, 
reheating the universe. 



producing particles and radiation. This is called a reheating phase of the universe, as 4> converts 
its energy density into heat by the particle production. A reheating temperature amounts to on 
the order of fF^mpi) 1 / 2 . While a precise value of reheating temperature depends upon models, it is 
typically on the order of 10 14-16 GeV ~ 10 27-29 K. Note that the smallness of the observed CMB 
anisotropy implies that <fi is coupled to other particles only very weakly, i.e., < H, giving lower 
reheating temperature. After the reheating, radiation dominates the universe, and the Big-bang 
scenario describes the rest of the cosmic history. 

A class of inflation models with the potential sketched in figure 2A is called the chaotic inflation, 
for which V{(p) oc (j) n ( Linde, 1983 ). Until now, this model has remained the most successful 
realization of inflation with broad applications (Linde, 1990). 



2.2 Quantum Fluctuations 

Inflation predicts emergence of quantum fluctuations in early universe. As soon as the fluctua- 
tions emerge from a vacuum, the exponential expansion stretches the proper wavelength of the 
fluctuations out of the Hubble-horizon scale, H . After leaving the horizon, the fluctuation am- 
plitude does not change in time; on the contrary, it stays constant in time with characteristic r.m.s. 
amplitude, |0| rms ~ H/(2tt). 

After inflation, as the universe decelerates, the fluctuations reenter the Hubble horizon, seeding 
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matter and radiation fluctuations in the universe. Figure 2.2 summarizes the evolution of character- 
istic length scales: the Hubble-horizon scale (H" 1 ), the COBE DMR-scale fluctuation wavelength, 
and the galaxy-scale fluctuation wavelength. 

We estimate H -1 during inflation as follows. The scalar-field fluctuations produce CMB fluc- 
tuations of order H/m p \. Using the DMR measurement flSmoot et al., 1992 ), AT/T ~ 10 -5 , we 



obtain H ~ 10 _5 m p i, or H^ 1 ~ 10 5 m~^ ~ 10~ 28 cm. Since H -1 stays nearly constant in time 
during inflation, this value represents the horizon scale throughout inflation approximately. After 
inflation, H^ 1 grows as H^ 1 (a) oc a 2 in the radiation era, and oc a 3 / 2 in the matter era. 

DMR probes a present-day fluctuation wavelength on the order of 3 Gpc ~ 10 28 cm. By 
comparing the reheating temperature, ~ 10 27-29 K, with the present-day CMB temperature, 2.73 K, 
one finds that the universe has expanded by a factor of ao/a r h ~ 3 x 10 26-28 since the reheating; 
thus, the DMR scale corresponds to a proper wavelength of 0.3 — 30 cm at the reheating epoch 
(the number could be more uncertain). Here, ao is the present-day scale factor, while a r h is the 
reheating epoch. 

The galaxy-scale fluctuations have the linear comoving wavelength on the order of 1 Mpc ~ 3 x 
10 cm. The galaxy-scale fluctuations have left the horizon later than the DMR-scale fluctuations: 
«rh/ a gai ~ 10 25 , while a r h/admr ~ 10 28 . Here, a ga i and ad mr are the scale factors at which the 
galaxy- and DMR-scale fluctuations leave the horizon, respectively. 

These ratios are often calculated with e-folding numbers, N = ln(a r h/a). For the galaxy- 
and DMR-scale fluctuations, we have N ga x = ln(a r h/a ga i) ~ 58, and iVdmr = l n («rh/«dmr) ~ 64. 
Moreover, using equation ( |2.1[ ), we obtain t ga i — td mi = H^ 1 (N^ mi — N ga \) ~ 10~ 38 s; thus, inflation 
generates the fluctuations on the DMR scales down to the galaxy scales almost instantaneously. 

2.2.1 Quantization in de Sitter spacetime 

A basic idea behind quantum-fluctuation generation in inflation is well described by the second 
quantization of a massive-free scalar field in unperturbed de Sitter spacetime, for which a{t) = 
aQe H ^~ to ' with H independent of t. In this system, the problem is exactly solvable, and finite mass 
captures an essential point of generating a "tilted" fluctuation spectrum, as we will show in the 
next subsection. In this subsection, we describe the quantization procedure in de Sitter spacetime, 
following Birrell and Davies (1982| ). 



Quantization in curved spacetime is generally complicated, as there is no unique vacuum state to 
define a ground state of quanta, even for inertial observers who detect no particles in the Minkowski 
vacuum (the vacuum state for the quantum field theory in Minkowski spacetime). Fortunately, in 
the spatially-flat Robertson- Walker metric, there is a prescription for quantization, largely because 
of the metric being conformal to the Minkowski metric, g^ u = a?"q^ u , or more specifically 

ds 2 = -dt 2 + a 2 {t)5 ij dx i dx j = a 2 (r) (-dr 2 + S ij dx i dx j ^j , (2.3) 

where r = J a~ 1 (t')dt' is called the conformal time. In this metric, there is a reasonable vacuum 
state for particles whose comoving frequencies, u>k, are higher than the conformal expansion rate, 
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Figure 2.2: Physical Length Scales in Cosmic Evolution 
Representative physical length-scale evolution in the cosmic history. The thick lines draw the 
evolution of fluctuation wavelengths, 2ira/k, while the thin line draws the evolution of the 
Hubble-horizon scale, H^ 1 (a), where a is the expansion scale factor. H^ 1 (a) stays constant 
during inflation, grows as a 2 in the radiation era, and grows as a 3 / 2 in the matter era. There 
are characteristic scale factors: ao is the present-day, a cq is the matter-radiation equality, and 
a r h is the reheating epoch. The DMR-scale fluctuations leave the horizon at ad mr , while the 
galaxy-scale fluctuations leave at a ga i . These scale factors are related to each other roughly 
as follows: a /a cq <~ 10 4 , a cq /a r h ~ 10 24 , a r h/a ga i <~ 10 25 (e 58 ), and o r h/fldmr ~ 10 28 (e 64 ). 
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dlna/dr. We should differ the conformal expansion rate, d\na/dr, from the expansion rate, H 
dlna/dt. They are related through dlna/dr = aH. In de Sitter spacetime, we obtain 



T = T + 



1 - exp [-H(t - t )\ 
a H 



T + 



1 



1 



ciqH J cl{t)H 



(2.4) 



For simplicity, we set the first term vanishing, so that r = — [a(r)H]~ ; thus, r lies in — oo < r < 
for — oo <t<oo(0<a< oo). From now on, we will use dots for conformal-time derivatives: 
x = dx/dr. 

We start quantizing a scalar field, c/>(x, r), by expanding it into creation, at, and annihilation, 



a k) operators, which satisfy the commutation relation 

d 3 k 



^(3) (k- k'). We have 



(x,r) 



a k ^(T-)e* k ' x + a^(r)e 



ikx 



(2tt) 3 / 2 

The canonical commutation relation between cj> and the conjugate momentum, = a 2 (r)<p: 



[0(x,r),7r 9 i(x / ,r)] = a 2 (r) 



x',r 



^ 3 )(x-x'), 



(2.5) 



(2.6) 



gives a normalization condition on ^(r), a 2 ((fkfik ~ Pk^k) = 

The normalization condition motivates our using a new mode function given by Xk = o^fe? which 
satisfies a new normalization condition, XfcXfc — X%Xk = i- If Xk has a positive frequency mode with 
respect to the conformal timelike Killing vector (d/dr for our metric), i.e., Xk = —i^kXki then the 
condition gives Xfc( r ) = (2tiJ] t )~ 1 ' 2 e~ iu3kT , a ground state in the Minkowski vacuum. Since Xk gives 
the closest analogy to the Minkowski vacuum state, we will use Xk more frequently than ipk. 



The Klein-Gordon equation for a massive-free scalar field, g a ^ 
motion for Xfc( r )> 

Xk(r) + [/c 2 + m 2 (r)] X k(r) = 0, 
where m 2 (r) is the time-dependent effective mass, 



>;/3 



2^,2 



gives equation of 



(2.7) 



m 2 (r) 



m 



2F> 2 (r) = mV(r) - 



(2. 



We thus find that the Hubble parameter effectively reduces m 2 by 2H 2 . This time-dependent, 
negative contribution to m^r) is the effect of de Sitter spacetime, which is not Minkowski but 
curved. 

Fortunately, there is an exact solution to the Klein-Gordon equation ([2.7|): 



Xk(r) 



dH® {-kr) + C 2 H, 



(2) 



-kr) 



(2.J 



where C\ and C2 are integration constants, and v 2 = 9/4 — m 2 /H 2 . H„ (x) is a Hankel function 

/9 s ) r (-\\ 1 * 

of the first kind; H v (x) = H„ (x) . We have negative sign in front of r to recall that r lies in 
-00 < r < 0. 
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How do we determine the integration constants, C\ and C2? In other words, how do we normalize 
our mode function properly? Since we know how to quantize a scalar field in the Minkowski vacuum, 
we should find a mode function that matches the Minkowski positive frequency mode, Xfc( r ) = 
(2oo k )~ 1/2 e~ iu}kT ; however, we cannot find an unique positive frequency mode valid throughout 
inflation, as the time-dependent spacetime creates particles. 

Instead, we define a vacuum state in the in state, the remote past, r — ► —00. Using an 
asymptotic form of the Hankel function, 



cxp 



7T 



(2.10) 



we obtain Xk in the in state (r 



-00) from equation ( gjj ) 
~2 



Xk{T 



-OO 



ie ikT + C 2 e ikT 



(2.11) 



Here, we have neglected the contribution from m 2 /H 2 compared with — kr in the exponent. The 
second term has a negative frequency, so that C2 = 0. The first term with C\ = gives 
Xk( T ) = {2k)~ 1 l 2 e~ lkl ~ ', the Minkowski positive frequency mode with uj^ = k, and thus in the in 
state, the solution describes a ground state of a massless field in the Minkowski vacuum. 
Using the solution for x/c( T )> we obtain a solution for </?fc(r), 



and </>(x, r) becomes 
</>(x,r) 



<Pk{r) 



d 3 k 



2a(r) 



(i) 



(2.12) 



V 7 — 7TT 

2a(r) 7 (2tt) 3 / 2 



ikx 



(2.13) 



Since all the k modes in the integral are independent of each other, the nearly infinite sum of those 
modes makes <j> obey Gaussian statistics almost exactly, because of the central limit theorem; thus, 
the two-point statistics specify all the statistical properties of <f>. This is a generic property of the 
ground-state quantum fluctuations. 

The annihilation operator, ak, annihilates the vacuum state defined in the in state: Sk| 0i n ) = 0. 
In this vacuum, we calculate amplitude of ground-state (ft fluctuations as 



i: 



<^(x, t)0(x,t) 



k 2 dk 2 
■^rlWkWI 



8vra 2 (r) 



k 2 dk 



(2.14) 



Since we probe a limited range of k observationally, we also use the fluctuation spectrum in a log- 
arithmic k range, A 2 (fc), which represents variance of fluctuations at a given comoving wavelength 

27T/C- 1 , 



k 3 
2^2 



A K k ) = 7^ \<Pk(r) 



-kh 



H^i-kr) 



8vra 2 (r) 



E-(-kT) 3 \HU(-kr) 



(2.15) 
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where we have used a(r) = —(Ht)^ 1 in the last equality. Let us recall that v 2 = 9/4 — m 2 /H 2 , 
and r lies in — oo < r < 0. 

Here, we have the quantization of 4> completed, and formally calculated the fluctuation spec- 
trum. These results are exact, and valid on all scales. In the next subsection, we study the solution 
on super-horizon scales, where inflation produces observationally relevant fluctuations. 



2.2.2 Scale-invariant fluctuations on super-horizon scales 

Equation ( |2.13| ) describes a quantum massive-free scalar field, (f), in the unperturbed de Sitter 
spacetime on all scales. As the universe expands exponentially, —kr = k/(aH) quickly becomes 



very small; the mode leaves the Hubble-horizon scale, H . Figure ~L2 shows that the fluctuations 
on the observationally relevant scales should have left the horizon during inflation. Hence, the 
behavior of <j> on super-horizon scales is practically important. In this subsection, we study the (j) 
fluctuation spectrum on super-horizon scales. 



In equation ( 2.15 ), using an asymptotic form of the Hankel function, 



vr \2 

we obtain the <ft fluctuation spectrum on super-horizon scales, 



, (2.16) 



3-2u 

(2.17) 



Gtt, 


J 2 2 2,-3 


r m - 


2 


t k \ 






|r(3/2)_ 




\aH) 



or the spectral index, 



where v 2 = 9/4 - m 2 /H 2 . We have used r = -(aH)' 1 and vr = T 2 (l/2) = 4r 2 (3/2). 

One finds that v = 3/2 is a special point, for which the spectrum is independent of k, i.e., 
scale invariant, A 2 (/c) = H 2 /{2ir) 2 . This happens when we assume m 2 /H 2 <C 1, and thus v = 
3/2 - m 2 /(3H 2 ) + 0(m 4 /H 4 ), which gives 

( H\ 2 ( k \ 2m2 /(3H 2 ) 

dlnAj _ 2m 2 
^lnT " MP' (2 ' 19) 
Since m 2 <C H 2 , the spectrum is almost scale invariant, giving the <j> fluctuations characteristic 
r.m.s. amplitude, |<^| rms = H/(2ir). Finite mass makes the spectrum slightly "blue", the power of 
k being positive. 

The above assumption, m 2 <C H 2 , offers long-lasting inflation that makes the observable uni- 
verse flat and homogeneous; otherwise, (j) rolls down to a potential minimum too quickly, terminating 
inflation too early. In the inflationary regime, the Friedmann equation gives 

to 2 3 rn 2 i „ 

pl ~ 10" 2 . (2.20) 



H 2 4vr 



2. 2. Q UANTUM FL UCTUATIONS 



23 



log[A 2 (k)] 



blue 



red 



log(k) 



Figure 2.3: Color of Spectrum 

A sketch of the fluctuation spectrum, A 2 (fc), which represents the fluctuation power in a 
logarithmic k range. The solid line plots a scale-invariant spectrum, the dotted line plots a 
"red" spectrum whose spectral index is negative, and the dashed line plots a "blue" spectrum 
whose spectral index is positive. 



Hence, for a massive-free scalar field to drive inflation, the mass cannot be comparable to the 
Hubble parameter, and the fluctuation spectrum is almost exactly scale invariant. 

The argument until now has assumed the exact de Sitter spacetime in which H is constant in 
time, and neglected perturbations in the metric. As a result, we have obtained a blue spectrum 
whose spectral index is < 2m 2 /3H 2 < 1. In realistic inflation models, however, none of the 
above assumptions apply: H decreases slowly in time, and the metric is perturbed. In the next 
subsection, we will show that both the effects give a tilted "red" spectrum, for which the power of 
k is negative of order — 10 -2 . Figure 2.2 sketches what blue, scale- invariant, and red spectra look 
like. 



2.2.3 Tilted "red" spectrum 

If H decreases in time, and the metric is linearly perturbed, a scalar field acquires a negative 
effective mass-squared, Am 2 < 0. It also modifies the mass of Xki m x( T ) (Eq.(|2.8[)), as 



m 2 (r) — > (m 2 + Am 2 ) cl 2 (t) - (2.21) 



We thus expect that Am 2 modifies the fluctuation spectral index, for the mass determines the 
spectral index through equation (|2. 19 ) . This parameterization, Am 2 , may be useful to understand 
what physical effect is responsible for the spectral index. 

How large is Am 2 ? The derivation of Am 2 is quite involved, but Mukhanov et al. (1992) show 
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that the exact form of m 2 (t) 



H d 2 (0/H) _ f d 2 V 
6 dr 2 ~ I d6 2 



+ 



'■aJ a(T) 



2 2, 

m cfr a ( 



(2.22) 



Here, since -ff and <p change in time t only very slowly, we have retained the terms on the order of 
—dH/dt = 4:irm~ 2 (dtp / dt) 2 or larger. This approximation is called the slow-roll approximation, and 
we evaluate the approximation explicitly in appendix [A|. For a massive-free field, V{4>) = m 2 (j) 2 /2, 
by comparing equation ( 2. 22] ) with (2.21), we find Am 2 = 9dH/dt < 0. One can show that the 3 
of the 9 comes from the effect of H changing in time, and the 6 comes from the effect of the metric 
perturbations. 

By using the first-order slow-roll approximation, we obtain a relation between dH/dt and the 
potential slope, dV/d<j), dH/dt rj — -^{dV/dcj)) 2 = — \m 2 . It thus follows that the total effective 
mass-squared becomes negative, m 2 s = m 2 + Am 2 
a negative spectral index, 



-2m < 0, and we have a "red" spectrum or 



dink 



2<r 
3H 2 



Am 2 
3IP 



< 0. 



(2.23) 



Let us summarize what has made the spectral index negative. The spectral index is determined 
by the mass, or the effective mass, of a scalar field. In inflation, the intrinsic mass, m 2 , is over- 
compensated by the induced mass from gravitational effects, Am 2 , which is negative: m 2 -\- Am 2 < 
0. As a result, the fluctuation spectrum becomes "red". Actually, any power-law potential of the 
form V((f>) = A(fi n with A positive gives a negative m 2 s through m 2 s = — (1 + n / 2)nA(f> n ~ 2 . Here, 
we have used dH/dt w —^{dV/d<j)) 2 = -\n 2 A^ n ~ 2 . 

For a generic scalar field with an arbitrary potential, we find 



d 2 V 



9- 



dH 
~~dt 



6- 



JH 

IF 



3H 



dP^/dt 7 



d(f>/dt 



and the spectral index 



din A 2 



% _ dH/dt d 2 (f>/dt 2 



dink 



H 2 



H{d(j)/dt) ' 



(2.24) 



(2.25) 



which agrees with Liddle and Lyth (1992 ). 

CMB experiments have shown that a scale-invariant fluctuation spectrum fits the data well (dc 



Bernadis et al., 2000 ; Hanany et al., 2000 ). Combining all the CMB experiments to date, Wang 
et al. (2001| ) show that a slightly red spectrum, the power of k being ~ —0.07, fits the data even 
better, while the error of the fit is still of order 0.1. Accurate measurement of the spectral index 



constrains the shape of V{4>) through equation (2.25), as V(<fi) determines the time variation of H 
and 4>; thus, it potentially discriminates between different inflation models. 

Although the CMB experiments almost exclude possibility of the massive-free scalar field with 
sizable m 2 /H 2 ~ 1 dominating the matter and radiation fluctuations in the universe (in terms of the 
spectrum index), the field may produce some of the fluctuations that are non-Gaussian. In this case, 
4> does not drive inflation, so that m 2 can be comparable to H 2 ; the field rolls down to a potential 
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minimum quickly, and oscillates about the minimum, producing non-Gaussian isocurvature density 
fluctuations, bp oc 4> 2 flLinde and Mukhanov, 1997 ). Even if the isocurvature fluctuations are sub- 



dominant in the universe, they could produce non-Gaussian temperature fluctuations in CMB. 
Measuring non-Gaussianity in CMB thus potentially probes particle physics in inflation. 

2.2.4 Emergence of classical fluctuations 

Until now, we have considered generation of quantum fluctuations in inflation, and derived a fluctu- 
ation spectrum on super-horizon scales. We then expect that the fluctuations seed observed CMB 
anisotropies and large-scale structures in the universe; but, how? Inflation generates quantum fluc- 
tuations, not classical fluctuations. How can the quantum fluctuations make classical objects like 
galaxies seen today? 

Generation of classical fluctuations, or quantum-to-classical transition of fluctuations, during 
inflation, has been in debate (Calzatta and Hu, 1995; Matacz, 1997a| ,|b|; |Kicfcr et al., 199*8 ). In this 



subsection, we describe a possible mechanism to produce classical fluctuations in inflation. Our 



approach is partly close to Kiefer et al. (199 



A basic idea behind our classical-to-quantum transition mechanism is to approximate a field, (f> 
(Eq,([2.5|)), with the sum of the long- wavelength (super-horizon) modes and the short-wavelength 
(sub-horizon) modes. This mode separation is unique, as the comoving Hubble-horizon scale, 
(aH)' 1 , has been a characteristic length scale in the solutions for mode functions, <Pk(i~) (Eq. fl2.12j )). 



In the super-horizon limit, k <C aH, equations ( |2.12[) and ( |2.16[ ) give <pl<$:aH = ~<fik<s:aH- Using 
this in equation ( [2.5D , we obtain 



ik-x 



0(x,t) ~ / d 3 k (a k - a^ k ) v9 fc « a ^(T)e i 
Jk<aH v 7 

+a - V)/ -*(y>~<» + ^™)- iw 

Jk>aH (2Try/ 2 v2k v 7 

We name the first term (pch, the second term (J)qm- The corresponding conjugate momenta are 
ti"cl = cl 2 4>cl and 7tqm = a 2 <pQM, respectively. The first term, (ftch, becomes just the Fourier 
transform, while the second term, </>qmj becomes an ordinary ground state in the Minkowski vacuum 
except for a _1 (r) in the front. To see what happens to these terms in the context of the quantum 
field theory, we calculate the canonical commutation relations of (frcL and <^qm- 
For 4>qm, we find 

lk>aH faY 

for |x — x'| < (aH)' 1 ; thus, i^qm is a quantum field inside the horizon, and a<^>QM is identical to a 
quantum field in the Minkowski vacuum. For (ficL, we nn d the commutation relation vanishing, 



[^QM(x,r),7r Q M(x',r)] =i / __ e *-(*-*') = (x - x') (2.27) 

Jk>aH (27TP 



[0CL(x,r),7r C L(x / ,r)] = 0; (2.28) 

thus, 4>cl is no longer a quantum field, but a classical field. In other words, once smoothing out 
the fluctuations inside the horizon, we are left with the classical fluctuations. Notice that there has 
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been no explicit decoherence mechanism in the system. The exact solution to the Klein-Gordon 
equation in the exact de Sitter spacetime naturally yields the classical fluctuations on super-horizon 
scales. 

Without smoothing out the sub-horizon fluctuations, however, eft remains quantum. The com- 
mutation relation of (j) (Eq. (|2.6| )) is exact as long as all the k modes equally contribute to <f>. If 
nothing happens to the sub-horizon fluctuations, then the quantum, not classical, fluctuations will 
reenter the horizon after inflation. Yet, practically speaking, this argument may not be relevant 
for our actual observations because of the following reason. Consider a fluctuation wavelength just 
leaving the horizon at the end of inflation, at which the proper horizon size is H" 1 ~ 10 -28 cm. 
As the universe has expanded by a factor of order 10 27 since the end of inflation, we may find 
the fluctuation wavelength to be ~ 0.1 cm in the present universe. Now we ask: "does this-size 
fluctuation affects classicality of galaxy-scale fluctuations?" We may answer "no", if we assume 
that these substantially different-scale fluctuations have undergone different physics; if so, quantum 
coherence should have disappeared. 

Even right after inflation, the super-horizon modes and the sub-horizon modes have undergone 
different physics. At the end of inflation, the reheating begins, and scalar fields decay into particles 
and radiation, thermalizing the universe. By causality, the thermalization process occurs only 
inside the horizon; thus, the reheating affects the sub-horizon fluctuations differently from the 
super-horizon fluctuations, breaking quantum nature of <p. In other words, in equation ( 2.26 ), the 
second term may have disappeared during the reheating, while the first term may remain unaffected. 
As a result, the reheating effectively smoothes out the sub-horizon scale fluctuations, and makes a 
quantum-to-classical transition possible. 



2.3 Linear Perturbation Theory in Inflation 

In the previous section, we have followed generation of scalar-field fluctuations in the unperturbed de 
Sitter spacetime, i.e., no perturbations in the metric. Scalar-field fluctuations, however, perturb the 
stress-energy tensor, and produce metric perturbations. Since the metric perturbations regulate the 
matter and radiation fluctuations that we observe today, we must include the metric perturbations 
in the analysis, and follow the evolution. In this section, we explore the linear perturbation theory 
in inflation, which includes perturbations to the metric and a scalar field. Our notation follows 
Bardeen (1980| ). 



We use a linearly perturbed conformal Robertson- Walker metric of the form, 

ds 2 = a 2 (r) {-(1 + 2AQ)dr 2 - 2BQ i drdx i + [(1 + 2H L Q) 5 tj + 2H T Q ij ] dx i dx j \ . (2.29) 

Here, all the metric perturbations, A, B, H^, and H^, are <C 1, and functions of r. The spatial 
coordinate dependence of the perturbations is described by the scalar harmonic eigenfunctions, Q, 
Qi, and Qij, that satisfy S^Q^j = —k 2 Q, Qi = —k~ 1 Q i i, and Qij = k~ 2 Q^j + ^5{jQ. Note that 
Qij is traceless: 8 ll Qij = 0. Kodama and Sasaki (1984 ) use different symbols, Y, Yi, and Yy, for 
Q, Qi, and Qij, respectively. 
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The four metric-perturbation variables are not entirely free, but some of which should be fixed 
to fix our coordinate system before we analyze the perturbations. The choice of coordinate system 
is often called the choice of gauge, or the gauge transformation; we will describe it later. 

2.3.1 Fluid representation of scalar field 

The metric perturbations enter into the stress-energy tensor perturbations, 5Tjf. We expand a 
scalar field into its homogeneous mean field, <j>(r), and fluctuations about the mean, d<j}(r)Q{x). 
The energy density and pressure fluctuations are given by 

Sp</,Q = -<5T ° = [a- 2 (J)6<j> - A/> 2 ) + V >4> 5^] Q, (2.30) 

5T k r / \ i 

SP4>Q = "jp = [ a ~ 2 - M 2 ) - Q. (2.31) 

The energy flux, Zf, gives the velocity field, v^Qi, 

(P4> + P4>) K - B )Qi = T i = [^ k5 ^j Qi- ( 2 - 32 ) 

Using p ( f,+P(f > = a~ 2 (p 2 , we obtain v^ — B = k(fi~ 1 54>; thus, 5<fi is directly responsible for the fluid's 
peculiar motion. The anisotropic stress, — p^Sj, is a second-order perturbation variable for a 
scalar field, being negligible. 

When we choose our coordinate system so as B = (a fluid element is comoving with the origin 
of the spatial coordinate), we have 5<p vanishing, 5(f) = 0. This coordinate is called the comoving 
gauge, and we write the scalar-field fluctuations in this gauge as 5(j) com = 0. 

Since we have only one degree of freedom, a scalar field, in the system, Spj,, Spj,, and vj, are 
not independent of each other. Nevertheless, this fluid representation is useful, as the cosmological 
linear perturbation theory has been developed as the general relativistic fluid dynamics. We can 
plague these fluid variables into the well-established general relativistic fluid equations, and see 
what happens to the metric perturbations. While we do not use those fluid equations explicitly in 
the following, but solve equation of motion for a scalar field (Klein-Gordon equation) directly, the 
fluid equations give the same answer. 

2.3.2 Gauge-invariant perturbations 

In the previous subsection, we have seen that scalar-field fluctuations vanish in the comoving gauge 
in which B = v^; thus, a choice of gauge defines perturbations. For the scalar-type perturbations 
that we are considering, the gauge transformation is 

r — ► t' = t + T(t)Q(x), (2.33) 
x l — ► x' 1 = x i + L(r)Q,(x), (2.34) 

where T and L are <C 1. Accordingly, scalar-field fluctuations, 5(f), transform as 

8<j>{r) — > 54>(t') = 54>{t) - 4>{t)T(t). (2.35) 
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Hence, if we choose T = <fi~ 1 8(f), then we obtain 8(f> = 0. This choice of T defines the comoving 
gauge, S(j) com = 0. In this way, we find different values for the perturbation variables in different 
gauges. 

So, what gauge should we use? Unfortunately, there is no answer to the question: "what gauge 
should we use?" . Although there is no best gauge in the world, depending on a problem that we 
intend to solve, we may find that one gauge is more useful than the other, or vice versa. As long 
as we fix the gauge uniquely, and understand what gauge we are working on clearly, no problems 
occur. 

In practice, however, problems occur when one author understands its own gauge, but does not 
understand the other author's gauge. Since there is no best gauge in the world, different authors 
may use different gauges, and may disagree with each other because of their misunderstanding of 
the gauges. In other words, one author's calculation on amplitude of 8(f> may disagree with the 
other's calculation, if they are using different gauges. The author using the comoving gauge sees 
8(f) = 0, but others may see 8(f) ^ 0. 

One way to overcome this undesirable property is to make perturbation variables invariant 
under the gauge transformation, and let them represent gauge-invariant perturbations. As an 
example, consider a new perturbation variable ( pVIukhanov et al., 1992| ), 

u = 8(P-^(h l + ijT T ) . (2.36) 



One can prove this variable gauge invariant, u = u, using equation (|2.35|) and 



H L + = H L + ^H T - aHT. (2.37) 

Actually, Hl + represents perturbations in the intrinsic spatial curvature, 7Z, as it is the scalar 
potential of the 3-dimension Ricci scalar: 8^R = a~ 2 k 2 TZQ, where 1Z = Hi, + \Hi. While u 
reduces to 8(f) in the spatially flat gauge (TZ = 0), or to — {(f)/aH)lZ in the comoving gauge {8(f) = 0), 
its value is invariant under any gauge transformation. Any authors should agree upon the value of 
u. 

For the physical interpretation of u, we may name u "scalar-field fluctuations in the spatially 
flat gauge" or "intrinsic spatial curvature perturbations in the comoving gauge". Either name 
describes the physical meaning of u correctly. The physical meaning of u depends upon what gauge 
we are using; however, the most important point is that the value of u is independent of a gauge 
choice. In this sense, u can be a "common language" among different authors. 

Bardeen et al. (198f) use a similar gauge-invariant variable to u, 

(=-^ u = 1Z-^8(f), (2.38) 



that reduces to TZ in the comoving gauge, or to —{aH / <f))8(f) in the spatially flat gauge. This variable 
helps our perturbation analysis not only because of being gauge invariant, but also being conserved 
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on super-horizon scales throughout the cosmic evolution. We will show this property in the next 
subsection. 

Using gauge invariance of u or £, we obtain a relation between 5(f> in the spatially flat gauge, 
<5</>fl at , and TZ in the comoving gauge, TZ com , as 

clH 

^com = p^flat- (2-39) 

It is derived from « com = tifl a t, or £ com = Cflat- As we have seen in the previous section, 5(f) obeys 
Gaussian statistics to very good accuracy because of the central limit theorem, that is, 5(f) is the 
sum of the nearly infinite number of independent modes (Eq.( |2.13| )). Since TZ com is linearly related 
to <5</>flat) ^com also obeys Gaussian statistics in the linear order; however, as we will show in the 
next section, non-linear correction to this linear relation makes lZ CO m weakly non-Gaussian. 

The spatial curvature perturbation, TZ, is more relevant for the structure formation in the 
universe than the scalar-field fluctuation, 5(f), itself, as TZ regulates the matter density and velocity 
perturbations through the Poisson equation. Actually, TZ reduces to the Newtonian potential inside 
the horizon. Since the quantum fluctuations generate 5(f), we expect it to generate TZ through 
K com = —(aH/<ft)5(f)fi a t. This is naively true, but may sound tricky. In the next subsection, we will 
show how to calculate TZ generated in inflation more rigorously. 



2.3.3 Generation of spatial curvature perturbations 

To calculate the intrinsic spatial curvature perturbation, TZ = + ^Ht, that is generated in 
inflation, we need to track its evolution equation, and figure out how it is related to 5(f). We 
will show in this subsection that TZ is actually more than related to 5(f); it is almost equivalent to 
5(f>. We can track the evolution of 5(f) and TZ simultaneously, using the gauge-invariant variable, u 
(Eq.Q). 

Mukhanov et al. (1992) show that au = a5(f> — (<j)/H)TZ obeys the same Klein-Gordon equation 
as we have used in the previous section (Eq.([T7|)), 

Xk+ [k 2 + m 2 x (T)] Xk = 0, (2.40) 

where Xk is the mode function that expands au (see Eq.( |2.5| )). It thus follows that a5(f> and (<f)/H)TZ 
obey the same equation, and our argument on the quantum-fluctuation generation during inflation 
in the previous section applies to u as well, u being quantum fluctuations means that it also obeys 
Gaussian statistics very well because of the central limit theorem (see Eq.( 2.13| ) and the text after 
equation) . 



We consider the Klein-Gordon equation for au on super-horizon scales. As equation (|2.5| ), we 
expand au into the mode functions, \k = acpk, where Xk and ipk are exactly the same functions 
that we have used in the previous section. 

We give a slightly different expression for the solution, emphasizing its time dependence on 
super-horizon scales. Taking the long- wavelength limit, k 2 <C m 2 {r), and using the exact form of 
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m^ir) (Eq. Q2.22| )) ( [Mukhanov et al., 1992| ), we obtain the Klein-Gordon equation on super-horizon 
scales, 

Hd 2 (d>/H) , A . 

There is an exact solution to this equation, 

H clH „ „ f H 2 , . . 

TXfc = = d + C 2 -j^dr, (2.42) 

(ft (ft J (ft 

where C\ and C2 are integration constants independent of r. The second term is a decaying mode as 
J (ft~ 2 H 2 dr = J a~ i (d(ft/dt)~ 2 H 2 dt, and thus (H/(ft)xk remains constant in time on super-horizon 
scales. This implies that £ = —{aH/<j))u also remains constant in time on super-horizon scales. 
Note that £ obeys Gaussian statistics in the linear order, as it is related to a Gaussian variable, it, 
linearly; however, as we will show in the next section, non-linear correction to this linear relation 
makes £ weakly non-Gaussian. This statement is equivalent to that we have made on TZ com . 

The solution obtained here for £ is valid throughout the cosmic history regardless of whether a 
scalar field, radiation, or matter dominates the universe; thus, once created and leaving the Hubble 
horizon during inflation, £ remains constant in time throughout the subsequent cosmic evolution 
until reentering the horizon. The amplitude of £, i.e., C±, is fixed by the quantum-fluctuation 
amplitude derived in the previous section (Eq.(|jJ|)), 



fc 3 , „ , 2 / aH\ 2 2 / aH 2 



)■- 


H 2 




2ir(d(j)/dt) 



(2.43) 



This is the spectrum of £, A 2 (k), on super-horizon scales. While we have neglected the k dependence 
of the spectrum here, the spectral index of £ is the same as of (f> (Eq.( |2.25| )), 

dlnA 2 _ idH /dt 2 d 2 Hdt 2 

dink H 2 H(d(p/dt)' v ' ; 

A 2 (k) gives the primordial curvature-perturbation spectrum. This is very important prediction 
of inflation, as it directly predicts the observables such as the CMB anisotropy spectrum and 
the matter fluctuation spectrum. Strictly speaking, A. 2 (k) reduces to the curvature-perturbation 
spectrum in the comoving gauge. 

To summarize, the quantum fluctuations generate the gauge-invariant perturbation, u, that 
reduces to either <5</>fl a t or (cp/aH)TZ COTn depending on which gauge we use, either the spatially 
flat gauge or the comoving gauge. Hence, 5^fl a t and ((ft / aH)TZ com are essentially equivalent to 
each other. The benefit of u is that it relates these two variables unambiguously, simplifying the 
transformation between Scft^ and TZ com . This is a virtue of the linear perturbation theory; we do 
not have this simplification when dealing with non-linear perturbations for which we have to find 
non- linear transformation between Scft^ and TZ com . The non-linear transformation actually makes 
K com weakly non-Gaussian, even if dcp^t 1S exactly Gaussian. We will see this in the next section. 

Here, we have the generation of the primordial spatial curvature perturbations completed. In 
the next subsection, we will derive the CMB anisotropy spectrum. 
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2.3.4 Generation of primary CMB anisotropy 

The metric perturbations perturb CMB, producing the CMB anisotropy on the sky. Among the 
metric perturbation variables, the curvature perturbations play a central role in producing the 
CMB anisotropy. 

As we have shown in the previous subsection, the gauge-invariant perturbation, £, does not 
change in time on super-horizon scales throughout the cosmic evolution regardless of whether a 
scalar field, radiation, or matter dominates the universe. The intrinsic spatial curvature perturba- 
tion, TZ, however, does change when equation of state of the universe, w = p/p, changes. Since £ 
remains constant, it may be useful to write the evolution of 1Z in terms of £ and w; however, 1Z is 
not gauge invariant itself, but £ is gauge invariant, so that the relation between 1Z and £ may look 
misleading. 

Bardeen (1980) has introduced another gauge-invariant variable, $ (or $h in the original nota- 
tion), which reduces to 1Z in the zero-shear gauge, or the Newtonian gauge, in which B = = H^. 
<I> is given by 

Here, the terms in the parenthesis represent the shear, or the anisotropic expansion rate, of the 
r = constant hyper surf aces. While <& represents the curvature perturbations in the zero-shear 
gauge, it also represents the shear in the spatially flat gauge in which 1Z = 0. Using <£, we may 
write £ as 

aH aH ( H T \ 

c=K _ T ^ = # __^__ij, (2 . 46) 

where the terms in the parenthesis represent the gauge-invariant fluid velocity. 

Why use <3?? We use <3? because it gives the closest analogy to the Newtonian potential, for $ 
reduces to 1Z in the zero-shear gauge (or the Newtonian gauge) in which the metric (Eq.( 2T29| )) 



becomes just like the Newtonian limit of the general relativity. It thus gives a natural connection 
to the ordinary Newtonian analysis. 

The gauge-invariant velocity term, v — k~ l Hi, differs £ from In other words, the velocity and 
<3? share the value of Q. Since a fraction of sharing depends upon equation of state of the universe, 
w = p/p, the velocity and $ change as w changes. £ is independent of w. 

The general relativistic cosmological linear perturbation theory gives the evolution of $ on 



super- horizon scales ( Kodama and Sasaki, 1984 ) 



$ = ^^C, (2-47) 
for adiabatic fluctuations, and hence $ = |C in the radiation era (w = 1/3), and $ = |C i n t ne 



matter era (w = 0). $ then perturbs CMB through the so-called (static) Sachs- Wolfe effect (Sachs 
and Wolfe, 1967j ). 



The Sachs-Wolfe effect predicts that CMB that resides in a $ potential well initially has an 
initial adiabatic temperature fluctuation of AT/T = [2/3(1 + if)]^ and it further receives an 
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Figure 2.4: The Static Sachs-Wolfe Effect 

The static Sachs-Wolfe effect predicts that CMB that resides in a $ potential well has 
an initial adiabatic temperature fluctuation of AT/T = |<E> in the matter era. It further 
receives an additional fluctuation of — $ when climbing up the potential at the decoupling 
epoch. In total, we observe AT/T = 



additional fluctuation of — <]? when climbing up the potential at the decoupling epoch. In total, the 
CMB temperature fluctuations that we observe today amount to 



AT 
~T~ 



3(1 + w) 



$ _ $ 



l + 3w 

3 + 3u> 



l + 3w 
5 + 3w 



C- 



(2.48) 



Figure 2.4 sketches the static Sachs-Wolfe effect. 



For isocurvature fluctuations, initial temperature fluctuations in a potential well are given by 
— $ in both the radiation era and the matter era; thus, total temperature fluctuations amount to 
AT/T = — <£ — $ = — 2<3?. By definition of the isocurvature fluctuations, <3? is initially zero, but there 
exist non-vanishing initial entropy fluctuations. As the universe evolves, the entropy fluctuations 
create and hence the temperature fluctuations. 

At the decoupling epoch, the universe has already been in the matter era in which w = 0, 



so that we observe adiabatic temperature fluctuations of AT/T 
fluctuation spectrum of the Sachs- Wolfe effect, Ag W (fc), is 



A 2 



H 2 



10vr(#/dt) 



= -k, and the CMB 



(2.49) 



where H is the Hubble parameter during inflation. While we have not shown the k dependence of 
the spectrum here, the spectral index is given by equation (2.44). By projecting the 3-dimension 
CMB fluctuation spectrum, A| W (A;), on the sky, we obtain the angular power spectrum, C\ ( |Bond 
and Efstathiou, 1987|) , 



47T 



c 2 



sw J 



I - n)/2] r [I + (n - l)/2] 
r [(n + 3)/2] r [I + (5 - n)/2] : 



(2.50) 



where To and Td oc denote the present day and the decoupling epoch, respectively, and n = 1 + 
[din A 2 (k)/dln k] is a spectral index which is conventionally used in the literature. If the spectrum 
is exactly scale invariant, n = 1, then we obtain Cf w = [1(1 + l)]" 1 6Cf w . 
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Equation ( |2,50D provides a simple yet good fit to the CMB power spectrum measured by COBE 
DMR. The spectrum comprises two parameters, C| w and n. [Bennett et al. (1996| ) find n = 1.2±0.3 



and C| w = 7.9t?;2 x 10~ u . When fixing n = 1, they find C^' 11 " 1 = (9.3 ± 0.8) x HT 11 . The 
measured C\ is thus consistent with CMB being scale invariant, supporting inflation. Moreover, 

1 /2 

it implies that H 2 / (dtp / dt) ~ ^C| w ^ ~ 1CP 5 , constraining amplitude of the Hubble parameter 
during inflation. 

On the angular scales smaller than the DMR angular scales, the Sachs-Wolfe approximation 
breaks down, and the acoustic physics in the photon-baryon fluid system modifies the primordial 
radiation spectrum ( Peebles and Yu, 1970 ; Bond and Efstathiou, 198^ ). To calculate the modifica- 



tion, we have to solve the Boltzmann photon transfer equation together with the Einstein equations 
The modification is often described by the radiation transfer function, gTi(k), which can be cal- 
culated numerically with the Boltzmann code such as CMB FAST (Seljak and Zaldarriaga, 1996) 
Using gTi{k), we write the CMB power spectrum, C/, as 



coo AU 

d = 4tt / - Al(k)g 2 Tl (k). (2.51) 



o 



k 



Note that for the static Sachs- Wolfe effect, adiabatic fluctuations give gTi(k) = [k(jQ — Td ec )], 
while isocurvature fluctuations give g^i{k) = — 2ji [k(To — Td ec )]. Here, we have used A|(A;) rather 
than A| w (/c) or A 2 (k), following the literature. The literature often uses the $ power spectrum, 
P${k), to replace A\{k); the relation is A|(A;) = (2tt 2 )^ 1 k 3 P^(k) . A|,(A:) is called the dimensionless 
power spectrum. 

If $ were exactly Gaussian, then Ci would specify all the statistical properties of $, which are 
equivalent to those of £• Since C is related to a Gaussian variable, u, through ( = —(aH/(p)u, in the 
linear order £ also obeys Gaussian statistics; however, the relation between £ and u becomes non- 
linear when we take into account non-linear perturbations. As a result, and hence becomes 
non-Gaussian even if u is exactly Gaussian, yielding non-Gaussian CMB anisotropies. In the next 
section, we will analyze non-linear perturbations in inflation. 

Using the second-order gravitational perturbation theory, Pyne and Carroll (1996j ) derive the 



second-order correction to the relation between AT and (Eq.( |2.48| )). It gives AT/T = + 
0(l)& 2 ; thus, even if $ is Gaussian, AT becomes weakly non-Gaussian. 



2.4 Non-linear Perturbations in Inflation 

In the previous section, we have shown that the quantum fluctuations generate the gauge-invariant 
perturbation, u = 5(j) — (cft/aH)lZ, and u obeys Gaussian statistics very well because of the central 
limit theorem. Another gauge-invariant variable, £ = —(aH/(p)u = 7Z — (aH/(p)5(j), which remains 
constant in time outside the horizon, also obeys Gaussian statistics in the linear order, as it is 
related to a Gaussian variable, u, linearly. 

In the non-linear order, however, the situation may change. In the relation between £ and u, 
the factor in front of u, aH/<fi, is also a function of (f>, and it may produce additional fluctuations 
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like \d{aH / 4>) / d(j>\5<t>. Suppose that TZ in the comoving gauge (S(f) com = 0), 1Z com , is an arbitrary 
function of a scalar field: 1Z com = f((j>). Note that lZ com is equivalent to Q. By perturbing <p as 
4> = 4>o + 54>Ha.ti where <5</>fl a t is a scalar- field fluctuation in the spatially flat gauge (7£fl a t = 0), we 
have 

= /w,o) + (If) %at + 5 (0) ^ L + ° (*0 • (2 - 52) 

By comparing this equation with the linear-perturbation result (Eq. Q2.3S| )), 1Z com = —{aH / 4>)5(f)^ at , 
we find /(0q) = 0, df/dcj) = -(aH/j>), and 




(2.53) 



thus, even if <50fl at is exactly Gaussian, TZ com , and hence £, becomes weakly non-Gaussian because 
of 

^flat or ^ ne higher-order terms. While the treatment here may look rather crude, we will show 
in this section that the solution (Eq. (|2.53| )) actually satisfies a more proper treatment of non-linear 
perturbations in inflation. 

In the linear regime, we have the gauge-invariant perturbation variable that characterizes the 
curvature perturbations as well as the scalar-field fluctuations, and the single equation that describes 
the perturbation evolution on all scales. In the non-linear regime, however, we cannot make such 
great simplification. Since the Einstein equations are highly non-linear, fully analyzing non-linear 
problems is technically very difficult. Hence, we need a certain approximation. 



2.4.1 Gradient expansion of Einstein equations 

In inflation, there is an useful scheme of approximation, the so-called anti-Newtonian approximation 
(Tomita, 1975, 1982; Tomita and Deruelle, I99"2| ), or later called the long- wavelength approximation 
(Kodama and Hamazaki, 1998; Sasaki and Tanaka, 1998) or the gradient expansion method ( |Salopek 
and Bond, 1990| ; Salopek and Stewart, 1992; Nambu and Taruya, 1996, 1998| ). 

The approximation neglects higher-order spatial derivatives in the Einstein equations as well as 
in the equations of motion for matter fields, and is equivalent to taking a long-wavelength limit of 
the system. The equation system is further simplified if we set the shift vector zero in the metric. 
Once neglecting higher-order spatial derivatives and the shift vector, one finds that the shear decays 
away very rapidly. 

We use the metric of the form 



ds 2 



-N 2 dr 2 + ^g ij dx i dx j , 



(2.54) 



where iV is the Lapse function, and ^g%j describes the 3- metric. We have set the shift vector zero; 
it corresponds to B = in the linearized metric ( |2.29 ). We perturb N and ^ gij non-linearly. While 
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the full treatment of non-linear evolution of the system is highly complicated, by neglecting higher- 
order spatial derivatives, we reduce the Einstein equations to rather simplified forms flSalopek and 



Bond, 19901) : 



H 2 



a 13 a. 



8vrG 



Hj 



H + 3H Z + aVaij 



3 

-4vrG 
-8ttG 



2N 2 



N 



2N 2 



ft + V(</>) 



ft - v{4>) 



ma) 



0. 



(2.55) 

(2.56) 
(2.57) 
(2.58) 



H is the inhomogeneous Hubble parameter which defines the isotropic expansion rate of the 
r = constant hypersurfaces, 

' (2.59) 



1 

6iV 



or ^ (jij = 2NH^ gij. By introducing the inhomogeneous scale factor, a(x, r), given by H = 
a/ (Na), we can write the 3- metric as ^gij = a 2 (x, r)^ij(x), where 7ij(x) is a function of the spatial 
coordinate only. This form of ^gij satisfies equation ( |2.59| ). We have obtained the simplified form 
for the 3-metric because of setting the shift vector zero. 

a) is the shear which quantifies the anisotropic expansion rate, 

1 



H (s) 



9ij- 



(2.60) 



It follows from the traceless-part equation ( p. 58 ) and NH = a/a that the shear decays rapidly 
as the universe expands: a) oc a" 3 ; thus, we neglect the shear terms in the Einstein equations 
henceforth. 

By neglecting the shear term in the trace-part equation ( p. 57 ), and substituting the Fried- 
mann equation ( 2.55|) for H 2 , we obtain H = —&irGN~ l ft. Comparing this equation with the 
momentum constraint equation ( 2.56| ) without the shear term, = —A-KGN~ l cj)(j) t i, we find 
H(x,t) = H (0(x, r)), and the scalar-field momentum 



1 



f dH 



(2.61) 



N AttG V 

Substituting this for the kinetic term in the Friedmann equation ( [2.55 ), and neglecting the shear 
term, we finally obtain a closed evolution equation for H (0), 

1 (dH\- 
12vrG \~d^J 



H 2 (cP) 



+ — -V {(/>). 



(2.62) 



From this equation, H((f>) may be solved as H((p,I), where I is an integration constant which 
parameterizes the initial condition. This equation fully describes the evolution of the system in- 
cluding non-linear perturbations. Note that H(cft) depends upon the spatial coordinate through 



x,r . 
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A perturbation to H is given by 5H = (dH j d(f>) j 8<f> + (dH/ 81)^51. For the latter term, by 
differentiating equation ( 2.62| ) with respect to / for a fixed (ft, we find (dH/dl)^ oc a -3 ; thus, 
it decays very rapidly during inflation, giving 5H = (dH/d(f>) j 5(f). Hence, the comoving gauge, 
5(f) = 0, coincides with the constant Hubble parameter gauge, 5H = 0. Sasaki and Tanaka (199§| ) 
also observe this property from a different point of view, and find that this property holds for 
multiple scalar-field system as well. 



2.4.2 Generation of non-linear curvature perturbations 

Our goal in this subsection is to find a non-linear relation between TZ CO m and 5</>fl a t> following 
[Salopek and Bond (1990| ). In the absence of the shear, TZ obeys 7Z = 5 (NH) = 5 (dlna/dr), and 
hence 7Z is equivalent to fluctuations in the inhomogeneous scale factor, 5 In a. 

We calculate TZ com by perturbing In a non-linearly in the comoving gauge. Since (f> is homo- 
geneous in the comoving gauge, we choose <j) as a time coordinate: r = (f>. We find the Lapse 
function 

4vrG 

" = sm (2li3> 

for this time coordinate from equation ( |2.61| ) with setting 0=1. 

To calculate the scalar-field fluctuations in the spatially flat gauge, #0fl a t, we need to generate 
quantum fluctuations first; however, since we are solving the equation system on super-horizon 
scales only, we cannot calculate quantum fluctuations of 5(f) within the current framework. Instead, 
we assume that 5(f) on super-horizon scales is provided by the small-scale quantum fluctuations 
that are stretched out of the horizon by inflationary expansion. We use the linear perturbation 
theory to calculate the fluctuation amplitude at the horizon crossing, and provide 5(f) as initially 
linear, Gaussian fluctuations. We then calculate 50fl a t on the In a = constant hyper surf aces. While 
this treatment may sacrifice a virtue of the current framework which does not assume linearity 
of perturbations, non-linearity of the scalar-field fluctuations may also be incorporated into the 
analysis with the so-called stochastic inflation approach ( |5tarobinsky, 19*86 ; Balopek and Bond, 



1991). Using this, Gangui et al. (1994 ) show that the non-linearity makes 50fl a t weakly non- 
Gaussian. 

After all, our goal is to relate 5(f>(liia) to 5 In a((f>). In other words, we transform the perturba- 
tions on the In a = constant hypersurfaces to the ones on the (f) = constant hypersurfaces. We do 
this as follows ( Salopek and Bond, 1990| ), and figure (|2.5D shows the following process schematically. 



First, imagine — In a plane on which we transform the perturbations. We then mark a point 
on the plane with (</>o, In ao), and draw a short line from this point to {<pQ + 5(f), In ao) in parallel to 
the (p axis. This line represents 5(f>(x, In ag) — <K X > m a o) — 0o> i- e -> 4 1 perturbations on a In a = In ao 
hypersurface. Next, using evolution equation of In a, we evolve the line until it coincides with the 
(f) = (f)Q line, namely 

In a = In [a(x, (f) + 5(f))] — ► In [a(x, O )] . (2.64) 
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Figure 2.5: Hypersurface Transformation 
Hypersurface transformation of the scalar-field fluctuations in the spatially flat gauge, 
<50flat = 5<fi(lnao), into the spatial curvature perturbations in the comoving gauge, 
T^com = 51na((j>o). (1) To represent 5<f)(lnao), draw a short line from (0 o ,l na o) to 
(4>o + <50, In «o) in parallel to the <fr axis. (2) By evolving In a, evolve this line un- 
til it coincides with the <f> = line. (3) Measure the line length, <51na(0o)- 
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The evolution equation of In a is dlna/dr = NH. By integrating this equation, we obtain 

In [o(x, O )] - lnoo = (*° d<f> N ^ H ^\ ( 2 .65) 

J0O+5</>(x,lna o ) 

where the Lapse function, N, is given by equation ( |2.63| ) for <f> being the time coordinate. Finally, 
we measure the line length, 5 In [a(x, cfo)] = In [a(x, 4>o)] — lnao, to obtain 7£ C om- 

By noting <5(/>fl a t = S(f)(x, In ao), we find a non-linear relationship between 7£ com and <50flat> 



/ # 


'dlnH' 


-l 


ho 







(2.66) 

This non- linear relation should be compared with the crude estimate from the linear analysis, 
equation ( |2.53| ). By expanding the non- linear relation into Taylor series with respect to <5</>fl a t> we 
find that the crude estimate agrees with the non-linear result on the term-by-term basis. Note 
that a in equation ( 2.53|) should be replaced by N, as it has used the conformal time as the time 



coordinate, i.e., N = a. 

2.4.3 Generation of weakly non-Gaussian adiabatic fluctuations 

We have shown that non-linearity in inflation creates weakly non-Gaussian curvature perturbations 
outside the horizon from Gaussian quantum fluctuations inside the horizon. By expanding the non- 
linear relation between 7Z com and S(j>^ up to the second order, we obtain a non-linear curvature 
perturbation, lZ com = T^om + ^Sn) where 

(2.67) 
(2.68) 




Using equation ( |2.46 ), we obtain a non- linear Newtonian potential, = 3>l + &nl, where <3?l 



l^-com an d ^nl = f^-com i n the radiation era, and <I>l = f^-com an d $nl = f^com i n the matter 
era. Note that these formulae should be modified when we include non-linear effects of the stochastic 
inflation approach ( Gangui et al., 1994| ), which add extra terms to the formulae. 



In the following chapters, we will focus on measuring <I>nl using non-Gaussian CMB temperature 
fluctuations. We parameterize the amplitude of <J>nl with a non-linear coupling parameter, /nl, as 

$nl (x) = /nl [*£ (x) - (x))] . (2.69) 

How big is /nl? If a scalar field rolls down on a potential, V(4>), slowly, then din H/d(j) ~ 
^dlnV/dcj). We thus obtain 

/NL " ~2^G [-^- J ~ "487G I ~W~ J • (2 - 70) 
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As an example, consider a power-law potential, V(cp) oc 4> n . We find 

f - 5n m2pl in- 2 r9?n 
/NL -48^^~ 10 • (2 - 71) 

On the other hand, we will show in chapter || that |/nl| should be larger than of order unity 
to produce detectable non-Gaussian signals in CMB. Even for the fastest motion of 4> allowed in 
slowly-rolling single-field inflation, we find |/nl| < 5/2. Therefore, it is quite hard for this class 
of inflation models to produce detectable CMB non-Gaussianity. In other words, any detection of 
cosmological non-Gaussian signals in CMB constrains slow-roll inflation models very strongly, or 
may favor multiple scalar-field models or isocurvature fluctuations. 



Chapter 3 



Angular n-point Harmonic Spectrum 
on the Sky 

The angular n-point correlation function, 

(/(n 1 )/(n 2 ).../(n n )), (3.1) 
is a simple statistic characterizing a clustering pattern of fluctuations on the sky, /(n). Here, the 



bracket denotes the ensemble average, and figure 3.1 sketches the meaning. If the fluctuation is 
Gaussian, then the two-point correlation function specifies all the statistical properties of /(n), 
for the two-point correlation function is the only parameter in Gaussian distribution. If it is not 
Gaussian, then we need higher-order correlation functions to determine the statistical properties. 

Yet simple, one disadvantage of the angular correlation function is that data points of the 
correlation function at different angular scales are generally not independent of each other, but 
correlated: two-point correlation at 1 degree is correlated with that at 2 degrees, and so on. This 
property makes a detailed statistical analysis and interpretation of the data complicated. 

Hence, one finds it more convenient to expand /(n) into spherical harmonics, the orthonormal 
basis on the sphere, as 

oo I 

/(*) = £ E ai m Yi m (n), (3.2) 

l=Q m=-l 

and then to consider the angular n-point harmonic spectrum, {ai irril ai 2m2 . . . az nm „)- While a/ m for 
m ^ is complex, reality of /(n) gives a;_ m = a* m (— l) m , and thus the number of independent 
modes is not 41 + 1, but 21 + 1. 

Especially, the angular two-, three-, and four-point harmonic spectra are called the angular 
power spectrum, bispectrum, and trispectrum, respectively. For a Gaussian field, the angular spectra 
at different angular scales, or at different Z's, are uncorrelated. Even for a non-Gaussian field, they 
are reasonably uncorrelated as long as the non-Gaussianity is weak. Moreover, since the spherical 
harmonic is orthogonal for different Z's, it highlights characteristic structures on the sky at a given 
I. In other words, even if the angular correlation function is featureless, the angular spectrum may 
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universe 1 



Figure 3.1: Ensemble Average of Angular Correlation Function 

A schematic view of the ensemble average of the n-point angular correlation function, 
/(ni)/(na)/(n3) . . . /(n„). We measure it on each universe, and then average it over many 
universes. 

have a distinct structure, for inflation predicts a prominent peak in the angular power spectrum, not 
in the angular correlation function. In this chapter, we study statistical properties of the angular 
n-point harmonic spectra. 

3.1 Statistical Isotropy of the Universe 

In reality, we cannot measure the ensemble average of the angular harmonic spectrum, but one 
realization such as a« imi a/ 2m2 ■ . .ai nrrin , which is so noisy that we want to average it somehow to 
reduce the noise. 

We assume statistical isotropy of the universe from which it follows that our sky is isotropic 
and has no preferred direction. Isotropy of CMB justifies the assumption. The assumption readily 
implies that one can average the spectrum over rrii with an appropriate weight, as m; represent an 
azimuthal orientation on the sky. The average over rrii enables us to reduce statistical error of the 
measured harmonic spectra. 

How can we find the weight? One finds it as a solution to statistical isotropy, or rotational 
invariance of the angular correlation function on the sky, 

(D/(nOZ)/(n 2 ) . . . Df(h n )) = </(fii)/(ft 2 ) • • • f(n n )) , (3.3) 

where D = D(a,j3, 7) is a rotation matrix for the Euler angles a, j3, and 7. Figure |3.2| sketches 
the meaning of statistical isotropy. Substituting equation ( |3.2| ) for /(n) in equation ( |3.3| ), we then 
need rotation of the spherical harmonic, DY[ m (h). It is formally represented by the rotation matrix 
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our sky 



Figure 3.2: Statistical Isotropy of Angular Correlation Function 

A schematic view of statistical isotropy of the angular correlation function. As long as its 
configuration is preserved, we can average /(ni) . . . /(n n ) over all possible orientations and 
positions on the sky. 



element, D^, m (a, (3, 7), as ( [Rotenberg et al., 1959| ) 



I 



DY lm (n)= J2 D m'm Y lm'(n). (3.4) 



m'=—l 



The matrix element, D~i m = (l,m' \D\l,m), describes finite rotation of an initial state whose 
orbital angular momentum is represented by I and m into a final state represented by I and ml '. 
Finally, we obtain the statistical isotropy condition on the angular n-point harmonic spectrum: 

(a kmi a hm2 ■ ■ ■ ai nmn ) = ^ (o, im / a hm > 2 . . . a /nm ; ) D { ^ mi . . . . (3.5) 

all m' 



Using this equation, Hu (2001 ) has systematically evaluated appropriate weights for averaging 
the angular power spectrum (n = 2), bispectrum (n = 3), and trispectrum (n = 4), over azimuthal 
angles. Some of those may be found more intuitively; however, this method allows us to find the 
weight for any higher-order harmonic spectrum. In the following sections, we derive rotationally 
invariant, azimuthally averaged harmonic spectra for n = 2, 3, and 4, and study their statistical 
properties. 

3.2 Angular Power Spectrum 

The angular power spectrum measures how much fluctuations exist on a given angular scale. For 
example, the variance of ai m for I > 1, (a| m a* m ), measures the amplitude of fluctuations at a given 
I. 
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Generally speaking, the covariance matrix of ai m , \ a>i 1 mt a i 2m2 /> * s no ^ necessarily diagonal. It 
is, however, actually diagonal once we assume full sky coverage and rotational invariance of the 
angular two-point correlation function, as we will show in this section. The variance of ai m thus 
describes the two-point correlation completely. 

Rotational invariance (Eq.( |3.5| )) requires 

( fl f.»i a U) = E ( a ^s a U' 2 ) jD s ) mi Z) !n> 2 ( 3 - 6 ) 

to be satisfied, where we have used the complex conjugate for simplifying calculations. From this 
equation, we seek for a rotationally invariant representation of the angular power spectrum. Suppose 
that the covariance matrix of a/ m is diagonal, i.e., ^Oi imi a^ m2 ^> = (C^) $i 1 i 2 $m 1 m2- Equation ( |3.6[) 
then reduces to 

("iimi^mj) = ( C h) Kh X! ^Smi^mj = ( C h) S hh S mim 2 - (3.7) 

Thus, we have proven (C;) rotationally invariant. Rotational invariance implies that the covariance 
matrix is diagonal. 



3.2.1 Estimator 



Observationally, the unbiased estimator of (C/) should be 

1 ' 1 ( 2 ' * \ 

1 l \ m=l / 

III 



1 



21 + 1 



m=l 

The second equality follows from a/_ m = a* m (— l) m , i.e., ai_ m a[_ m = a; m a ; * m , and hence we average 
2Z + 1 independent samples for a given £. It suggests that fractional statistical error of Cj is reduced 
by y/ 1/(21 + 1). This property is the main motivation of our considering the azimuthally averaged 
harmonic spectrum. 

We find it useful to define an azimuthally averaged harmonic transform, e;(n), as 

/ i 

ei ^ = V2lTi ^ a imYi m (n), (3.9) 

' rn=—l 

which is interpreted as a square-root of C\ at a given position of the sky, 

^ ef(n) = Q. (3.10) 

e;(n) is particularly useful for measuring the angular bispectrum ( Spergel and Goldberg, 1999| ; |Ko-| 
|matsu et al., 2001b| ) (chapter |5|) , trispectrum (chapter |6|) , and probably any higher-order harmonic 
spectra, because of being computationally very fast to calculate. This is very important, as the 
new satellite experiments, MAP and Planck, have more than millions of pixels, for which we will 
crucially need a fast algorithm of measuring these higher-order harmonic spectra. 
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3.2.2 Covariance matrix 

We derive the covariance matrix of Q, (C/CV) — (Ci) {Ci>}, with the four-point function, the trispec- 
trum. Starting with 

(QCV) = (2/ + l)(2l' + 1) ^ ( aima *i™ a i'™' a *i'rn') > ( 3 - n ) 
we obtain the power spectrum covariance matrix 



(CjCf) - (Oi) (Cj.) = ^fr + (a + 1) ' (2l , + 1) I> 



mO-lm, a l 'm ' a V 'm' 



in in' 



where (azmaf m az'm/af/ m /) c is the connected four-point harmonic spectrum, the connected trispec- 
trum, which is exactly zero for a Gaussian field. It follows from this equation that the covariance 
matrix of Ci is exactly diagonal only when ai m is Gaussian. (T 1 ^ 2 (L)\ is the ensemble average of 



the angular averaged connected trispectrum, which we will define in § |3.4| (Eq.( |3.24| )). 

Unfortunately, we cannot measure the connected T//,,(0) directly from the angular trispectrum 
(see § 0). We will thus never be sure if the power spectrum covariance is precisely diagonal, as long 
as we use the angular trispectrum. We need the other statistics that can pick up information of the 
connected T/,^(0), even though they are indirect. Otherwise, we need a model for the connected 
trispectrum, and use the model to constrain the connected TjL(0) from the other trispectrum 
configurations. We will discuss this point in chapter |6|. 

There is no reason to assume the connected T//,,(0) small. It is produced on large angular 
scales, if topology of the universe is closed hyperbolic ( Inoue, 2001b| ). In appendix |B|, we derive 



an analytic prediction for the connected trispectrum produced in a closed hyperbolic universe. On 
small angular scales, several authors have shown that the weak gravitational lensing effect produces 
the connected trispectrum or four-point correlation function ( |Bernardeau, 1997 ; Zaldarriaga and 



[Seljak, 1999^ Zaldarriaga, 2000 ); Hu (2001] ) finds that the induced off-diagonal terms are negligible 



compared with the diagonal terms out to / ~ 2000. 

If the connected trispectrum is negligible, then we obtain 



(QQ,) - (Q) (Q,) ^ ^j^-Su, (3.13) 



The fractional error of C/ is thus proportional to y/ 1/(21 + 1), as expected from our having 21 + 1 
independent samples to average for a given /. The exact form follows from C\ being x 2 distribution 
with 21 + 1 degrees of freedom when a\ m is Gaussian. If a/ m is Gaussian, then its probability density 
distribution is 

P(aim) = ew[ ~f J{2m] . (3.14) 
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Figure 3.3: Angular Bispectrum Configuration 

We use this distribution to generate Gaussian random realizations of a\ m for a given (Cj). First, 
we calculate (Ci) with the CMBFAST code (Seljak and Zaldarriaga, 1996) for a set of cosmological 



1/2 

parameters. We then generate a realization of a/ m , a; m = e (C7) ' , where e is a Gaussian random 
variable with the unit variance. 



3.3 Angular Bispectrum 

The angular bispectrum consists of three harmonic transforms, ai imi ai 2rn2 ai 3rn3 . For Gaussian ai m , 
the expectation value is exactly zero. By imposing statistical isotropy upon the angular three-point 
correlation function, one finds that the angular averaged bispectrum, Bi ± i 2 i 3 , given by 



{ahm 1 ai 2 m2 a hm 3 ) = {Blihh) (3.15) 




satisfies rotational invariance (Eq,( |3.5| )). Here, the matrix denotes the Wigner-3j symbol (see 
appendix |B|). Since l\, I2, and I3 form a triangle, Bi x i 2 i 3 satisfies the triangle condition, \h — lj\ < 
Ik < h + lj for all permutations of indices. Parity invariance of the angular correlation function 



demands h + h + h = even. Figure 13^ sketches a configuration of the angular bispectrum. 

The Wigner-3j symbol, which describes coupling of two angular momenta, represents the az- 
imuthal angle dependence of the angular bispectrum, for the bispectrum forms a triangle. Sup- 
pose that two "states" with (li,mi) and (12,1112) angular momenta form a coupled state with 
(^3, 771,3). They form a triangle whose orientation is represented by mi, 7772, and 7773, with satisfying 
mi + "72 + 7773 = 0. As we rotate the system, the Wigner-3j symbol transforms m's, yet preserving 
the configuration of the triangle. Similarly, rotational invariance of the angular bispectrum de- 
mands that the same triangle configuration give the same amplitude of the bispectrum regardless 
of its orientation, and thus the Wigner-3j symbol describes the azimuthal angle dependence. 

The proof of (B^i^) to be rotationally invariant is as follows. Substituting equation ( 3.15| ) for 
the statistical isotropy condition (Eq. (|3.5|) ) for 77 = 3, we obtain 



3 m-s I 
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Y (a hm > a hm > a w \ D {h ) D {l2 ) D {h ) 

/ / \ U" 1 ! t2"»2 '3 m 3 / mjmi m^rr^ 171^3 



all m' 



(■Biiiais) Y f J/ 
all m' V m l 



*2 

to 2 m' 3 



x y , , , 2 (2L + 1)Z^ *^ ' 

= <^ 3 > f ^ /2 ' 3 V (3-16) 
V mi 7712 rri3 y 

In the second equality, we have reduced -P^/ -P^m? to Djffi^ , using equation ( p. 17] ). In the 
third equality, we have used the identity ( Rotenberg et al., 1959 ), 

/ h h h \ ( h h L \ = Si 3L S m ' 3 M' 
J^, { m[ m' 2 m> 3 ) { m[ m' 2 M' ) " 2L + 1 

3.3.1 Estimator 

To obtain the unbiased estimator of the angular averaged bispectrum, B^i^, we invert equa- 
tion ( p. 15 ) with the identity ( |3,17| ), and obtain 

B hhh = Y ( 1 2 3 ) a hm 1 ai 2 m 2 ai 3m:i . (3.18) 



all m 



mi 77t2 m3 



We can rewrite this expression into a more computationally useful form. Using the azimuthally 
averaged harmonic transform, ej(n) (Eq. (|3.9| )), and the identity ( Rotenberg et al., 1959| ), 



-i 



h h h \ t h h h \ / W 

mi m 2 m 3 J \ J \j (2Z X + 1) (2l 2 + 1) (2i 3 + 1 

d 2 h 



Y limi (h)Yi 2m2 (ii)Y hm3 (h), (3.19) 



we rewrite equation ( |3.18| ) as 



Bhiah = y q q q J y ^ e h (h)e l2 (h)e l3 {h). (3.20) 

This expression is computationally efficient; we can quickly calculate e;(n) with the spherical har- 
monic transform. Then, the average over the full sky, / cZ 2 n/(47r), is done by the sum over all 
pixels divided by the total number of pixels, iV -1 Y^f , if all the pixels have the equal area. Note 
that the integral over n must be done on the full sky even when a sky-cut is applied, as e/(n) 
already encapsulates information of partial sky coverage through ai m , which may be measured on 
the incomplete sky. 
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3.3.2 Covariance matrix 



We calculate the covariance matrix of Bi x i 2 i 3 , provided that non-Gaussianity is weak, (-8^2/3) ~ 0. 
Since the covariance matrix is a product of six a/ m 's, we have ^2 ■ 4C2/3! = 15 terms to evaluate, 
according to the Wick's theorem; however, using the identity ( Rotenberg et al,, 1959| ), 



(-1)' 



-1) 



-m 



V2TTT 



Si'o, 



(3.21) 



and assuming none of Ts zero, we find only 3! = 6 terms that do not include (0/ 



i a ljrrij 



but 



include only (a; iWi a|. mj y non- vanishing. Evaluating these 6 terms, we obtain ( Luo, 1994 ; Heavens, 



1998| ; |Spergel and Goldberg, 1999| |Gangui and Martin, 2000| ) 



B hhh B l'M 



all mm 1 



n '2 '3 
771,1 fl2 ^3 



777i 



(C h ) (C h ) (C l3 ) 



v v v v v v v v v 

X'l'2'3 1 £ l 3 l l l 2 I £ l 2 l 3 l l 

°hhh + °hhh + °' 




a limi a l2m2 a l3m 3 a l'm', a l' r ,mL a l',m'o 



+ r_l\h+h+l 3 f £13 2 + £213 + A21 
'li2'3 1 V / \^ il'2'3 il'2'3 il'2'3 / 



(3.22) 



1' V V 

where Milil 



Si^Si^S^i' , and so on. Hence, the covariance matrix is diagonal in the weak 
non-Gaussian limit. The diagonal terms for li ^ and h + 12 + h 



even are 



B 



hhh 



(Qi) (Q 2 ) {Ci 3 ) (1 + 25 h i 2 5i 2 i 3 + S h i 2 + 6i 2 i 3 + 5i 3h ] 



(3.23) 



The variance is amplified by a factor of 2 or 6, when two or all Vs are same, respectively. 

We find that equation ( |3.23j ) becomes not exact on the incomplete sky, where the variance 
distribution becomes more scattered. Using simulated realizations of a Gaussian sky, we have 
measured the variance on the full sky as well as on the incomplete sky for three different Galactic 
sky-cuts, 20°, 25°, and 30°. Figure 3.4 plots the results; we find that equation ( |3.23 ) holds only 
approximately on the incomplete sky. 



3.4 Angular Trispectrum 



The angular trispectrum consists of four harmonic transforms, ai irni ai 2m2 ai 3m3 ai im4 . |Hu (20011 ) 
finds a rotationally invariant solution for the angular trispectrum as 



{a hmi a l2m2 a l3m3 a hm4 ) = V(j 1 J J „ J(Y W) • (3-24) 




One can prove this solution, (rj^(L)J, rotationally invariant by similar calculations to those 
proving the angular bispectrum to be so. By construction, 1%, I2, and L form one triangle, while 
Z3, Z4, and L form the other triangle in a quadrilateral with sides of l\, I2, Z3, and I4. L represents 
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Figure 3.4: Variance of Angular Bispectrum 

Histograms of variance of the angular bispectrum for h < h < h up to a maximum 
multipolc of 20. There are 466 modes. These are derived from simulated realizations 
of a Gaussian sky. The top-left panel shows the case of full sky coverage, while the 
rest of panels show the cases of incomplete sky coverage. The top-right, bottom-left, 
and bottom- right panels use the 20°, 25°, and 30° Galactic sky-cuts, respectively. 
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Figure 3.5: Angular Trispectrum Configuration 



a diagonal of the quadrilateral. Figure 3.5 sketches a configuration of the angular trispectrum. 
When we arrange l±, I2, h, and I4 in order of l\ < 1% < Z3 < I4, L lies in max(/2 — h,h — h) < 
L < min(Zi + ^2^3 + I4). Parity invariance of the angular four-point correlation function demands 
h + h + L = even and Z3 + I4 + L = even. 

The angular trispectrum generically consists of two parts. One is the unconnected part, the 
contribution from Gaussian fields, which is given by the angular power spectra ( |Hu, 2001] ), 

34 / unconnected 
h+h, ' 

V ^'i "T" J -;v.'"*3 -r v \w 3 / u iih u kii4 u ^> 

(3.25) 



(_l)^^3 v / (2Zl + 1)(2Z g + 1} (Cja) 6 hh 6 l3h 6 L0 
+ (2L + 1) (CjJ (C l2 ) U-l) h+h+L S h i 3 6i ali + <5yA* 3 



For l\ < I2 < Z3 < I4, the unconnected terms are non-zero only when L = or l\ = h = h = h- 
We have numerically confirmed that our estimator given below (Eq.( |3.30| )) accurately reproduces 
the unconnected terms (Eq.( |3.25| )) on a simulated Gaussian sky. 

The other is the connected part whose expectation value is exactly zero for Gaussian fields; 
thus, the connected part is sensitive to non-Gaussianity. When none of Vs are same in 
one might expect the trispectrum to comprise the connected part only; however, it is true only on 
the full sky. The unconnected terms on the incomplete sky, which are often much bigger than the 
connected terms, leak the power to the other modes for which all Z's are different. We should take 
this effect into account in the analysis. 



3.4.1 Estimator 

Inverting equation (3.24), we obtain the unbiased estimator of Tj^(L) ( ]Hu, 2001 ), 



all m M \ 



m\ 771-2 M 



h h L 
7773 7774 — M 



X O/xmi G>; 2 rri2 ^l^m^ O'l^rrn ■ 



(3.26) 
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Note that this expression includes both the connected and the unconnected terms. 

We find that this estimator has a special property for L = (which demands l\ = I2 and 

l lah 

power spectrum estimators, C^C^, 



h = h)- The trispectrum estimator for these configurations, T} 1 } 1 ^), reduces to a product of two 



^hh^ ~ ( m _m, f) ) ( mo _ m „ fl ) ' imi ^ ~ mi a/ 3"l 3 % -m 3 





= (-ir +i3 V(2Zi + l)(2Z 3 + l)a i Q 3 , (3.27) 

where C\ = (2Z+1) -1 Y^ m aj m a* m . We have used the identity, equation (|3.21| ), and a;_ m = (— \) m a\ m 
in the second equality. From this equation, one may assume that T}^(0) coincides with the 
unconnected terms for L = (see Eq.( 3.25[ )), 



(^3 (°))_ ctcd = (-l)' 1+ V(2Zi + l)(2Z 3 + l) W) W + 2 (C/i) 2 <W (3-28) 

They are, however, different for non-Gaussian fields, because of the power spectrum covariance, 
equation ( 3.12j ). By taking the ensemble average of T^(0), and substituting equation ( 3.12j ) for 
(C/ 1 C; 3 ), we find a rather trivial result: 

(^t(O)) = (-l) il+i3 V /(2/i + l)(2/3 + l)<C Jl C j3 ) 

= (-l)h+h^2h + 1)(2I 3 + 1) (Qx) (Ci s ) + 2<<7, 1 > 2 «U + (t££(0) 

= («(0))_ cted + («(0)) c - (3-29) 

Hence, T^(0) contains information not only of the unconnected trispectrum, but also of the 
connected trispectrum. 

Unfortunately, we cannot measure the connected part of T^(0) directly from the angular 
trispectrum because of the following reason. To measure the connected terms, we have to subtract 
the unconnected terms from the measured trispectrum first. Since we are never able to measure 
the ensemble average of the unconnected terms (Eq. fl3.28|) ), we estimate them by using estimated 
power spectrum, C\. If we subtract the estimated unconnected terms, oc C; 1 Q 3 , from measured 
T/ 1 / 1 ^), then it follows from equation ( |3~27j ) that lj£(0) vanishes exactly: lj£(0) = 0; thus, 
(0) has no statistical power of measuring the connected terms. 

For practical measurement of the angular trispectrum, we rewrite the trispectrum estimator 
given by equation ( 0.26 ) with the azimuthally averaged harmonic transform, e;(n) (Eq. (|3.9f )). We 
find that the following form is particularly computationally efficient: 

^ = ^7 E Ozm, (3-30) 



2L + 1 M=-L 



where t 1 ^ is given by 




J d 2 n[e h (h)e h (h)]Yl M (h). (3.31) 
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Since t^ M * s the harmonic transform on the full sky, we can calculate it quickly. This method 
makes measurement of the angular trispectrum computationally feasible even for the MAP data in 
which we have more than millions of pixels; thus, the methods developed here can be applied not 
only to the COBE DMR data, but also to the MAP data. 



3.4.2 Covariance matrix 

We calculate the covariance of the trispectrum in the weakly non-Gaussian limit. Since the trispec- 
trum covariance comprises eight a; m 's, the total number of terms is sCVeCV 4C2/ '4! = 105 according 
to the Wick's theorem. The full calculation will be a nightmare, for we have to deal with 



(2L + 1)(2L' + 1) £ S (-1 

all mm' MM' 

L \( h k 

M J \ TO3 7TT-4 
* 

&l\m\ 0-l2m2 ^Izmz ^l^m^l 1 



,M+M' 





(3.32) 



We reduce this intricate expression to much more simplified forms for some particular configura- 
tions. For L,L' 7^ terms, thanks to the identity ( p.21| ), only 4! = 24 terms that do not include 
(ai i m i a >ljmj S } but include only (ai i m i a i :jm ^ are non- vanishing. For L = V = terms, the triangle 
conditions in a quadrilateral demand h = I2 and I3 = I4 (see figure |3~5|) . As we have shown, these 
configurations have no statistical power of measuring the connected trispectrum of interest. Hence, 
we evaluate L, L' ^ terms in the following. 

Evaluating 24 L, L' 7^ terms is still a headache; however, for l\ < 1% < Z3 < I4, we have only 8 
terms left: 



V V 



(2L + l)(C h )(C l2 ) (C l3 )(C u ) 



+(-1) 



xhfalsU 
<J,i 11 11 11 

l l l 2 l 3 l 4 

h+l 2 +L 



+ °l'l 



,' V V 



+ (-1 



^1+^2+^3+^4 



I V V V V ' I' VVV I ' \ 

\ '2134 '4312/ 



( shhl-ili 

1 Uu U II II 

\ '2'l'4 l 3 

\l 3 +h+L 



+ 5} 



I1I2I3I4, 



V V V V 

'4 3 2*1 



/ shhteh 
\ I' I' I' I' 

\ l \ l 2 4 3 



4 _|_ fihfohh} 



V V V V I 

3 4 2 1 / 



(3.33) 



V V V V 

where W^sU = Kl'Jhl^hl^hl'^ and so on - Using parity invariance, h + l 2 + L = even and 
/3 + Ia + ^ = even, we find the covariance matrix diagonal. Thus, the diagonal terms for L / and 
h < h < ^3 < ^4 are simplified very much as 



' hi a 



(2L + 1) (C h ) (C l2 ) (C h ) (C U ) (1 + 5 hh + S hU + 5 w Ak) 



(3.34) 



This result is strictly correct only on the full sky; the incomplete sky makes the variance distribution 
much more scattered. Figure |3.6| plots the variance on the full sky as well as on the incomplete sky. 
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Figure 3.6: Variance of Angular Trispectrum I 

Histograms of variance of the angular trispectrum for L =/= and l\ < I2 < h < h, for 
which the unconnected terms vanish on the full sky. There are 16,554 modes, up to a 
maximum multipole of 20. The meaning of the panels is the same as in figure |3.4|. 
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For the rest of configurations for which the unconnected terms vanish, L / 0, I2 = I3, and 
l\ 7^ I4, the covariance matrix is no longer diagonal in L,L' ( |Hu, 2001 ). Figure |3.7| plots the 
numerically evaluated variance on the full sky as well as on the incomplete sky. The variance 
divided by (2L + 1) (C^) (C/ 2 ) (Q 3 ) (Cj 4 ) is no longer an integer, but more scattered than that for 
L ^ and l\ < I2 < I3 < £4 even on the full sky. 



3.5 Power Spectrum and Bispectrum on the Incomplete Sky 

Incomplete sky coverage destroys orthonormality of the spherical harmonics on the sky. The degree 
to which orthonormality is broken is often characterized by the coupling integral (Peebles, 1980| ), 

Wjjw = / d 2 n W(n)Y^ (n) Y Vm , (n) = / d 2 n Y lm (n) K,, m , (n) , (3.35) 

where W(n) is zero in a cut region otherwise 1, and o b s denotes a solid angle of the observed sky. 
When Wu 'mm' 7^ fiu'fimm') the measured harmonic transform of the temperature anisotropy field, 
ai m , becomes a biased estimator of the true harmonic transform, a|™ 6 , through 

00 V 

O-lm = X! X! a \'m'Wu<mm'- (3.36) 
J'=0 m'=-V 

Hence, we must correct our estimators of the power spectrum and the bispectrum for the bias 
arising from incomplete sky coverage. 

First, we derive a relationship between the angular power spectrum on the incomplete sky and 
that on the full sky. Taking the ensemble average of the estimator of the power spectrum, the 
pseudo-C; (|Wandelt et al., 1998| , |2Q0QQ , Q = (21 + l)' 1 £ m \a tm \ 2 , we have 

1 

21 + 



1 . 

~ " i' mm' 



1 c tme £ W d 2 n W W Y *m (") ^'m' (n) / d 2 m VF (m)l, m (m) y; m , (m) 

m i'm' 

1 c true £ /• d 2 A W{m * m (n) /" d 2 m VF(m)y im (m) 5^ (n - m) 

T J- ™ ./ J 



21 , . 

m 

2a 



Q trUe /^^(n)^(l) 



= cf uc ^. (3.37) 

In the second equality, we have taken C]f ue out of the summation over I', as \Wui mm i\ 2 peaks very 
sharply at I = I', and Cjf ue varies much more slowly than |W/// mm '| 2 in I'. This approximation is 
good for nearly full sky coverage. In the third equality, we have used J2i'm' Yi'm' ^* m ' (™) = 
(n — rh). In the forth equality, we have used J2m ^Cm (p) Y\ m (A) = ^^-Pi(n ■ m). The result 
indicates that the bias amounts approximately to a fraction of the sky covered by observations. 
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Figure 3.7: Variance of Angular Trispectrum II 

Histograms of variance of the angular trispectrum for L ^ 0, I2 = h, and l\ ^ U, for 
which the unconnected terms vanish on the full sky. There are 4,059 modes, up to a 
maximum multipole of 20. The meaning of the panels is the same as in figure |3.4|. 



3.5. POWER SPECTRUM AND BISPECTRUM ON THE INCOMPLETE SKY 



55 



Next, we derive a relationship between the angular bispectrum on the incomplete sky and that 
on the full sky. We begin with 



{ai imi ai 2m2 ai 3m3 ) - ^ \ a z;m; a Z^ a |m^/ ^J^mim^ ^M^m^si^m^ • (3.38) 
all I'm' 

Rotational and parity invariance of the bispectrum implies the bispectrum given by 
\ a liirii Q-Z 2 m 2 a hm s ) = kMs / d 2 n Y l 

imi (A) ^Z 2 m 2 ^3m3 (A) > (3.39) 

where fez^^ is an arbitrary real symmetric function, which is related to the angular averaged 
bispectrum, Bi x i 2 i 3 , When b^f^ varies much more slowly than the coupling integral, we obtain 

(a hmi ai 2m2 a hm3 ) = ^ ^ / d 2 fi y^ w , (A) Y^ m , (A) y| m , (A) 

all i' all m' 

x J d 2 ^ W^Y^ (Ai)y^ TOl (Ai) 
x I d 2 n 2 W(n 2 )y,, m , (n 2 )y^ m2 (n 2 ) 
x J d 2 n 3 W(n 3 )y^ m ^ (n3)y^ m3 (n 3 ) 
~ tt 3 / W(m* imi (n) ^ m2 (A) (A) . (3.40) 

Then, we calculate the angular averaged bispectrum, Bi x i 2 i 3 (Eq. (|3.15| )). By convolving equa- 
tion ( p. 18 ) with the Wigner-3j symbol and using the identity ( |3.19| ), we obtain 

/r \ ~ At™e / 4vr ( h l 2 k \ 1 



x 

all m 



(2Zi + l)(2i 2 + l)(2J 3 + l) v o o o ; 
imi ( m ) y 2 m 2 

(m) y am3 (m) 

ill m 

x | d 2 A ty(n)y,* mi (n) y,* m2 (n) y^ ma (n) 

true /(2Zi + l)(2Z 2 + l)(2Z 3 + l) / /i / 2 Z 3 \ 1 



4vr ^ J 

I IT l^ W{il)Ph (* -A) A, -A) 



/, true , 



\2h + 1)(2Z 2 + 1)(2Z 3 + 1) / h h h \ ^obs 



4vr 4t 



otrue ^obs AT\ 



where we have used the identity, 

2 



£Y fl 1 (x)fl 2 (x)fl 3 (x)=( J q q j • (3-42) 
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Thus, the bias for the angular bispectrum on the incomplete sky is also approximately given by a 
fraction of the sky covered by observations. 



Chapter 4 

Theoretical Predictions for the CMB 
Bispectrum 



In inflation, quantum fluctuations of scalar fields generate the matter and the radiation fluctuations 
in the universe (chapter ^). In the stochastic inflationary scenario of [Starobinsky (1986 ), quantum 
fluctuations decohere to generate classical fluctuations. There are two potential sources of non- 
Gaussianity in this inflationary model: (a) non-linear coupling between the classical inflaton field 
and the observed fluctuation field, and (b) non-linear coupling between the quantum noise field 
and the classical fluctuation field. Salopek and Bond (1990| , |1991| ) have studied the former; Gangui 
et al. (1994) have studied the latter. 

Calzatta and Hu (1995) and Matacz (1997a| ,b|) present an alternative treatment of the decoher- 
ence process that leads to different results for the primordial density perturbations from [Starobiri- 



|sky (1986 ). Matacz's treatment makes similar predictions for the level of non-Gaussianity to the 
Starobinsky's treatment ( Matacz, 1997a| Jb|). These studies conclude that in a slow-roll regime, 
fluctuations are Gaussian; however, features in a inflaton potential can produce significant non- 
Gaussianity ( [Kofman et al., 1991 ). 

Previous work on the primary non-Gaussianity has focused on very large angular scales, the 
COBE scale, where the temperature fluctuations trace the primordial fluctuations. For MAP and 
Planck] however, we need the full effect of the radiation transfer function. In this chapter, we 
develop a formalism for doing this, and then present numerical results. Both the formalism and 
the numerical results are main results of this chapter. We also discuss how well we can separate 
the primary bispectrum from various secondary bispectra. 

This chapter is organized as follows. In § 4T, we define the angular bispectrum, the Gaunt 
integral, and a new quantity called the reduced bispectrum, which plays a fundamental role in 
estimating physical properties of the bispectrum. In § 4.2, we formulate the primary bispectrum that 
uses the full radiation transfer function, and presents numerical results for the primary bispectrum 
and skewness. In § 4.2, we calculate secondary bispectra from the coupling between the Sunyaev- 
Zel'dovich effect and the weak lensing effect ( fSpergel and Goldberg, 1999| ; Goldberg and Spergel 



1999; Cooray and Hu, 2000), and from extragalactic radio and infrared sources. In 



we 
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study how well we can measure each bispectrum, and how well we can discriminate between those 



bispectra. § 4.5 is devoted to further discussion and our conclusions in this chapter. 



4.1 Reduced Bispectrum 

We expand the observed CMB temperature fluctuation field, AT(n)/T, into the spherical harmon- 



ics, 



^AT(n) 



a lm = J d 2 h^^YC m (h), (4.1) 
where the hats denote unit vectors. The CMB angular bispectrum is given by 

B hhh 3 = ^ a hm 1 a l2rn2 a hm'j,) > (4-2) 

and the angular averaged bispectrum is (Eq. (|3.18| )) 

B hhh = E { h h h ) B hkh ma > ( 4 - 3 ) 

aU m V ^1 ™2 m 3 j hhh 

where the matrix is the Wigner-3j symbol. The bispectrum, B™^ 2 ™ 3 , satisfies the triangle con- 
ditions and parity invariance: mi + m<i + 7713 = 0, l\ + I2 + h = even, and \h — lj\ < Ik < h + lj 
for all permutations of indices. It implies that B™^ 2 ™ 3 consists of the Gaunt integral, G™l™^ mz > 
defined by 



'(2Zi + l)(2Z 2 + l) (2^3 + 1) / h h h\ h h h 




Air \ / \ mi 777,2 7713 



(4.4) 



G™i^ 2Tn3 is real, and satisfies all the conditions mentioned above. 

Rotational invariance of the angular three-point correlation function implies that Bi x i 2 i z is writ- 
ten as 

where 6/ 1 i 2 ; 3 is an arbitrary real symmetric function of l\, I2, and Z3. This form, equation (|4.5| ), is 
necessary and sufficient to construct generic H™^ 2 ™ 3 under rotational invariance; thus, we will use 
6; 1 ; 2 ; 3 more frequently than S™^ 2 ™ 3 in this chapter, and call this function the reduced bispectrum, 
as 6; 1 ; 2 j 3 contains all physical information in B™l™ 2in3 . Since the reduced bispectrum does not 
contain the Wigner-3j symbol, which merely ensures the triangle conditions and parity invariance, 
it is easier to calculate physical properties of the bispectrum. 

We calculate the angular averaged bispectrum, B^i^, by substituting equation ( [i~5| ) into ( |t.3| ). 



! {2h + l)(2Z 2 + l)(2/ 3 + l 




B hhh = \i ^ I n n n iKhh, (4-6) 
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where we have used the identity, 



V- ( h h h \ Qmi m 2 m 3 = / (2Zl + l)(2/ 2 + l)(2k + 1) ( h h h\ 

£ln \ mi m2 m 3 ) hhh "V 4vr 1^000^' l ' j 

Alternatively, one can define the bispectrum in the flat-sky approximation, 

(a(h)a(h)a(l 3 )) = (2vr)V 2 ) (h + 1 2 + 1 3 ) B(h, 1 2 , 1 3 ), (4.8) 

where 1 is a two-dimensional wave vector on the sky. This definition of £?(li,l 2 ,l3) reduces to 
equation witn th e correspondence, G™l 2 ™ 3 2m3 -> (2vr) 2 (5 (2) (h + 1 2 + 1 3 ), in the flat-sky limit 

(|Hu, 20001) . Thus, we have 

h^iz w -B(li, 1 2 , 13) (flat-sky approximation) . (4.9) 

This fact motivates our using the reduced bispectrum, b^jjjg, rather than the angular averaged 
bispectrum, B^i^. Note that fej^g is similar to Bi^k defined by Magueijo ( pVlagueijo^ lOOOD ; the 
relation is b hhh = \/4nvB hhh . 



4.2 Primary Bispectrum and Skewness 



4.2.1 Model of the primordial non-Gaussianity 

If primordial fluctuations are adiabatic scalar fluctuations, then 



a li. 



4ir(-i) 1 



d 3 k 

(2^)3 



$(k) 5T ^)^(k), 



(4.10) 



where $(k) is the primordial curvature perturbation in Fourier space, and gTi(k) is the radiation 
transfer function. ai m takes over the non-Gaussianity, if any, from ^(k). Although equation ( 4.10 ) 
is valid only if the universe is flat, it is straightforward to extend this to an arbitrary geometry. 
We can calculate the isocurvature fluctuations similarly by using the entropy perturbation and the 
proper transfer function. 

In this chapter, we explore the simplest weak non-linear coupling case, 



$(X) = $ L (x)+/ NL $ 2 (X) 



(4.11) 



in real space, where $l(x) denotes a linear Gaussian part of the perturbation, and (<3?(x)) = is 
guaranteed. Henceforth, we call /nl the non-linear coupling parameter. This model is based upon 
slow-roll inflation; [Salopek and Bond (1990 , 1991 ) and |Gangui et al. (1994 ) have found that /nl 
is given by a certain combination of slope and curvature of an inflaton potential (<3? 3 = — 2/nl in 
Gangui et al. (1994| )). pangui et al. (1994 ) have found that |/nl| ~ 10 -2 for quadratic and quartic 
potential models. We have reproduced these results in chapter |2[ 
In Fourier space, we decompose 3>(k) into two parts, 



$(k) = $ L (k) + $ NL (k), 



(4.12) 
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and accordingly we have 



"Zm ~r u lm > 



where 3>NL(k) is a non-linear curvature perturbation defined by 

dh 



$NL(k) = /nl 



/ (2^ <I>L(k + P) ^(P) " (2^) 3 ^ (3) (k) (*£(x) 



(4.13) 



(4.14) 



One can confirm that (<&(k)) = is satisfied. In this model, a non-vanishing component of the 
<3?(k)-field bispectrum is 

($ L (k 1 )$ L (k 2 )$ NL (k 3 )) = 2(27r)V 3 )(k 1 +k 2 + k 3 )f NL P t ,(k 1 )P„(k 2 ), (4.15) 

where P&(k) is the linear power spectrum given by 

($ L (ki)$ L (k 2 )) = (27T) 3 P <s> (k 1 )5^(k 1 + k 2 ). (4.16) 

We have also used 



<<& L (k + p)$J,(p)> = (27r) 3 P*(p)<f( 3 >(k), 



and 



(2^)~ 3 / d 3 kP$(fc). 



(4.17) 



(4.18) 



Substituting equation ( 4.10| ) into ([4.2]), using equation (|4.15| ) for the <I>(k)-field bispectrum, and 
then integrating over angles ki, k3, and k3, we obtain the primary CMB angular bispectrum, 



L L NL \ , / L NL L \ , / NL L I 

H 1 m 1 u 'l2m2 U 'hm i I ' \ "Zimi u Z2"i2 u 'hm 3 I < \ Uj l 1 m\ u 'l2m2 U, l 



L 

3 "13 



(,NL/,„^L /„\iX 



r 2 dr 







6fc(r)6fe(r)C(0 + 6fc(r)CW(»- 



NL/ 



where 



#(r) 



7T 70 

2 



^dkP^gvikMkr), 
k 2 dkf NL g T i(k)ji(kr). 



Note that 6p(r) is dimensionless, while 6^ L (r) has a dimension of L -3 . 

One confirms that equation ( |4.5| ) holds; thus, the reduced bispectrum, b^i^ (Eq. 
primordial non-Gaussianity is 



(4.19) 

(4.20) 
(4.21) 

|)), for the 



.primary 
°hhh 



r 2 dr 



NL/ 



+C(0^W^ 3 (r) 



(4.22) 



We can fully specify bf*}™ Ty by a single constant parameter, /nl, for the CMB angular power 
spectrum, C7j, will precisely measure cosmological parameters ( Pond et al., 1997| ). We stress again 
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that this is the special case in the slow-roll limit. If the slow-roll condition is not satisfied, then 
/nl = /nl(&i, &2, ^3) at equation Q4.15| ) ( Gangui et al., 1994| ), |Wang and Kamionkowski (2000| ) 
have developed a formula to compute -Bz 1 i 2 ; 3 from generic <I>(k)-field bispectrum. Our formula 
(Eq,( 4~19|) ) agrees with theirs, given our form of the <£(k)-field bispectrum (Eq.( 4~15| )), 

Even if inflation produces Gaussian fluctuations, the general relativistic second-order pertur- 
bation theory produces /nl ~ 0(1) ( Pyne and Carroll, 1996|) . For generic slow-roll models, these 
terms dominate the primary non-Gaussianity. 



4.2.2 Numerical results of the primary bispectrum 

We calculate the primary CMB bispectrum (Eqs.( ^T^ )-( [4.22 )) numerically as follows. We compute 



the full radiation transfer function, gTi(k), with the CMBFAST code (Seljak and Zaldarriaga, 1996), 
assuming a single power-law spectrum, P$(k) oc k n ~ 4 , for the primordial curvature fluctuations. 
After doing the integration over k (Eqs.( |4T20[ ) and ( 4.21] )) with the same algorithm as of CMBFAST, 
we do the integration over r (Eq.flOg)), r = c(r - r), where r is the conformal time, tq is the 
present-day value. In our model, ctq = 11.8 Gpc, and the decoupling occurs at er* = 235 Mpc at 
which the differential visibility has a maximum. Our ctq includes radiation effects on the expansion 
of the universe; otherwise, ctq = 12.0 Gpc. Since the most of the primary signal is generated at 
t#, we choose the r integration boundary as c(tq — 2r*) < r < c(tq — O.lr*). We use a step-size 
of O.lcr*, as we have found that a step size of O.OIct* gives very similar results. As a cosmological 
model, we use the scale-invariant standard CDM model with J7 m = 1, £7a = 0, Ob = 0.05, h = 0.5, 
and n=l, and with the power spectrum, P$(k), normalized to COBE ( [Bunn and White, 1997 ). 
Although this model is almost excluded by current observations, it is still useful to depict basic 
effects of the transfer function on the bispectrum. 



Figure |4j] shows bf(r) (Eq. fl4.20| )) and &f L (r) (Eq. ([4.21 )) for several different values of 



We find that bf(r) and Q look very similar to each other in shape and amplitude at I ^ 100, 
although the amplitude in the Sachs-Wolfe regime is different by a factor of —3. This is because 
C\ oc P$(k)g^[(k), while bf(r) oc P$(k)gTi(k), where g^i = —1/3. We also find that b\{r) has a 
good phase coherence over wide range of r, while the phase of bf L (r) in the high-/ regime oscillates 
rapidly as a function of r. This strongly damps the integrated result (Eq. (|4.19|) ) in the high-/ 
regime. The main difference between C\ and b\{f) is that b\{r) changes the sign, while Cj does not. 

j^L ~ 10~ 10 Mpc -3 . The most signal 



Looking at figure |4.1| , we find I b\ ~ 2 x 10 and bf f } 
coming from the decoupling, the volume element at t* is r^Ar 



(10 4 ) 2 x 10 2 Mpc 3 ; thus, we 



estimate an order of magnitude of the primary reduced bispectrum (Eq. (|4.22|) ) as 



2r 2 Ar, (l 2 b\f h f^ x 3 



r 4 x2x 10 -17 /] 



NL- 



(4.23) 



NL f-1 
NL 



6(r — r*) (see Eq.( 4.26 )), r 2 Ar*6^ L / N L ~ 1- This rough estimate agrees with 



Since b? L f } 

the numerical result below (figure 

Figure shows the integrated bispectrum (Eq. (|4.1S| )) divided by the Gaunt integral, 
which is the reduced bispectrum, bf^™ Ty . While the bispectrum is a 3-d function, we show different 



mim2m3 
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Figure 4.1: Components of Primary CMB Bispectrum 

This figure shows bf(r) (Eq.(|4~20|)) and bf L (r) (Eq.(p~2l[)), the two terms in our calcula- 
tion of the primary CMB angular bispectrum, as a function of r. Various lines in the top 
panel show [l(l + l)6^(r)/2?r] x 10 10 , where r = c (r - r), at r = 0.4,0.6,0.8, 1.0, 1.2, 1.4, 
and 1.6 x r* (decoupling time); [&J* l (0/nl] x are snown m ^ ne bottom panel. To 
is the present-day conformal time. Note that ctq — 11.8 Gpc, and ct» = 235 Mpc 
in our cosmological model chosen here. The thickest solid line in the top panel is the 
CMB angular power spectrum, [1(1 + l)Ci/2ir] x 10 10 . Ci is shown for comparison. 



4.2. PRIMARY BISPECTRUM AND SKEWNESS 



63 



1-d slices of the bispectrum in this figure. We plot 

h(h + l)h(h + 1) (af^AmAms) to 2 " 3 )" 1 
clS cl function of I3 in the top panel, while we plot 

h(h + l)h(h + 1) (al imi almA L m3 ) te 2 " 3 )" 1 

in the bottom panel. We have multiplied each b\(r) which contains P$(k) by 1(1 + l)/(27r) so 
that the Sachs- Wolfe plateau at I3 10 is easily seen. We have chosen l\ and I2 so as (ii, Z2) = 
(9, 11), (99, 101), (199, 201), and (499,501). We find that the (h,l 2 ) = (199,201) mode, the first 
acoustic peak mode, has the largest signal in this family of parameters. The top panel has a 
prominent first acoustic peak, and strongly damped oscillations in the high-/ regime; the bottom 
panel also has a first peak, but damps more slowly. Typical amplitude of the reduced bispectrum 
is ^ 4 ^" mary /NL ~ 10 -17 , which agrees with an order of magnitude estimate (Eq. Q4.23 )). 



Our formula (Eq.( |4.22| )) and numerical results agree with Gangui et al. (1994| ) in the Sachs- 
Wolfe regime, where g^i(k) « — ji(kr*)/3, and 

€T y « -6/nl (cf^C™ + CTCT + Cf 2 W Cf 3 W ) (SW approximation). (4.24) 

Each term is in the same order as equation ( f4.22| ). Here, Cf w is the CMB angular power spectrum 
in the Sachs-Wolfe approximation, 

O POO 

Cf w = — / k 2 dkP^(k)jf(kn). (4.25) 

97T Jo 

In deriving equation ( [4.24 ) from ( f4.22| ), we have approximated bf L (r) (Eq.( 4.2l| )) with 



bf L (r) ™(-ff)~ jf ^dkMkrMkr) = -^r; 2 5(r - r m ). (4.26) 

The Sachs- Wolfe approximation (Eq.( f4.24| )) is valid only when 1%, I2, and I3 are all smaller 
than ~ 10, for which |Gangui et al. (1994|) gives ~ —6 x 10 -20 in figure 4^. We stress again 
that the Sachs-Wolfe approximation gives a qualitatively different result from our full calculation 
(Eq.( [4.22| )) at U ^ 10. The full bispectrum does change the sign; the approximation never changes 
the sign because of using Cf w . The acoustic oscillation and the sign change are actually great 
advantages, when we try to separate the primary bispectrum from various secondary bispectra. We 
will study this point later. 

4.2.3 Primary skewness 

The skewness, S3, given by 
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Figure 4.2: Primary CMB Bispectrum 
The primary angular bispectrum (Eq.( 4.19| )), divided by the Gaunt integral, G^l™ 2 ™ 3 
(Eq.(fO|)). The bispectrum is plotted as a function of l 3 for (h,h) =(9,11), (99,101), 
(199,201), and (499,501). Each panel plots a different 1-dimensional slice of the bispec- 
trum. The top panel plots l 2 (h + l)h(h + 1) (< h m ArnA 3 rn 3 ) /nl Km^V 
while the bottom panel plots h{h + 1)Z 2 (Z 2 + 1) (a^o^ogj^) f£ Kffi"*)' 1 l^f- 
Note that we have multiplied the bispectrum in each panel by a factor of 10 19 . 
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is the simplest statistic characterizing non-Gaussianity. We expand S3 in terms of Bi x i 2 i 3 (Eq.Q4.3j)), 



1 (2k + 1 ) (2Z 2 + 1) (2Z 3 + 1 



4"7T 



2vr 2 ^ 
xKl 2 l z W h Wi 2 W h , 



An 




B hhh W h W h W h 



(4.28) 



where W\ is the experimental window function. We have used equation (4J3) to replace -B; 1 / 2 z 3 by 
the reduced bispectrum, b\ x i 2 i z , in the last equality. Since I = and 1 modes are not observable, 



we have excluded them from the summation. Throughout this chapter, we consider a single-beam 
window function, W t = e -^+i)/(2^), where a h = FWHM/ v / 81n2. Since ( ^ 2 ^ J & Wa / 3 is 



e -i(/+i)/(2o- b ) > w here a h = FWHM/\/81n2. Since 
symmetric under permutation of indices, we change the way of summation as 



E 

2<hhh 



6 E 

2<h</ 2 <^3 



(4.29) 



This reduces the number of summations by a factor of ~ 6. We will use this convention henceforth. 

The top panel of figure plots S3(< I3), which is S3 summed up to a certain Z3, for FWHM 
beam sizes of 7°, 13', and 5'.5. These values correspond to COBE, MAP, and Planck beam sizes, 



respectively. Figure |4.3| also plots the infinitesimally thin beam case. We find that MAP, Planck, 
and the ideal experiments measure very similar S3 to one another, despite the fact that Planck and 
the ideal experiments can use much more number of modes than MAP. The reason is as follows. 



Looking at equation ( |4.28[) , one finds that S3 is a linear integral of bij,^ over thus, integrating 
oscillations in bf^™ Iy around zero (see figure |4.2[ ) damps the non-Gaussian signal on small angular 
scales, / ^ 300. Since the Sachs- Wolfe effect, no oscillation, dominates the COBE scale anisotropy, 
the cancellation on the COBE scales affects S3 less significantly than on the AL4Pand Planck scales. 
Planck suffers from severe cancellation in small angular scales: Planck and the ideal experiments 
measure only the same amount of S3 as MAP does. As a result, measured S3 almost saturates at 
the MAP resolution scale, / ~ 500. 

We conclude this section by noting that when we can calculate the expected form of the bispec- 
trum, then it becomes a "matched filter" for detecting the non-Gaussianity in data, and thus much 
more powerful tool than the skewness in which the information is lost through the coarse-graining. 



4.3 Secondary Sources of the CMB Bispectrum 

Even if the CMB bispectrum were significantly detected in the CMB map, the origin would not 
necessarily be primordial, but rather would be various secondary sources such as the Sunyaev- 
Zel'dovich (SZ) effect ( p^el'dovich and Sunyaev, 1969| ), the weak lensing effect, and so on, or fore- 
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Figure 4.3: Primary Skewness 

The top panel shows the primary CMB skewness (Eq.( 4.2S| )) summed up to a cer- 
tain l 3 , -S 3 (< h)f^l x 10 15 . The bottom panel shows the error of S 3 (Eq.( [l.6:j| )) 
summed up to I3, crs 3 (< ^3) x 10 15 . The solid line represents the zero-noise ideal 
experiment, while the dotted lines show COBE, MAP, and Planck experiments. 
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ground sources such as extragalactic radio sources. To isolate the primordial origin from the others, 
we have to know the accurate form of bispectra produced by secondary and foreground sources. 



4.3.1 Coupling between the weak lensing and the Sunyaev— Zel'dovich effects 

The coupling between the SZ effect and the weak lensing effect produces an observable effect in the 
bispectrum ( Goldberg and Spergel, 1999| ; Cooray and Hu, 20001 ) . We expand the CMB temperature 



field including the SZ and the lensing effect as 



AT (A) AT P (n + VG(n)) + AT sz (n) 



T T T 

AT p (n) /AT P (A)\ ^, an AT sz (n) , , 

~ -^^+Vf T { J • V9(n) + ^-L, (4.30) 

where P denotes the primary anisotropy, 0(n) is the lensing potential, 

G(n) = -2 / * dr^— -$(r, nr), (4.31) 
Jo rr* 

and SZ denotes the SZ effect, 

AT sz (n) , s . . 
= V(n)j v , (4.32) 

where j v is a spectral function of the SZ effect (Zel'dovich and Sunyaev, 1969). y(h) is the Compton 



y-parameter given by 

„ W = w /±Sj£*2.-» (r) , (4.33) 

» = ^Sg^ - « >< -Ve' (^) (m tih) ^ 



where 



T p = p gas T e /p gas is the electron temperature weighted by the gas mass density, the overline denotes 
the volume average, and the subscript means the present epoch. We adopt /x^ 1 = 0.88, where 
fi^ 1 = n e j (pgas/rrip) is the number of electrons per proton mass in the fully ionized medium. Other 
quantities have their usual meanings. 

Transforming equation ( 4.30| ) into harmonic space, we obtain 



aim = af m + X! X! (-l) m 5 u ™ 

I'm' l"m" 



a Vm i^ V i m ii + a lm 



P i \ 1 \ / -i\m+m'+m"/->—mm'm" 

a im + 2^ 2^ v _i J ywi" 

I'm' l"m" 

J(l> + !)_/(/ + 1) +/"(/" + !) ( sz 



a Z'-m'©r"-m" + a imJ (4.35) 
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Figure 4.4: SZ-lensing Coupling 

A schematic view of the SZ-lensing coupling bispectrum. One of the three CMB photons, 
which are decoupled at the last scattering surface (LSS), penetrates through a SZ cluster, 
changing its temperature, and coming toward us. As the other two photons pass through 
near the SZ cluster, they are deflected by the lensing effect, changing their propagation 
directions, and coming toward us. As a result, the three photons are correlated, generating 
three-point correlation, the bispectrum. 



where C/™ 1 ™ 2 " 13 is the Gaunt integral (Eq.([D|)). Substituting equation ( [4.35|) into flip, and usin 



hhfe 

the identity, Gi^i~ ma ~ m3 = G™lJ!^ m3 , we obtain the bispectrum, 



hfofa 



y hhh 



+5 permutations] 



sz 



rr pi* n sz 



in :i 



(4.36) 



The form of equation ( |4.5| ) is confirmed; the reduced bispectrum fr^/~|g ns includes the terms in the 
square bracket. 

While equation ( 4.36| ) is complicated, we can understand the physical effect producing the SZ- 



lensing bispectrum intuitively. Figure 4.4 shows how the SZ-lensing coupling produces the three- 
point correlation. Suppose that there are three CMB photons decoupled at the last scattering 
surface (LSS), and one of these photons penetrates through a SZ cluster between LSS and us; the 
energy of the photon changes because of the SZ effect. When the other two photons pass through 
near the SZ cluster, they are deflected by the gravitational lensing effect, changing their propagation 
directions, and coming toward us. What do we see after all? We see that the three CMB photons 
are correlated; we measure non-zero angular bispectrum. The cross-correlation strength between 
the SZ and lensing effects, ^0j*3 m3 a f 3 m 3 )' thus determines the bispectrum amplitude, as indicated 
by equation ( f4,36| ). 

Goldberg and Spergel (1999) have derived (®l m a^^j, assuming the linear pressure bias model 



( Persi et al., 1995| ), T p = T p b gas 5, and the mean temperature evolution, T p ~ T p q(1 + z) , for 
z < 2, which is roughly suggested by recent hydrodynamic simulations ( Pen and Ostriker, 1999| ; 
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Rcfrcgicr et al,, 2000| ; Springel et al,, 200l[ ), They have obtained 



Pi* n sz 



3ft, 



Pq, k 



r(z) 



(4.37) 



where D(z) is the linear growth factor. Simulations without non-gravitational heating (Refregier 
et al., 200q ; [Springel et al., 2001| ) suggest that T p0 ~ 0.2 - 0.4 keV and 6 gas ~ 5 - 10; analytic 
estimations give similar numbers ( [Refregier et al., 2000 ; Zhang and Pen, 200 In the pressure bias 
model, free parameters except cosmological parameters are T p q and 6 gas ; however, both actually 
depend upon cosmological models ( Refregier et al., 2000]) . Since Z 3 (@/* m a/ S m) ~2x 10 _10 j I/ T /9 o&gas 



QGoldberg and Spergel, 1999j ; |Cooray and Hu, 2000|) and l 2 Cf ~ 6 x 10" 



10 



■ ~ 2 x 10 
we have 



isz— lens 
b lll 



e; 



sz\ 

lm a lm / 



x 5/2 



r 3 x 3 x 10- 19 j u T pQ 



(a mi 



where T p q is in units of 1 keV, and bi x i 2 i 3 
(Eq.(fO|)). Comparing this with equation ( 4.23j ), we obtain 

/nl 



m 2 m 3 \ ' 
hh ) 



^primary 

lsz— lens 
°lll 



r 1 X 10 



v ju^ p0 b gas 



(4.38) 

is the reduced bispectrum 
(4.39) 



This estimate suggests that the SZ-lensing bispectrum overwhelms the primary bispectrum in small 
angular scales. This is why we have to separate the primary from the SZ-lensing effect. 

While the pressure bias model gives a rough estimate of the SZ power spectrum, more accurate 
predictions exist. Several authors have predicted the SZ power spectrum analytically using the 
Press-Schechter approach ( pole and Kaiser, 1988 ; Makino and Suto, 1993; [Atrio-Barandela and 
Miickct, 199S; Komatsu and Kitayama, 1999; Cooray, 2000| ) or the hyper-extended perturbation 
theory (Zhang and Pen, 2001 ). The predictions agree with hydrodynamic simulations well (Refregier 
|et al., 2000 ; Scljak et al., 2001; Springel et al., 2001; Rcfrcgicr and Teyssier, 2000| ; Komatsu et al. 



2001a). While a big uncertainty in the predictions has been lying in a phenomenological model 
which describes the SZ surface brightness profile of halos, Komatsu and Seljak (2001 ) have proposed 
universal gas and temperature profiles, and predicted the SZ profile relying on the physical basis. 
Hence, using those universal profiles, one can improve the analytic prediction for the SZ power 
spectrum. We expect that the universal profiles describe the SZ profile in the average sense; on the 
individual halo-to-halo basis, there could be significant deviation from the universal profile, owing 
to substructures in halos ( [Komatsu et al., 2001c ). 



4.3.2 Extragalactic radio and infrared sources 

The bispectrum from extragalactic radio and infrared sources whose fluxes, F, are smaller than a 
certain detection threshold, -Fd, is simple to estimate, when we assume the Poisson distribution. 
Toffolatti et al. (199§| ) have shown that the Poisson distribution is good approximation. The Poisson 
distribution has the white noise spectrum; thus, the reduced bispectrum (Eq.( |4.5| )) is constant, 
tfihh = ^ PS = cons ^ an t) an d we obtain 

Ki\z m3 =Gwr 3bP8 > ( 4 - 4 °) 
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where 

« F d ) = 5 3 (z) J"" dFF*^ = 9 "(x)^^Fl (4.41) 

Here, dn/dF is the differential source count per unit solid angle, and n(> Fd) = Jp 3 dF(dn/dF). 
We have assumed a power-law count, dn/dF oc F~@~ 1 , for (3 < 2. The other symbols mean 
x = hu/k B T ~ (V/56.80 GHz)(T/2.726 K) -1 , and 

, (lie) 2 /sinhx/2\ 2 1 / T \ ~ 3 /sinhx/2\ 2 

9( - X)= Jk^Tf\ : x* ) ~ 67.55 MJy sr^ 1 \2J26KJ ) ' ( ' 



Using the Poisson angular power spectrum, C ps , given by 

,dn__ 2 1 
dF =9 [ -~' 2 -[3 



C^« F d ) ee g\x) I ^ dFF 2 ^ = g\x) r ^^Fl (4.43) 



o 



we can rewrite 6 ps into a different form, 



b ps (< F d ) = (2 3 _^ /2 [n(> Fa)]" 1 / 2 [C ps (< F d )] 3/2 . (4.44) 



Tottolatti et al. (1998) have estimated n(> F d ) ~ 300 sr" 1 for F d ~ 0.2 Jy at 217 GHz. This 



Fi corresponds to 5a detection threshold for the Planck experiment at 217 GHz. Refregier et al 



(2000D have extrapolated their estimation to 94 GHz, finding n(> Fi) ~ 7 sr 1 for Fj ~ 2 Jy, 



which corresponds to MAP 5a threshold. These values yield 



C ps (90 GHz, < 2 Jy) ~ 2 x HT 16 , (4.45) 
C ps (217 GHz, < 0.2 Jy) ~ 1 x 10" 17 . (4.46) 

Thus, rough estimates for 6 ps are 

6 ps (90 GHz,< 2 Jy) ~ 2 x 10~ 25 , (4.47) 
6 ps (217 GHz, < 0.2 Jy) ~ 5 x 10" 28 . (4.48) 

While we have assumed the Euclidean source count ((3 = 3/2) for definiteness, this assumption 
does not affect an order of magnitude estimate here. 

As the primary reduced bispectrum is oc F 4 (Eq.(4.23)), and the SZ-lensing reduced bispectrum 
IS oc F 3 (Eq.flP8D), the point -source bispectrum rapidly becomes to dominate the total bispectrum 
in small angular scales: 

~ r 4 xl0 7 f L /NL _ 9 , V (4.49) 



5ps V^ ps /10- 25 /' 

isz-lens / • Tp 7 \ 

2m r 3 x io 6 i vI f° bg f. . (4.50) 

For example, the point sources overwhelm the SZ-lensing bispectrum measured by M^4Pat I ^ 100. 
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What the SZ-lensing bispectrum and the point-source bispectrum look like? Figure [D^ plots the 
primary, the SZ-lensing, and the point-source reduced bispecta for the equilateral configurations, 
I = l\ = li = I3. We have plotted l 2 (l + l) 2 bm/(2ir) 2 . We find that these bispecra are very different 
from each other in shape on small angular scales. It thus suggests that we can separate these three 
contributions on the basis of shape difference. We study this point in the next section. 



4.4 Measuring Bispectra 



4.4.1 Fisher matrix 



In this section, we study how well we can measure the primary bispectrum, and how well we can 
separate it from the secondary bispectra. Suppose that we fit the observed bispectrum, , 
by theoretically calculated bispectra, which include both the primary and secondary sources. We 
minimize x 2 defined by 



X' 



E 

2<h<h<l 3 



{ B hltl 3 Ei A i B hhh 



(4.51) 



(1/2/3 



where % denotes a component such as the primary, the SZ and lensing effects, extragalactic sources, 
and so on. We have removed unobservable modes, I = and 1. 

As we have shown in chapter ||, the variance of the bispectrum, afi is the six-point function 
of a\ m ( Luo, 1994 ; Heavens, 1998 ). When non-Gaussianity is weak, we calculate it as ( (Spcrgcl and 
Goldberg, 19991 ; |Gangui and Martin, 2000| ) 



°hhh = \ B hhh 



{ B hhh) 



(4.52) 



where A/ 1 ; 2 ; 3 takes values 1, 2, or 6 when all Vs are different, two are same, or all are same, 
respectively. C\ = C\ + Cf 1 is the total CMB angular power spectrum, which includes the power 



spectrum of the detector noise, C ; N . We calculate Q N analytically following Knox (1995 ) with the 
noise characteristics of relevant experiments. We do not include C\ from secondary sources, as they 
are subdominant compared with the primary C\ and Cf for relevant experiments. Including C\ 
from extragalactic sources (Eqs.(E~45p or flP6| )) chan ges our results less than 10%. 
Taking dx 2 /dAi = 0, we obtain normal equation, 



E 



B (i) B U) 



2<h<h<h 



hhh 



A: 



E 

2<h<h<h 



Dobs r(*) 
D hhh D hhh 



a 



(4.53) 



hhh 



We then define the Fisher matrix, as 



1.1 



E 



hhh hhh 



2</i</ 2 </3 

2 

7T 



a 



h hh 



2<h<h<h 



h + 




hhh hhh 



hhh 



(4.54) 
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Figure 4.5: Equilateral Reduced Bispectra 
Comparison between the primary (solid line), the SZ-lensing (dashed line), and the point- 
source (dotted line) reduced bispectra for the equilateral configurations, I = l\ = fo = 
l 3 . We have plotted [l 2 {I + l) 2 b m / {2tt) 2 \ x 10 16 , which makes the Sachs-Wolfe plateau 
of the primary reduced bispectrum on large angular scales, I ^ 10, easily seen. 
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where we have used equation ( |4.6| ) to replace Bi ± i 2 i 3 by the reduced bispectrum, fe^jjjg (see Eq.([4.5D 
for definition). Since the covariance matrix of A{ is F^j , we define the signal-to-noise ratio, (S/N)i, 
for a component i, the correlation coefficient, rj-,-, between different components i and j, and the 
degradation parameter, di, of (S/N)i due to rjj, as 



f 



-l 



FZ 



di 



F^F- 1 



F ■ F~ 



(4.55) 

(4.56) 
(4.57) 



Note that does not depend upon the amplitude of bispectra, but the shape. We have defined 
di so as di = 1 for zero degradation, while di > 1 for degraded (S/N)i. Spergel and Goldberg 



(1999 ) and Cooray and Hu (200C| ) have considered the diagonal component of F~j . We study all 
the components to study the separatability between various bispectra. 

We estimate an order of magnitude of S/N as a function of a certain angular resolution, I, as 

' I I I \ 2 
/ 



follows. Since the number of modes contributing to S/N increases as I 3 / 2 , and I 3 
0.36 x I, we estimate (S/N)i ~ (Fa) 1 / 2 as 



_L /3 /2 /3 /2 

3vr 




v lll 



fib® x 4 x 10 12 , 



(4.58) 



-10 



where we have used l 2 Ci ~ 6 x 10 

Table pQl tabulates Fij, while table [Q| tabulates F-" 1 ; table |Q| tabulates (S/N)i, while table [O 
tabulates di in the diagonal, and in the off-diagonal parts. 



4.4.2 Measuring primary bispectrum 



Figure 4J3 shows the signal-to-noise ratio, S/N. The top panel shows the differential S/N for 
the primary bispectrum at Inl^ interval, [d(S/N) 2 /dln I3] ^ 2 f^, and the bottom panel shows the 
cumulative S/N, (S/N)(< ls)f^, which is S/N summed up to a certain Z 3 . We have computed 
the detector noise power spectrum, Cf, for COBE four-year map QBennett et al., 1996|) , MAP 90 
GHz channel, and Planck 217 GHz channel, and assumed full sky coverage. Figure [O] also shows 
the ideal experiment with no noise: Cf 1 = 0. Both [d(S/N) 2 /din I3] and (S/N)(< Z3) increase 
monotonically with Z3, roughly cx Z3, up to 1% ~ 2000 for the ideal experiment. 

Beyond I3 ~ 2000, an enhancement of the damping tail in C/ because of the weak lensing 

1 /2 

effect dSeljak, 1996|) stops [d(S/N) 2 /dlnl 3 ] 1 , and hence (S/N)(< l 3 ), increasing. This leads to 
an important constraint on the observation; even for the ideal noise-free, infinitesimally thin beam 
experiment, there is an upper limit on the value of S/N ^ 0.3/nl- For a given realistic experiment, 
[d(S/N) 2 /dlnl 3 ] 1/2 has a maximum at a scale near the beam size. 
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Table 4.1: Fisher Matrix 



Fisher matrix, F^ (Eq.(4.54)). i denotes a component in the first row; j denotes a component 



in the first column. T p q is in units of 1 keV, b^l = 6 ps /10 



-2", 



and &27 



6P S /10 



-27 



CORE 


primary 


SZ-lensing 


point sources 


primary 


4.2 x 10- 6 / 2 L 


-4.0 x lO" 7 h L juT p0 b sas 


-1.0 x 10" 9 /nlC 


SZ-lensing 




1.3 x lO" 7 (j u T p0 b gas ) 2 


3.1 x 10- 10 j v Tpob gas bf b 


point sources 






1.1 x 10- 12 (^) 2 


MAP 


primary 


3.4 x 10- 3 / 2 L 


2.6 x 10~ 3 f NL jvT p0 b gas 


2.4 x 10" 3 f NL bf 5 


SZ-lensing 




0.14 (>T p0 6 gas ) 2 


0.31 jvT p0 b gas bf 5 


point sources 






5.6 (fog) 2 


Planck 


primary 


3.8 x 10- 2 / 2 L 


7.2 x 10" 2 f NL j»T p0 b gas 


1.6 x 10" 2 fa L bf 7 


SZ-lensing 




39 (j„T p0 b gas ) 2 


5-7 juT p ob gas b2j 


point sources 






2.7 x 10 3 (b%) 2 



Table 4.2: Inverted Fisher Matrix 



Inverted Fisher matrix, F- 1 . The meaning of the symbols is the same as in table 4.1 



CORE 


primary 


SZ-lensing 


point sources 


primary 
SZ-lensing 
point sources 


3.5 x 10 5 / NL 2 


1.1 x 10 6 (hduTpob^)- 1 
3.1 x 10 7 {jvT p0 b gas )- 2 


1.3 x 10 7 (/nlC)" 1 

-7.8 X 10 9 (^Tpoftgas^s)" 1 

3.1 x 10 12 (6^)- 2 


MAP 


primary 
SZ-lensing 
point sources 


3.0 x 10 2 / NL 2 


-6.1 (fahjvT p0 b gas ) 1 
8.4 UuT p0 b gas )- 2 


0.21 (/NL^r 1 
-0.46 (jvTpob^bH)- 1 
0.21 (6^)~ 2 


Planck 


primary 
SZ-lensing 
point sources 


26 /nl 


-4.9 x lO" 2 {hWuTpob^y 1 
2.6 x 10~ 2 {j u T pQ b gas )- 2 


-5.7 x lO" 5 (hLb^)- 1 
-5.4 x lO" 5 (juTpob^b^y 1 
3.7 x 10- 4 (bf 7 )- 2 
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Table 4.3: Signal-to-noise Ratio 



Signal-to-noise ratio, (S/N)i (Eq.(4.55)), of detecting the bispectrum. i denotes a compo- 



nent in the first row. The meaning of the symbols is the same as in table 4.1 





primary 


SZ-lensing 


point sources 


COBE 
MAP 
Planck 


1.7 x 10" 3 / NL 

5.8 x 1(T 2 / NL 
0.19 /nl 


1.8 x 10" 4 \UT p0 b gas 
0.34 \j u \T p0 b gas 

6.2 \ju\Tpobgaa 


5.7 x 10" 7 bf 5 
2.2^ " 
52 ftg 



Table 4.4: Signal Degradation and Correlation Matrix 
Signal degradation parameter, di (Eq.(4.57)), and correlation coefficient, ry (Eq.(4.56)) 



matrix, i denotes a component in the first row; j denotes a component in the first column. 
di for i = j, while for i ^ j . 



COBE 


primary 


SZ-lensing 


point sources 


primary 


1.46 


0.33 sgn(>) 


1.6 x 10" 2 


SZ-lensing 




3.89 


-0.79 sgn(>) 


point sources 






3.45 


MAP 


primary 


1.01 


-0.12 8 gn(j„) 


2.7 x 10~ 2 


SZ-lensing 




1.16 


-0.35 sgn(>) 


point sources 






1.14 


Planck 


primary 


1.00 


-5.9 x 10~ 2 sgn(>) 


-5.8 x 10~ 4 


SZ-lensing 




1.00 


-1.8 x 10~ 2 sgn(>) 


point sources 






1.00 
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Figure 4.6: Signal-to-noise Ratio 

The predictions of the signal-to-noise ratio, S/N, for COBE, MAP, and Planck exper- 
iments (see Eq.( 4.55| )). The differential S/N at ln?3 interval is shown in the upper 
panel, while the cumulative S/N up to a certain I3 is shown in the bottom panel. 
Both are in units of /nl- Solid line represents the zero-noise ideal experiment, while 
dotted lines show the realistic experiments mentioned above. The total (S/N)f^ are 
1.7 x 10~ 3 , 5.8 x f0~ 2 , and 0.19 for COBE, MAP, and Planck experiments, respectively. 
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For COBE, MAP and Planck experiments, the total (S/N)fal are 1.7 x 1CT 3 , 5.8 x 1(T 2 , and 
0.19, respectively (see table 4.3). To obtain S/N > 1, we need /nl > 600, 20, and 5, while the 
ideal experiment requires /nl > 3 (see table [h5|). We can also roughly obtain these values by 
substituting equation ( 4.23|) into (|4.5S| ), 



(4.59) 



~ / x 10~ 4 /nl. 

primary 

The degradation parameters, d pr i mary , are 1.46, 1.01, and 1.00 for COBE, MAP, and Planck 



experiments, respectively (see table 4^4), suggesting that MAP and Planck experiments will sepa- 
rate the primary bispectrum from the others with 1% or better accuracy; however, COBE cannot 
discriminate between them very well, as the primary and the secondary sources change monotoni- 
cally on the COBE angular scales. On the MAP and Planck scales, the primary bispectrum starts 
oscillating around zero, being well separated in shape from the secondaries that do not oscillate. 
This is good news for the forthcoming high angular resolution CMB experiments. 



4.4.3 Measuring secondary bispectra 

Signal-to-noise ratios of detecting the SZ-lensing bispectrum, (S/N) sz ^\ cns , in units of \j u \ T p ob gas 
are 1.8 x 10 -4 , 0.34, and 6.2 for COBE, MAP, and Planck experiments, respectively (see table JO), 
where T p o is in units of 1 keV. Using equations ( 4.58 ) and ( 4.38 ), we estimate (S/N) sz _\ ens roughly 

(-77) -l 2 x 10- 6 |j,|r p0 & g a S . (4.60) 

V * * / sz— lens 

Hence, (S/N) sz _[ ens increases with the angular resolution more rapidly than the primary bispectrum 
(see Eq.(p39p). Since \j u \ T p ob gSuS should be of order unity, COBE and MAP cannot detect the SZ- 
lensing bispectrum; however, Planck is sensitive enough to detect, depending on the frequency, i.e., 
a value of j u . For example, 217 GHz is insensitive to the SZ effect as j u ~ 0, while j u = —2 in the 
Rayleigh-Jeans regime. 

The degradation parameters, <i sz _i cns , are 3.89, 1.16, and 1.00 for COBE, MAP, and Planck 
experiments, respectively (see table 4.4); thus, Planck will separate the SZ-lensing bispectrum 
from the other effects. Note that (S/N) sz _\ ens values must be an order of magnitude estimation, 
for our cosmological model is the COBE normalized SCDM that yields as = 1.2, which is a factor 
of 2 greater than the cluster normalization for Q m = 1, and 20% greater than the normalization 
for Q m = 0.3 ( Kitayama and Suto, 199"^ ). Hence, this factor tends to overestimate ^@;* m op^^) 
(Eq.( [4.37 )) by a factor of less than 10; on the other hand, using the linear P$(k) power spectrum 
rather than the non-linear power spectrum tends to underestimate the effect by a factor of less 
than 10 at Z ~ 3000 (Cooray and Hu, 2000). Yet, our main goal is to discriminate between the 
shapes of various bispectra, not to determine the amplitude, so that this factor does not affect our 
conclusion on the degradation parameters, <ij. 

For the extragalactic radio and infrared sources, we estimate the signal-to-noise ratios as 
5.7 x 10- 7 (6 ps /10- 25 ), 2.2(&p s /10- 25 ), and 52(&p s /10- 27 ) for COBE, MAP, and Planck experi- 
ments, respectively (see table [4.3| ), and the degradation parameters, d ps , as 3.45, 1.14, and 1.00 
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(see table 4.4). Our estimate is consistent with Rcfregier et al. (2000| ). From equation ( f4,58| ), we 
find 

~ Z 5 x 10- 13 (J^) ; (4.61) 



N 



ps 



10 



-25 



thus, S/N of the point-source bispectrum increases very rapidly with the angular resolution. 

Although MAP cannot separate the Poisson bispectrum from the SZ-lensing bispectrum very 
well (see in table |4.4| ), the SZ-lensing bispectrum is too small to be measured by MAP anyway. 
Planck will do an excellent job on separating all kinds of bispectra, at least including the primary 
signal, SZ-lensing coupling, and extragalactic point sources, on the basis of the shape difference. 



4.4.4 Measuring primary skewness 

For the skewness, we define S/N as 



2 s? 



-) 

N a 



2 ' 

s 3 



(4.62) 



where the variance is dSrednicki, jggf ) 

- 1 dcosi 



4 3 



6 



-i 



[C(9)f 



(SsY 



6E 

2 i 



(2Zi + 1) (2l 2 + 1) (2Z 3 + 1) / h l 2 I 
(4tt) 3 \ 




xd 1 c h c h wf 1 w£Wf 3 . 



(4.63) 



In the last equality, we have used symmetry of the summed quantity with respect to indices 
(Eq. fl4.29| )), and removed unobservable modes, I = and 1. Typically as 3 ~ 10~ 15 , as as 3 ~ 
[C(0)P 

:< h) 



^ 10 15 , where C{9) is the temperature auto correlation function including noise. 
The bottom panel of figure 4.3 plots o"s 3 (< Z3), which is as 3 summed up to a certain Z3, 



for COBE, MAP, and Planck experiments as well as for the ideal experiment. Since CiWf 



Cie '(' +1 )°b + w 1 , where w 1 is the white noise power spectrum of the detector noise (Knox 



1995 ), w 1 keeps 0"s 3 (< Z3) slightly increasing with £3 beyond the experimental angular resolution 



scale, / ~ <7 b . In contrast, Ss(< I3) becomes constant beyond / ~ a h (see the top panel of fig- 



ure [4.3| ). As a result, S/N starts slightly decreasing beyond the resolution. We use the maximum 
S/N for calculating the minimum value of /nl above which the primary S3 is detectable; we find 
that /nl > 800, 80, 70, and 60 for COBE, MAP, Planck, and the ideal experiments, respectively, 
assuming full sky coverage. 

These /nl values are systematically larger than those for detecting B^i^ by a factor of 1.3, 4, 
14, and 20, respectively (see table |4~5| ) . The higher the angular resolution is, the less sensitive the 
primary S3 is to non-Gaussianity than i?/ 1 ; 2 ; 3 . This is because the cancellation effect on smaller 
angular scales because of the oscillation of damps S3. 



4.5. DISCUSSION AND CONCLUSIONS 



79 



Table 4.5: Detection Limit for the Non-linear Coupling Parameter 
The minimum non- linear coupling parameter, /nl, needed for detecting the primary non- 
Gaussianity by the bispectrum or the skewness with the signal-to-noise ratio greater than 
1. These estimates include the effects of cosmic variance, detector noise, and foreground 
sources. 

Experiments /nl (Bispectrum) /nl (Skewness) 
COBE 600 800 

MAP 20 80 

Planck 5 70 

Ideal 3 60 



4.5 Discussion and Conclusions 

In this chapter, using the full radiation transfer function, we have numerically computed the primary 
cosmic microwave background bispectrum (Eq.( |4.19|) ) and skewness (Eq. (|4.28| )) down to arcminutes 
angular scales. As the primary bispectrum oscillates around zero (figure |4~2^ ) , the primary skewness 
saturates at the MAP angular resolution scale, I ~ 500 (figure |4~3| ) . We have introduced the reduced 
bispectrum, bi x i 2 i 3 , defined by equation fl4.5| ), and found that this quantity is more useful to describe 
physical properties of the bispectrum than the angular averaged bispectrum, (Eq.( |4.3D ). 

Figure [O] compares the expected signal-to-noise ratio of detecting the primary non-Gaussianity 
based on the bispectrum (Eq.( 4.55 )) with that based on the skewness (Eq.( 4.62; )). It shows that 



the bispectrum is almost an order of magnitude more sensitive to the non-Gaussianity than the 
skewness. We conclude that when we can compute the predicted form of the bispectrum, it becomes 
a "matched filter" for detecting the non-Gaussianity in data, and thus much more powerful tool than 



the skewness. Table L5 summarizes the minimum /nl for detecting the primary non-Gaussianity 
using the bispectrum or the skewness for COBE, MAP, Planck, and the ideal experiments. This 
shows that even the ideal experiment needs /nl > 3 to detect the primary bispectrum. 

We have calculated the secondary bispectra from the coupling between the SZ effect and the 
weak lensing effect, and from the extragalactic radio and infrared sources. Only Planck will detect 
the SZ-lensing bispectrum, while both MAP and Planck will detect the extragalactic point-source 
bispectrum (table |4.3|) . 

We have also studied how well we can discriminate between the primary, the SZ-lensing cou- 
pling, and the extragalactic point-source bispectra. We have found that MAP and Planck will 
separate the primary from the other secondary sources with 1% or better accuracy. This con- 
clusion is due to the acoustic oscillation in the primary bispectrum that does not appear in the 
secondary bispectra. The SZ-lensing coupling and the extragalactic sources are well separately 
measured by Planck experiment, although COBE and MAP cannot discriminate between them 
(table |3). 

Our arguments on the ability to discriminate between various bispectra have been based upon 
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Figure 4.7: Bispectrum vs Skewness 

Comparison of the signal-to- noise ratio summed up to a certain I3, S/N(< for the 
bispectrum (top panel; Eq.( 4.55| )) and the skewness (bottom panel; Eq. (|4.62|) ). S/N(< I3) is 
in units of /nl- The dotted lines show COBE, MAP, and Planck experiments (dotted lines), 



while the solid line shows the ideal experiment. See table 4.5 for /nl to obtain S/N > 1 
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the shape difference, and thus have not taken into account the spectral difference in frequency 
space. [Tegmark and Efstathiou (1996 ) and Cooray et al. (2000) have shown that the multi-band 
observation is efficient for discriminating between the primary signal and the other foreground 
sources in the CMB power spectrum. Their scheme should be effective on the bispectrum as well, 
and will improve the accuracy of the foreground removal further. We thus expect that MAP and 
Planck will measure the primary bispectrum separately from the foregrounds. 

Simple slowly-rolling single- field inflation models predict /nl ~ 0(1O -2 ) ( [Balopek and BondJ 
1990] g99| ; IGangui et al., 1994D , while the second order perturbation theory predicts /nl ~ 0(1) 



(Pyne and Carroll, 1996); thus, significant detection of the primary bispectrum or the skewness with 
any experiments means that these inflation models need to be modified. According to our results, if 
the reported detection ( [Ferreira et al., 1998 ; Magueijo, 200C| ) of the bispectrum on the COBE map 
were cosmological in origin, then MAP and Planck would detect the primary bispectrum much more 
significantly. While Banday, Zaroubi and Gorski ( [Banday et al., 2000| ) have shown that the one of 
those detections ( Ferreira et al., 1998j ) is accounted for by the experimental systematic effects, the 
other flMaguei.jo, 2000|) is significant even after removing such the systematics. 

Although we have not discussed so far, spatial distribution of interstellar dust emissions is a 
potential source of microwave non-Gaussianity. While it is very hard to estimate the bispectrum 
analytically, we can use the dust template map compiled by Schlegel et al. (1998 ) to estimate 
the dust bispectrum. For example, we have found that the dimensionless skewness parameter, 
((AT) 3 ) / ((AT) 2 ) 3//2 , is as large as 51 on the template map. We have used the publicly available 
HEALPix-formatted ( Porski et al., 1998|) 100 //m map, which contains 12,582,912 pixels without 
sky-cut. The mean intensity in the map is 14.8 MJy sr _1 . Of course, this skewness is largely 
an overestimate for the real CMB measurement; we need to cut a fraction of the sky that con- 
tains the Galactic plane, and this will greatly reduce the non-Gaussianity. Nevertheless, residual 
non-Gaussianity is still a source of the microwave bispectrum, and has to be taken into account. 
Moreover, the form of the measured bispectrum on the dust map reflects the physics of interstellar 
dust, which is highly uncertain at present; thus, studying the interstellar dust bispectrum is a 
challenging field. 



Chapter 5 



Measurement of Bispectrum on the 
COBE DMR sky maps 



Several authors have attempted to measure non-Gaussianity in CMB using various statistical tech- 
niques (e.g., Kogut et al. 1996b); as yet no conclusive detection has been reported except for 
measurement of several modes of the normalized CMB bispectrum on the COBE Differential Mi- 
crowave Radiometer (DMR) sky maps ( [Ferreira et al., 1998| ; [Ma gueijo, 2000| ). The existence of 



non-Gaussianity in the DMR data is controversial. If the CMB sky were non-Gaussian, this would 
challenge our simplest inflationary model. 

The angular bispectrum, Bi t i 2 i 3 , is the harmonic transform of the three-point correlation func- 
tion. We carefully distinguish the normalized bispectrum, B\ x i 2 i z l (C^C^C^) 1 ' 2 , from the bispec- 
trum, B^i^. Ferreira et al. (1998|) have measured 9 equilateral (h = I2 = h) modes of the normal- 



1/2 

ized bispectrum, B^i^/ (C^C^C^) 1 , on the DMR map, claiming detection at h = h = h = 16. 



Their result has been under extensive efforts to confirm its significance and origin. Bromley and 



Tegmark (1999 ) claim that a few individual pixels in the DMR map are responsible for the most 



of the signal. Banday et al. (2000D have proposed an eclipse effect by the Earth against the COBE 



satellite as a possible source of the signal. Magueijo (2000| ) has measured other 8 inter-/ modes 



1 /2 

of the normalized bispectrum such as Bi_m + i/ (Ci_iCiCi+i) ' , and claims that scatter of the 



normalized bispectrum among 8 modes is too small to be consistent with Gaussian. Sandvik and 



Magueijo (200C| ) further report measurement of 24 other inter-Z modes for different lags in I, and 



conclude they are consistent with Gaussian. 

Thus, until now 41 modes of the normalized bispectrum have been measured on the DMR map. 
Here, we simply ask: "how many modes are available in the DMR map for the bispectrum?" The 
answer is 466, up to a maximum multipole of 20 that corresponds to the DMR beam size; thus, it is 
conceivable that the claimed detection of the normalized bispectrum at h = h = h = 16 would be 
explained by a statistical fluctuation, as 9 modes are expected to have statistical significance above 
98% out of 466 independent modes even if CMB is exactly Gaussian. In this chapter, we measure 
466 modes of the CMB bispectrum on the COBE DMR sky maps, testing the claimed detection 
of the bispectrum and non-Gaussianity. We take into account the covariance between these modes 
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due to the Galactic cut, which has not been done in the previous work. 

On the theoretical side, several predictions for the CMB bispectrum exist. Several authors QFalk 
|et al., 1993 ; Luo and Schramm, 1993 ; Gangui et al., 1994j ) have predicted the primary bispectrum 
(or equivalently three-point correlation function) on the DMR angular scales from slow-roll inflation 
models. In chapter |], we have extended the prediction down to arcminutes scales using the full 
radiation transfer function. 

In addition to the primary one, secondary sources in the low-redshift universe and foreground 
sources produce the bispectrum through their non-linearity. Luo and Schramm (1993j ) and Spergel 
and Goldberg (1999) have calculated the secondary bispectrum arising from non- linear evolution 
of gravitational potential; Goldberg and Spergel (1999| ) and pooray and Hu (200C| ) have calculated 
the one from the gravitational lensing effect coupled with various secondary anisotropy sources. 
In chapter ^j, we have calculated the foreground bispectrum from extragalactic radio and infrared 
point sources. While the bispectrum is not the best tool for detecting the signature of rare highly 
non-linear events, e.g., textures ( [Phillips and Kogut, 2001 ), it is sensitive to weakly non-linear 
effects. 

Having theoretical predictions is a great advantage in extracting physical information from 
measurement; one can fit a predicted bispectrum to the data so as to constrain parameters in a 
theory. Since the DMR beam size is large enough to minimize contribution from the secondary and 
the extragalactic foreground sources, the only relevant source would be the primary one. In this 
chapter, we fit a theoretical primary bispectrum (Komatsu and Spergel, 2001a) to the data. 

The Galactic plane contains strong microwave emissions from interstellar sources. The emissions 
are highly non-Gaussian, and distributed on fairly large angular scales. Unfortunately, predicting 
the CMB bispectrum from interstellar sources is very difficult; thus, we excise the galactic plane 
from the DMR data. We model the residual foreground bispectrum at high galactic latitude using 
foreground template maps. By simultaneously fitting the foreground bispectrum and the primary 
bispectrum to the DMR data for three different Galactic cuts, we quantify the importance of the 
interstellar emissions in our analysis. 

This chapter is organized as follows. In § [O], we define the angular bispectrum, and show how 
to compute it efficiently from observational data. In § 5.2, we study statistical properties of the 



bispectrum and the normalized bispectrum. We then measure the normalized bispectrum from 
the COBE DMR four-year sky maps (Bennett et al., 1996), testing Gaussianity of the DMR map. 
In § |5.3| , we fit predicted bispectra to the DMR data, constraining parameters in the predictions. 
The predictions include the primary bispectrum from inflation and the foreground bispectrum from 
interstellar Galactic emissions. Finally, § 5.4 concludes. 



5.1 Angular Bispectrum 

The CMB angular bispectrum consists of a product of three harmonic transforms of the CMB 
temperature field. For Gaussian fields, expectation value of the bispectrum is exactly zero. Given 
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statistical isotropy of the universe, the angular averaged bispectrum, B^^, is given by 

Bhhh = I I a hm 1 a l2m2 a hm- i i (5-1) 

Mm V mi ™ 2 m 3 J 

where the matrix denotes the Wigner-3j symbol, and harmonic coefficients, ai m , are given by 

a lm = f d 2 n ^P-Y? m (n) . (5.2) 

"'Slobs 1 

^obs denotes a solid angle of the observed sky. B^i^ satisfies the triangle condition, \k — lj\ < 
Ik < k + Zj for all permutations of indices, and parity invariance, h + h + h = even. 
We can rewrite equation Q5.1|) into a more useful form. Using the identity, 



-i 



h h h \ t h h h \ / (4vr)3 

mi m 2 m 3 J \ J \j (2Zi + 1) (2/ 2 + 1) (2Z 3 + 1) 

d 2 n 



/d n 
— ll imi (n)y i2m2 (n)yi 3m3 (n), (5.3) 



we rewrite equation (|5.1|) as 



5 iii 2 « 3 = q q q J e h (h)ei 2 (h)e l3 {h), (5.4) 

where the integral is not over rj b s , but over the whole sky; e;(n) already encapsulates the infor- 
mation of incomplete sky coverage through a\ m . Following chapter ||], we used the azimuthally 
averaged harmonic transform of the CMB temperature field, e/(n), 



. 21 + 

Similarly, we write the angular power spectrum, Q, as 



ez(n) = J ai m Y lm (h). (5.5) 



(',- l^efin). (5.6) 



Thus, e;(n) is a square-root of Ci at a given position of the sky. Equation (|5.4| ) is computationally 
efficient, as we can calculate e;(n) quickly with the spherical harmonic transform for a given I. 
Since the HEALPix pixels have the equal area flGorski et al., 1998|) , the average over the whole sky, 
/ d 2 n/(47r), is done by the sum over all pixels divided by the total number of pixels, iV -1 Y^i ■ 



5.2 Measurement of Bispectrum on the DMR Sky Maps 
5.2.1 The data 



We use the HEALPix-formatted (Gorski et al., 1998) COBE DMR four-year sky map, which con- 
tains 12,288 pixels in Galactic coordinate with a pixel size 1?83. We obtain the most sensitive sky 
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Table 5.1: Monopole and Dipole Subtraction 

The monopole, To, and the dipole, T\, anisotropics that have been subtracted from the DMR 
sky maps. We show the subtracted values for zero, 20°, 25°, and 30° cuts. The rightmost 
column shows the directions of the subtracted dipole in Galactic coordinate. 



beat | Tp [/jK] T\ [fiK] Ti/Ti (l,b) 

0° 1.40 63.1 (28?07, 2? 12) 

20° -70.3 26.2 (89P23, -4.°17) 

25° -72.4 26.9 (84P89, -5.°13) 

30° -73.4 28.6 (97P32, -6.°31) 



map to CMB by combining 53 GHz map with 90 GHz map, after coadding the channels A and B 
at each frequency. We do not subtract eclipse season time-ordered data; while Banday et al. (200C| ) 
ascribe the reported non-Gaussianity to this data, we will argue in this paper that the claimed 
detection of the normalized bispectrum at fa = fa = fa = 16 (Ferreira et al., 1998) can also be 
explained in terms of a statistical fluctuation. 

We reduce interstellar Galactic emissions by using three different Galactic cuts: 20°, 25°, and 
30° in Galactic latitude. Since we want to see how the different Galactic cuts affect the measured 
bispectrum, we use the three different Galactic cuts instead of the extended Galactic cut ( |Banday 
|et al., 1997 ), which is commonly used for analyzing the DMR sky maps. Then, we subtract the 
monopole and the dipole from each cut map, minimizing contaminations from these two multipoles 
to higher order multipoles through the mode-mode coupling. The coupling arises from incomplete 
sky coverage. This is very important to do, for the leakage of power from the monopole and the 
dipole to the higher order multipoles is rather big. We use the least-squares fit weighted by the 
pixel noise variance to measure the monopole and the dipole on each cut map. Table ^1] shows the 
amplitude of the subtracted monopole and the dipole for the different Galactic cuts. 

We measure the bispectrum, Bi t i 2 i 3 , on the DMR sky maps as follows. First, we measure a/ m 
using equation (|5.2j ). Then, we transform a\ m for —l<m<l into e;(n) through equation ( |5.5D , 



Finally, we obtain -B/ 1 ; 2 i 3 from equation (5.4), arranging fa, fa, and fa in order of fa < fa < fa, where 
the maximum fa is set to be 20. In total, we have 466 non-zero modes after taking into account 
\h — lj\ < fa < h + fa and fa + fa + fa = even. Measurement of 466 modes takes about 1 second of 
CPU time on a Pentium-Ill single processor personal computer. 



5.2.2 Monte— Carlo Simulations 

We use Monte-Carlo simulations to estimate the covariance matrix of the measured bispectrum. 
Our simulation includes (a) a Gaussian random realization of the primary CMB anisotropy field 
drawn from the CO-RE-normalized ACDM power spectrum, and (b) a Gaussian random realization 
of the instrumental noise drawn from diagonal terms of the COBE DMR noise covariance matrix 
( Lincweaver et al., 1994| ). For computational efficiency, we do not use off-diagonal terms as they 
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are smaller than 1% of the diagonal terms ( [Lineweaver et al., 1994] ) . 

We generate the input power spectrum, Cj, using the CMBFAST code ( [Seljak and Zaldarriagg 



1996J) with cosmological parameters fixed at O c dm = 0.25, fl\ = 0.7, Ob = 0.05, h = 0.7, and n = 1; 



the CMBFAST code uses the Bunn and White (1997 ) normalization 



In each realization, we generate a\ m from the power spectrum, multiply it by the harmonic- 
transformed DMR beam, G\ ( |Wright et al., 1994 ), transform G\a\ m back to a sky map, and add an 



instrumental noise realization to the map. Finally, we measure 466 modes of the bispectrum from 
each realization. We generate 50,000 realizations for one simulation; processing one realization 
takes about 1 second, so that one simulation takes about 16 hours of CPU time on a Pentium-Ill 
single processor personal computer. 

5.2.3 Normalized bispectrum 

The input power spectrum determines the variance of the bispectrum. Off-diagonal terms in the 
covariance matrix arise from incomplete sky coverage. When non-Gaussianity is weak, the variance 
is given by QLuo, 1994j ; |Heavens, 19981 ; ISpergel and Goldberg, 1999| ; |Gangui and Martin, 2000| ) 



(Bf lhl3 ) = (C h ) (C l2 ) (C h ) A hhl3 , (5.7) 

where A^/ 2 / 3 takes values 1, 2, or 6 for all Vs are different, two are same, or all are same, respectively. 
The brackets denote the ensemble average. 

The variance is undesirably sensitive to the input power spectrum; even if the input power 
spectrum were slightly different from the true power spectrum on the DMR map, the estimated 
variance from simulations would be significantly wrong, and we would erroneously conclude that 
the DMR map is inconsistent with Gaussian. It is thus not a robust test of Gaussianity to compare 
the measured bispectrum with the Monte-Carlo simulations. 

The normalized bispectrum, Bi x i^i z j (C^C^C^) 1 ^ 2 , is more sensible quantity than the bare 



bispectrum. Magucijo (1995 ) shows that the normalized bispectrum is a rotationally invariant 



spectrum independent of the power spectrum, as it factors out fluctuation amplitude in a/ m , which 

1/2 

is measured by . By construction, the variance of the normalized bispectrum is insensitive to 
the power spectrum, approximately given by A^; 2 ; 3 . 

One might wonder if the normalized bispectrum is too noisy to be useful, as the power spectrum 
in the denominator is also uncertain to some extent; however, we find that the variance is actually 
slightly smaller than A/ 1 / 2 ^ 3 . Figure ^]l] compares the variance of the normalized bispectrum, 
(Bf ihl3 / {C h C h C h )), with that of the bispectrum, (^ 2 lWa ) / ({C h ) (CQ (CQ). The top-left panel 
shows the case of full sky coverage. We find that the variance of the normalized bispectrum is 
precisely 1 when all l's are different, while it is slightly smaller than 2 or 6 when two Z's are 
same or all Vs are same, respectively. This arises due to correlation between the uncertainties in 
the bispectrum and the power spectrum, and this correlation tends to reduce the total variance 
of the normalized bispectrum. The rest of panels show the cases of incomplete sky coverage. 
While the variance becomes more scattered than the case of full sky coverage, the variance of the 
normalized bispectrum is still systematically smaller than that of the bare bispectrum. Thus, the 



5.2. MEASUREMENT OF BISPECTRUM ON THE DMR SKY MAPS 



87 



normalized bispectrum is reasonably sensitive to non-Gaussianity, yet it is not sensitive to the 
overall normalization of power spectrum. 

What distribution does the normalized bispectrum obey for a Gaussian field? First, even for a 
Gaussian field, the probability distribution of a single mode of Bi x i 2 i 3 is non-Gaussian, characterized 
by a large kurtosis. Figure |T2] plots the distributions of 9 modes of Bi t i 2 i 3 drawn from the Monte- 
Carlo simulations (solid lines) in comparison with Gaussian distributions calculated from r.m.s. 
values (dashed lines). We find that the distribution does not fit the Gaussian very well. Then, 

1 /2 

we examine distribution of the normalized bispectrum, B\ x i 2 i z j (C^C^C^) 1 . We find that the 
distribution is very much Gaussian except for l\ = l 2 = I3 = 2. Figure [53] plots the distributions of 
the 9 modes of the normalized bispectrum (solid lines) in comparison with Gaussian distributions 
calculated from r.m.s. values (dashed lines). The distribution fits the Gaussian remarkably well; 
this motivates our using standard statistical methods developed for Gaussian fields to analyze 
the normalized bispectrum. We could not make this simplification if we were analyzing the bare 
bispectrum. Furthermore, the central limit theorem implies that when we combine 466 modes the 
deviation of the distribution from Gaussianity becomes even smaller. 

Ferreira et al. (1998) claim detection of the normalized bispectrum at h = h = h = 16; 



Magueijo (2000 ) claims that the scatter of the normalized bispectrum for l\ = I2 — 1 and Z3 = I2 + 1 
is too small to be consistent with Gaussian. The former has analyzed 9 modes, while the latter has 
analyzed 8 modes. In the next section, we analyze 466 modes, testing the statistical significance 
of the non-Gaussianity with much more samples than the previous work. We calculate C/ from 
equation ( |5.6[ ), and then divide Bi 1 i 2 i 3 by {C^C^C^) 1 ^ 2 to obtain the normalized bispectrum. 



5.2.4 Testing Gaussianity of the DMR map 

We characterize statistical significance of the normalized bispectrum as probability of the measured 
normalized bispectrum being greater than those drawn from the Monte-Carlo simulations. We 
define the probability P as 



N ( /£ MR > If c 
P a = — K ' 



'dxF™ c {x), (5.8) 

where I a is the normalized bispectrum, A^ ota i = 50, 000 is the total number of simulated realizations, 
and a = 1, 2, 3, 4,..., 466 represent (h,l 2 ,k) = (2,2,2), (2,3,3), (2,2,4), (3,3,4),..., (20,20,20), 
respectively, with satisfying l\ < Z 2 < hi \h ~ h\ ^ h < k + lj, and h + h + h = even. F^f c (x) is 
the probability density distribution function (p.d.f ) of the simulated realizations for the normalized 
bispectrum, x = I^ C . The p.d.f is normalized to unity: J^° OD dxF^ lc (x) = 1; thus, P a lies in 
< P a < 1. 

By construction, the distribution of P a is uniform, if the DMR map is consistent with the 
simulated realizations, i.e., Gaussian. We give the proof as follows. By rewriting equation ( |5.8j ) as 
p a = /( /DMR ^ we ca i cu i ate the p . d . f of P ^ G(P a ), as 



/oo 
dy6[P a = f(\y\)]F™ R (y) 
-00 
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Figure 5.1: Variance of Normalized Bispectrum and Bare Bispectrum 

Variance of the normalized bispectrum, (Bf t t / (C; 1 C; 2 C; 3 )), in comparison with variance 
of the bare bispectrum, (Bf l t ) /((C; 1 ) (Cj 2 ) (Cj 3 )). These are derived from simulated 
realizations of a Gaussian sky. The top-left panel shows the case of full sky coverage, 
while the rest of panels show the cases of incomplete sky coverage. The top-right, bottom- 
left, and bottom- right panels use the 20°, 25°, and 30° Galactic cuts, respectively. 



:2. MEASUREMENT OF BISPECTRUM ON THE DMR SKY MAPS 



1 1 I 1 I 1 I 

■(2,2,2). 



1000 - 



1 1 I 1 I 1 I r 

(6,6,6) 



1 1 I 1 I 1 I r 
(10,10,10) 




-4-2 2 4 -4-2 2 4 -4-2 2 4 
B 111 /(<C 1 XC 1 XC 1 >) 1 / 2 

1 l 1 2 1 3 1 1 l Z X 3 



Figure 5.2: Distribution of Bispectrum 

Distribution of the bare bispectrum drawn from the Monte-Carlo simulations (solid 
lines). The 20° Galactic cut is used. B^i^/ ((C^) (C; 2 ) (C^)) 1 ^ 2 is plotted, where 
the brackets denote the ensemble average over realizations from the Monte-Carlo sim- 
ulations. The dashed lines plot Gaussian distributions calculated from r.m.s. val- 
ues. Each panel represents a certain mode of (I1J2J3) as quoted in the panels. 
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Figure 5.3: Distribution of Normalized Bispectrum 

Distribution of the normalized bispectrum (solid lines), Bi x i 2 i 3 / (G^C^Ci^) 1 ' 2 , in com- 
parison with Gaussian distributions calculated from r.m.s. values (dashed lines). The 
20° Galactic cut is used. The meaning of the panels is the same as in figure 5.2. 
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= f-\Pc 



F a MR (y) + Fr R (-y) 
^ DMR (y) + ^ DMR (- 



-y) 



(5.9) 



F^ c (y) + F^(-y) ' 

where F^ MR (y) is the p.d.f of the measured normalized bispectrum on the DMR map, y = 7^ >MR '. 
Our goal is to see if F^ MR (y) is consistent with the DMR data being Gaussian. It follows from 
equation (5J3) that G(P a ) = 1, when F^ M ^(y) = F^ c (y), regardless of the functional form of 
F~r (y)- In other words, the distribution of P a is uniform, if the distribution of the measured 
normalized bispectrum is the same as the simulated realizations. Since our simulation assumes the 
DMR map Gaussian, the P distribution, G(P), tests the Gaussianity of the DMR map. If the P 
distribution is not uniform, then we conclude the DMR data to be non-Gaussian. 

A Gaussian field gives equal number of modes in each bin of P. For example, it gives 46.6 modes 
in AP = 10% bin: 466 x G(P)AP = 466 x 0.1 = 46.6. If we detect the normalized bispectrum 
significantly, then we find that G(P) is not uniform, but increases rapidly as P increases. 

The top panel of figure (5^4] plots the P distribution for the three different Galactic cuts. We 
find that the distribution is uniform, and the number of modes in the bin (AP = 10%) is consistent 
with the expectation value for Gaussian fluctuations (46.6). 

To further quantify how well it is uniform, we calculate the Kolmogorov-Smirnov (KS) statistic 
for the P distribution in comparison with the uniform distribution. The bottom panel of figure 5^4 
plots the cumulative P distribution, for which we calculate the KS statistic. The probability of the 
distribution being uniform is 6.7%, 77%, and 52% for the three Galactic cuts, respectively. While 
we have confirmed that the normalized bispectrum at l± = I2 = I3 = 16 has P = 97.8% for the 20° 
cut and P = 99.3% for the 25° cut as similar to Ferreira et al. (1998 ), our result shows that it is 
well within statistical fluctuations. Thus, the properties of the normalized bispectrum of the DMR 
map are consistent with CMB being a Gaussian field. 



5.3 Model Fitting 

In this section, we fit predicted CMB bispectra to the measured normalized bispectrum. The 
predictions include the primary bispectrum from inflation and the interstellar foreground bispec- 
trum from the Galactic emissions. Then, we constrain a parameter characterizing the primary 
bispectrum. 



5.3.1 Primary bispectrum 



For the primary bispectrum from inflation, we consider weakly non-Gaussian adiabatic perturba- 
tions generated through non-linearity in slow-roll inflation. The simplest weak non-linear coupling 
gives 



$(x) = $ L (x) + / NL $E(x) - ($£(x) 



(5.10) 
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Figure 5.4: KS Test for Gaussianity with Bispectrum 

P is the probability of the CMB normalized bispectrum, 
Bhi 2 iJ {C h C h C h ) lr \ measured on the CORE DMR 53 + 90 GHz sky map, being 
larger than those drawn from the Monte-Carlo simulations. There are 466 modes in 
total. The thick dashed, solid, and dotted lines represent the three different Galac- 
tic cuts as quoted in the figure. The thin solid line shows the expectation value for 
a Gaussian field. The top panel shows the P distribution, while the bottom panel 
shows the cumulative P distribution, for which we calculate the KS statistic. The KS 
statistic gives the probability of the distribution being consistent with the expectation 
for Gaussianity as 6.7%, 77%, and 52% for the three Galactic cuts, respectively. 



P distribution (Eq. 
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where the square bracket denotes the volume average, and <&l( x ) is a linear Gaussian part of 
curvature perturbations. We call /nl the non-linear coupling parameter. 



Salopek and Bond (1990, 1991 ) and Gangui et al. (1994 ) show that slow-roll inflation gives this 



coupling; Pyne and Carroll (1996| ) shows that the second-order general relativistic perturbation 



theory gives this. The former predicts /nl as & certain combination of slope and curvature of a 
inflaton potential (<3?3 = — 2/nl in Gangui et al. (1994|) ). The latter predicts /nl ~ In 



chapter we have given the exact form of -Bi 1 z 2 z 3 f° r this model; we do not repeat it here, /nl is 
the parameter that we try to constrain by measuring the CMB bispectrum. 

Since the theoretical bispectrum assumes full sky coverage, we must correct it for the bias arising 
from incomplete sky coverage. We use an approximate correction factor for the bias, o bs/4vr, which 
we have derived in chapter 0. Moreover, the theoretical bispectrum must also be convolved with 



the DMR beam. We use the harmonic transform of the DMR beam, Gi, given in Wright et al. 



(1994 ). Hence, the observed bispectrum is related to the theoretical one through 



B?%l 3 = ^B^G h G h G l3 . (5.11) 

Note that o b s /47r = 1 — sin|6 cu t| for an azimuthally symmetric cut within certain latitude b CVL t, 
O obs /4vr = 0.658, 0.577, and 0.5 for |6 cut | = 20°, 25°, and 30°, respectively. 

5.3.2 Foreground bispectra from interstellar emissions 

Although we cut a fraction of the sky to reduce interstellar emissions from the Galactic plane, 



there should be some residuals at high Galactic latitude. |Kogut et al. (1996a|) have found signifi- 
cant correlation between COBE DMR maps at high Galactic latitude and COBE Diffuse Infrared 
Background Experiment (DIRBE) maps which mainly trace dust emission from the Galactic plane. 

The interstellar emissions are highly non-Gaussian. For example, the one-point p.d.f of the all- 
sky dust template map ( Schlcgcl et al., 1998|) is highly skewed. We find the normalized skewness, 



<(Ar) 3 )/((AT) 2 ) 3/2 ~ 51. Since these non-Gaussian emissions would confuse the parameter 
estimation of the primary CMB bispectrum, we take the effect into account. 

We estimate the foreground bispectra from interstellar sources by using two foreground template 
maps. One is the dust template map of [Bchlegel et al. (1998| ); the other is the synchrotron map of 
Haslam et al. (1981). Both maps are in the HEALPix format flGorski et al., 1998Q . 



We extrapolate the dust map to 53 GHz and 90 GHz with taking into account spatial variations 



of dust temperatures across the sky ( [Finkbeiner et al., 1999| ). We then cross-correlate the extrap 



olated maps with the DMR maps to confirm that the extrapolation is reasonable. We find that 
while the dust-correlated emission in the DMR 90 GHz map is consistent with the extrapolated 
dust emission, that in the DMR 53 GHz map is much larger than the extrapolated one. This is 



consistent with the anomalous microwave emission of Kogut et al. (1996a ). To take the excess 
emission into account, we multiply our extrapolated 53 GHz maps by factors of 3.66, 2.59, and 
2.45 for the 20° cut, the 25° cut, and the 30° cut, respectively. Note that we do not essentially 
need the correction for the excess emission, as it does not alter spatial distribution of the emission. 
Nevertheless, we do it for convenience of subsequent analyses. 
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We also extrapolate the synchrotron map to these two bands, assuming the spectrum of the 
source, T{y) oc v~ 2 ^. We do not need the extrapolation of the synchrotron template map either, as 
the extrapolation does not alter spatial distribution of the emission in contrast to the dust template 
map in which the extrapolation does alter it. We find no significant correlation between the DMR 
maps and the extrapolated synchrotron maps at both 53 GHz and 90 GHz. 

After coadding the extrapolated 53 and 90 GHz maps with the same weight as used for the 
DMR maps, we measure -B; 1 z 2 z 3 from the maps for the three different Galactic cuts, multiplying it 
by G^G^G^ to take into account the DMR beam. 

5.3.3 Constraints on non-linearity in inflation 

We simultaneously fit the primary, dust, and synchrotron bispectra to the measured bispectrum on 
the DMR map. We use the least-squares method based on a x 2 statistic defined by 

x 2 (/,) - E U MR - E W) (^L f ^ D ' MR - E f&) • 

aa' \ j J \ j J 

I 3 a is a model bispectrum divided by {icf^ c ^ (cf^ c ^j ^C^ C )) ^ ' wnere 3 represents a certain 
component such as the primary, dust, and synchrotron, fj is a fitting parameter for a component 
j, where /primary = /nl is the non- linear coupling parameter (Eq. (|4.1lD ). /dust and / syn c characterize 
amplitude of the foreground bispectra. 

C aa i is the covariance matrix of the normalized bispectrum calculated from the Monte-Carlo 
simulations: 

1 N 

Caa> EE j^-j £ (l^ ~ (l™ C )) (l^ ~ (l^)) , (5.13) 
i=l 

where iV = 50, 000 is the number of realizations. The bracket denotes an average over all realiza- 
tions, (^1^^ = N^ 1 Icx ■ Here, we have implicitly assumed the non-Gaussianity weak, so 
that we calculate the covariance matrix from Gaussian realizations. 

As we have observed in the previous section, distribution of I a is very much Gaussian; thus, 
X 2 {fj) should obey the x 2 distribution to good accuracy. Hence, minimizing x 2 {fj) with respect to 
fj gives the maximum-likelihood value of fj as a solution to the normal equation: 

h = E (^) 

i 

where 

^E'^L^- (5-15) 

aa' 

We estimate statistical uncertainties of the parameters using the Monte-Carlo simulations; we 
obtain parameter realizations by substituting for Iq MR in equation ( |5.14 ), 

Figure [5T5| plots the measured values of the non-linear coupling parameter, /nl, as well as the 
simulated realizations, for the three different Galactic cuts. The measured values are well within 



2J£ 



c 



-i 



rDMR 



(5.14) 
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the cosmic variance: we place 68% confidence limits on /nl as |/nl| < 1-6 X 10 3 , 1.7 x 10 3 , and 
2.0 x 10 3 , for the 20°, 25°, and 30° cuts, respectively. 



Figures and 15/7] plot constraints on /dust and /sync, respectively. There is no indication 
of either component contributing to the measured bispectrum significantly. The constraints are 
rapidly weakened as the Galactic cut is widened. 



5.3.4 Null test of the normalized bispectrum 



Using x 2 defined by equation ( 5.12 ), we can test Gaussianity of the DMR map. While the min- 
imization of x 2 (fj) gives constraints on the parameters, a value of x 2 (/j) t ens us goodness-of-fit; 
X 2 (0) tests a hypothesis of the bispectrum being zero. When x 2 (0) is either significantly greater 
or smaller than those drawn from the simulations, we conclude that the DMR map is inconsistent 
with zero bispectrum. 

X 2 (0) is similar to what several authors have used for quantifying statistical significance of non- 



Gaussianity in the DMR map QFerreira et al., 1998j ; [Magueijo, 2000| ; pandvik and Magueijo, 2000| ) 



They use only diagonal terms of the covariance matrix; however, the matrix is diagonal only on the 
full sky. As lack of sky coverage correlates one mode to the others, we should include off-diagonal 
terms as well. We did so in equation fl5.12|) . 

Figure |5.8| compares Xdmr(O) w ^h Xmc(O) f° r * ne different Galactic cuts. The measured values 
are Xdmr(O) = 475.6, 464.8, and 460.2 for the corresponding cuts, respectively, while (xmcW) = 
466. We find the probability of Xmc(°) bein § lar S er than Xdmr(°) to be p {Xmc > Xdmr) = 36.9%, 
49.7%, and 54.2%, respectively. 

We conclude that the DMR map is comfortably consistent with zero normalized bispectrum. 



We explain the claimed detection ( Ferreira et al., 1998 ) by a statistical fluctuation without invoking 



the eclipse effect of Banday et al. (2000| ). 



There is no evidence that the scatter of the normalized bispectrum is too small to be consistent 
with Gaussian, in contrast to the claim of Magueijo (2000| ) based on x 2 (0) derived from 8 modes. 



To clarify, our analysis does not reject the possibility that the CMB sky is non-Gaussian for only a 
small number of modes; however, in the absence of a theoretical motivation for limiting the analysis 
to a specific set of modes, we choose to treat all the bispectrum modes on an equal footing. [Sandvik 



and Magueijo (2000 ) claim that the non-Gaussianity found by Magueijo (2000| ) does not spread to 



other modes. This is consistent with our result. 

Incidentally, we plot in the figure |5,8| the x 2 distribution for 466 degrees of freedom, X466> m 
filled circles; we find that the distribution of Xmc(^) ^ s ver y similar to the x 2 66 distribution for a 
smaller cut as expected, while it becomes slightly broader for a larger cut for which the distribution 
of the normalized bispectrum at some modes deviates from Gaussian appreciably. Yet, we find 
that the 20° — 25° cuts reasonably retain the Gaussianity of the distribution of the normalized 
bispectrum. 
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Figure 5.5: Constraint on Non-linearity in Inflation 

DMR constraint on the non-linear coupling parameter, /nl, which characterizes non- 
linearity in inflation (Eq.( 4.11 )). The dashed, solid, and dotted lines represent the 
three different Galactic cuts as quoted in the figure. The thick vertical lines plot the 
measured values of /nl from the COBE DMR maps, while the histograms plot those 
drawn from the Monte-Carlo simulations for each cut. 68% confidence limits on /nl are 
|/nl| < 1-6 x fO 3 , 1.7 x 10 3 , and 2.0 x 10 3 for the three Galactic cuts, respectively. 
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Figure 5.6: Constraint on Dust Bispectrum 
DMR constraint on the amplitude of the interstellar dust bispectrum at high Galac- 
tic latitude, /dust- The meaning of the lines is the same as in figure 5.5 
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Figure 5.7: Constraint on Synchrotron Bispectrum 
DMR constraint on the amplitude of the interstellar synchrotron bispectrum at high 
Galactic latitude, / syllc . The meaning of the lines is the same as in figure |5.5| . 
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Figure 5.8: Null Test of Normalized Bispectrum 

Testing hypothesis of the normalized bispectrum, Bi t i 2 i 3 / (C/ 1 C/ 2 C; 3 ) 1 ^ 2 , being zero in 
the COBE DMR four-year 53 + 90 GHz sky map. The dashed, solid, and dotted lines 
represent the 20°, 25°, and 30° cuts, respectively. The thick vertical lines plot the 
measured x 2 (0)> while the histograms plot those drawn from the Monte-Carlo simu- 
lations. The filled circles plot the % 2 distribution for 466 degrees of freedom. 
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5.4 Discussion and Conclusions 

In this chapter, we have measured all independent configurations of the angular bispectrum on the 
CORE DMR map, down to the DMR beam size. Using the most sensitive sky map to CMB, which 
combines the maps at 53 and 90 GHz, we test the Gaussianity of the DMR map. 

We find that the normalized bispectrum, B^i^/ (C^C^C^) 1 ' 2 , gives more robust test of Gaus- 
sianity than the bare bispectrum, B^i^. We compare the measured data with the simulated 
realizations, finding the DMR map comfortably consistent with Gaussian. We explain the reported 
detection of the normalized bispectrum at h = h = h = 16 (Ferreira et al., 1998) by a statistical 



fluctuation as an alternative to the "eclipse effect" proposition made in Banday et al. (2000| ). 

We fit the predicted bispectra to the data, constraining the parameters in the predictions, which 
include the primary bispectrum from inflation and the foreground bispectra from interstellar dust 
and synchrotron emissions. We find that neither dust nor synchrotron emissions contribute to the 
bispectrum significantly. 

We have obtained a weak constraint on the non-linear coupling parameter, /nl, that charac- 
terizes non-linearity in inflation. We interpret the constraint in terms of a single-field inflation as 



follows. According to the analysis of non-linear perturbations on super horizon scales ( Salop ek and 



Bond, 1990), we can explicitly calculate /nl as 



5 f d 2 lnH \ 

/NL = ~2^G ) ' (5 - 16) 

where H is the Hubble parameter during inflation. When applying the slow-roll conditions to an 
inflaton potential V (</>), we have d In H/d(j) ~ (din V/d(j)) /2; thus, /nl is on the order of curvature of 
a slow-roll potential, implying that |/nl| should be smaller than 1 in slow-roll inflation. Therefore, 
the obtained constraint, |/nl| < 1-6 x 10 3 , seems too weak to be interesting; however, any deviation 
from slow-roll could yield larger |/nl|, bigger non-Gaussianity. 

The next generation satellite experiments, MAP and Planck, should be able to put more strin- 
gent constraints on /nl- In chapter we have shown that MAP and Planck should be sensitive 



down to |/nl| ~ 20 and 5, respectively. We find that the actual constraint from CORE (figure 5.5 ) 
is much worse than the estimate. This is partly due to different cosmology used for the model, 
but mainly due to incomplete sky coverage; the statistical power of the bispectrum at low multi- 
poles is significantly weakened by the Galactic cut. Since MAP and Planck probe much smaller 
angular scales, and their better angular resolution makes an extent of the Galactic cut smaller, 
the degradation of sensitivity should be minimal. Moreover, the improved frequency coverage of 
future experiments will aid in extracting more usable CMB pixels from the data. At this level of 
sensitivity, any deviation from slow-roll could give an interesting amount of the bispectrum, and 
MAP and Planck will put severe constraints on any substantial deviation from slow-roll. 

While we have explored adiabatic generation of the bispectrum only, isocurvature perturbations 
from inflation also generate non-Gaussianity flLinde and Mukhanov, 1997 ; [Peebles, 1997 ; pucher 



and Zhu, 19971 ) . They are in general more non-Gaussian than the adiabatic perturbations; it is 



worth constraining these models by the same strategy as we have done in this chapter. 



Chapter 6 

In Pursuit of Angular Trispectrum 



The angular trispectrum, the harmonic transform of the angular four-point correlation function, 
carries cosmological information which is independent of the power spectrum and the bispectrum. 

Several authors show that the weak lensing effect produces non-trivial connected CMB angular 
trispectrum or four-point correlation function on small angular scales ( Bernardeau, 1997 ; Zaldar-| 
riaga and Seljak, 1999 ; Zaldarriaga, 2000] ; Hu, 2001 ). On large angular scales, Inoue (2001a| ]b|) 



shows that the connected trispectrum is produced if topology of our universe is closed flat or closed 
hyperbolic. These effects do not produce the bispectrum, the angular three-point harmonic spec- 
trum, but the trispectrum. Hence, while we have not found significant bispectrum on the DMR 
map in chapter ||, we could find the trispectrum. 

The connected angular trispectrum contributes to the power spectrum covariance matrix. It 
increases the variance, and produces non-zero off-diagonal terms in the covariance matrix. Mea- 
suring the connected trispectrum thus constrains how much the connected trispectrum contributes 
to the power spectrum covariance. This is important to do. We have to know the power spectrum 
covariance matrix as accurate as possible, when we measure the power spectrum with better than 
1% accuracy, and determine many of cosmological parameters with better than 10% accuracy. 

So far, there has been no attempt to measure the angular trispectrum of the CMB anisotropy. 
In this chapter, we present the first measurement of the CMB trispectrum on the DMR map. We 
measure all independent terms of the angular trispectrum down to the DMR beam size, and test 
Gaussianity of the DMR data. 

This chapter is organized as follows. In 



6.1, we give our method of measuring the angular 



trispectrum from CMB sky maps. In § 6.2, we study statistical properties of the normalized 
trispectrum. In § |6.3| , applying the method to the DMR data, we measure the angular trispectrum. 
We then test Gaussianity of the DMR data with the normalized trispectrum. Finally, § S.4 concludes 
this chapter, 
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6.1 Angular Trispectrum 



The angular trispectrum comprises four harmonic transforms of the CMB temperature anisotropy 
field, ai imi ai 2m2 ai am3 ai im4 , where 



> AT (n) , 



ai m = d n — - — Y? m (n) . (6.1) 

•^obs 1 

^obs denotes a solid angle of the observed sky. A rotationally invariant angular averaged trispec- 
trum, T$*(L), takes the form (|Hu, 2001|) 



Tl£(L) = (2L + 1)££(-W W h t) 



k k L \ 

171% 771,4 — M J 

X O'hmi ®'l2m2 ^hms ^l^m^ i 

(6.2) 

where the Matrices denote the Wigner-3j symbol. By construction, l±, I2, and L form one tri- 
angle, while Z3, I4, and L form the other triangle in a quadrilateral with sides of l±, I2, h, and 



I4. L represents a diagonal of the quadrilateral. Figure 3.5 sketches a configuration of the an- 
gular trispectrum. When we arrange l\, I2, h, and I4 in order of l\ < I2 < I3 < I4, L lies in 
max(^2 — h, h — h) < L < min(/i + I2, I3 + I4). Hence, l\ = I2 and I3 = I4 for L = 0. Parity invari- 
ance of the angular four-point correlation function demands I1 + I2 + L = even and /3 + /4 + L = even. 
The angular trispectrum generically consists of two parts. One part is the unconnected part, 



the contribution from Gaussian fields, which is given by the angular power spectra ( Hu, 2001 ) 



+ (2L + 1) (C h ) (C l2 ) [{-l) l ^ +L 5 hh 5 l2h + S hU S hh ] . (6.3) 

For li < I2 < I3 < I4, the unconnected terms are non-zero only when L = or l± = I2 = h = h- 
We have numerically confirmed that our estimator given below (Eq.(6.4)) accurately reproduces 
the unconnected terms (Eq.(|6,3D) on a simulated Gaussian sky. 



The other part is the connected part whose expectation value is exactly zero for Gaussian fields; 
thus, the connected part is sensitive to non-Gaussianity. When none of Vs are same in 
one might expect the trispectrum to comprise the connected part only; however, it is true only on 
the full sky. The unconnected terms on the incomplete sky, which are often much bigger than the 
connected terms, leak the power to the other modes for which all Z's are different. We should take 
this effect into account in the analysis. 

Using the azimuthally averaged harmonic transform, e/(n) (Eq.(|3,9|)), we rewrite equation ( |6.2[ ) 
into a much more computationally efficient form, 



2L + 1 M=-L 



where 



^ S V^(o I d'n[e h (n)e l2 (n)]Yl M (n). (6.5) 
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Since t^M * s the harmonic transform on the full sky, we can calculate it quickly. This method 
makes measurement of the angular trispectrum computationally feasible even for the MAP data in 
which we have more than millions of pixels; thus, the methods developed here can be applied not 
only to the COBE DMR data, but also to the MAP data. 

For the DMR data for which the maximum I4 is 20, we have 21,012 non-zero trispectrum 
modes after taking into account the triangle conditions in a quadrilateral and parity invariance. 
Measurement of 21,012 modes takes about 5 second of CPU time on a Pentium-Ill single processor 
personal computer. 

We use Monte-Carlo simulations to estimate the covariance matrix of the angular trispectrum. 
We have described details of our simulations in chapter [|. We generate 5,000 realizations for one 
simulation; processing one realization takes about 5 second, so that one simulation takes about 
7 hours of CPU time on a Pentium-Ill single processor personal computer. In total, we run 6 
simulations, generating 30,000 realizations. 



6.2 Normalized Trispectrum 

In general, analytic forms of the angular trispectrum covariance, (T^^LjT^^L )^, are highly 



complicated ( Hu, 200~0| ); however, one can reduce them to rather simplified forms for h < I2 < 



h < h and L / (chapter |3|). The covariance of these terms on the full sky is diagonal, and the 
variance is given by equation ( 3.34| ), 



T !£( L )f) = ( 2L + 1) (Ch) (Ci 2 ) (C h ) (C h ) (1 + S hl2 + 5 hh + 6 hl2 6 l3h ) . (6.6) 

One finds that the variance of the trispectrum is very sensitive to the power spectrum normalization. 
A slight difference in the power spectrum normalization alters the variance of the trispectrum 
substantially. This makes a test for Gaussianity with the "bare" trispectrum very difficult, as it 
requires precise determination of the power spectrum normalization. 
We overcome the difficulty by normalizing the trispectrum as 

m . (6.7) 

This statistic, the normalized trispectrum, is analogous to the normalized bispectrum that we 
have used in chapter ||. As similar to the normalized bispectrum, the variance of the normalized 
trispectrum is insensitive to the power spectrum normalization, and systematically smaller than 
that of the bare trispectrum. The normalized trispectrum is thus reasonably sensitive to non- 
Gaussianity. 

Figure |0] compares the variance of the normalized trispectrum with that of the bare trispec- 
trum, for full sky coverage as well as for incomplete sky coverage. We have used l\ < I2 < h < U 
and L ^ terms, and calculated the variance from simulated realizations of a Gaussian sky. We 
confirm that the variance of the normalized trispectrum is systematically smaller than that of the 
bare trispectrum, and that the variance distribution becomes more scattered on the incomplete sky. 
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Figure 6.1: Variance of Normalized Trispectrum and Bare Trispectrum 
Comparison of the variance of the normalized trispectrum with that of the bare trispec- 
trum for Zi < l 2 < h < h and L / (group (a)). These are derived from simulated 
realizations of a Gaussian sky. The topdeft panel shows the case of full sky coverage, 
while the rest of panels show the cases of incomplete sky coverage. The top-right, bottom- 
left, and bottom- right panels use the 20°, 25°, and 30° Galactic cuts, respectively. 
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6.3 Testing Gaussianity of the DMR Map 



In this section, we test Gaussianity of the DMR data using the normalized trispectrum. Before 
performing the analysis, we should recall that even if CMB is exactly Gaussian, there are significant, 
non-zero unconnected trispectrum terms for L = or l\ = l 2 = I3 = Z4. We should analyze these 
configurations separately from the others for which the unconnected terms vanish. 

6.3.1 Classification of trispectrum configurations 

We divide the measured 21,012 modes into four groups: 

(a) l 2 £ l 3 , and L^O (16,554 modes) 

(b) l 2 = l 3 , h / Z 4 , and L / (4,059 modes) 

(c) l± = l 2 = l 3 = l 4 , and L ^ (209 modes) 

(d) L = (190 modes) 

Note that all the groups satisfy l\ < l 2 < h < h, max(/2 — h, U — h) < L < min(ii + l 2 , l 3 + I4), 
h + h + L = even, and I3 + U + L = even. 

The groups (a) and (b) are most sensitive to non-Gaussian signals, as for which the unconnected 
terms (Eq. (|6.3D ) vanish on the full sky. On the other hand, the groups (c) and (d) are dominated 
by the unconnected terms, and thus less sensitive to non-Gaussianity. These statements should, 
however, depend upon what kind of non-Gaussianity exists in the data. If non-Gaussian signals 
are subject to the groups (c) and (d) only, then the groups (a) and (b) become the least sensitive 
modes to the non-Gaussianity. 

On the incomplete sky, the unconnected terms also contaminate the groups (a) and (b) through 
the mode-mode coupling. We take this effect into account by using the Monte-Carlo simulations 
that use the same Galactic cut as the data. 

We discriminate between the group (a) and the group (b) in terms of the covariance matrix: 
the covariance matrix of the group (a) is diagonal, while that of the group (b) is not diagonal in 
L. We discriminate between the group (c) and the group (d) in terms of the statistical power: the 
group (d) has no statistical power of testing Gaussianity. The reason is as follows. For the group 
(d), the estimator given by equation (|]|) becomes ljjj(0) = {-l) h+h ^{2l l + 1)(2Z 3 + l)C h C h . 
The normalized bispectrum for the group (d) is thus just a pure number, 



This property holds regardless of Gaussianity. Even strongly non-Gaussian fields give exactly the 
same number. It thus follows from this result that we cannot measure the connected trispectrum 
for L = 0. 




(6.8) 



106 



CHAPTER 6. IN PURSUIT OF ANGULAR TRISPECTRUM 



This is unfortunate. It is the connected part of Tf^(0) that contributes to the covariance 
matrix of the power spectrum, 



(CC) - (Oi) (C,.) = + + (*(»)>„ • (6.! 



where ^T// ; ,(0)y is the ensemble average of the connected T// ; , (0). Even if we find the groups (a)-(c) 
consistent with Gaussianity, we can conclude nothing about the power spectrum covariance, unless 
we have a model for the connected trispectrum. We will discuss this point in § |6.4| . 

Henceforth, we analyze the groups (a)-(c) only, while we have used the group (d) to see if our 



code works properly. Our code reproduces equation (|6.8|) very well, and the numerical error is at 
most of order 10 -4 . 



6.3.2 Gaussianity test 

To quantify statistical significance of the measured trispectrum, we use a statistic P, which is 
the probability of the measured normalized trispectrum being greater than those drawn from the 
Monte-Carlo simulations: 



Pa 



N( J„ DMR > J™ C 



a 



-Wtotal 



(6.10) 



where J a denotes the normalized trispectrum, iVtotai = 30, 000 is the total number of the simulated 
realizations, and a represents a set of (Zi, I2, 13, 14, L)- The distribution of P a is uniform if the DMR 
map is consistent with Gaussian, for which there are equal number of modes in each bin of P. For 
example, when we calculate P a for all 21,012 a's, we expect a Gaussian field to give 210.12 modes 
in AP = 1% bin. If we detect the normalized trispectrum significantly, then the number of modes 
having higher P is much larger than the expectation value for a Gaussian field. 

In chapter [|, we have proven the P distribution uniform, if the DMR data are consistent with 
the simulated realizations (see Eq.(|5lj|)). We have also shown that this property holds regardless 
of the distribution function of J^ c . 

We calculate the KS statistic for the P distribution in comparison with the uniform distribution, 
to quantify how well the P distribution is uniform. We calculate the KS statistic for the groups (a)- 



(c) separately. Table 6A summarizes the KS-test results, and figures 6^2-6^ plot the cumulative P 
distribution, for which we have calculated the KS statistic. We find that the measured trispectrum 
is comfortably consistent with Gaussianity for all of the analyzed groups, (a)-(c). 

Since the groups (a) and (b) are zero for Gaussian fields, these groups provide the strongest 
constraint on generic non-Gaussian fluctuations. As we have done in chapter [5] for the angular 
bispectrum, if we have predictions for the CMB angular trispectrum (e.g., appendix |D|), then our 
measurement tests those predictions. 
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Figure 6.2: KS Test for Gaussianity with Trispectrum I 

Cumulative P Distribution (Eq. (|6.10 )), for which we calculate the KS statistic. P 
is the probability of the normalized trispectrum, T^^{L)/ [(2L + 1)C; 1 C; 2 C; 3 C; 4 ] 1 ^ 2 , for 
the group (a) (L ^ and h < l 2 < h < h), measured on the COBE DMR 
53 + 90 GHz sky map, being larger than those drawn from the Monte-Carlo simula- 
tions. There are 16,654 modes. From left to right panels, we use the 20°, 25°, and 30° 
Galactic cuts, respectively. The dashed lines show the expectation value for a Gaus- 
sian field. The KS statistic gives the probability of the distribution being consistent 
with Gaussianity as 5.4%, 12%, and 48% for the three Galactic cuts, respectively. 
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40 80 40 80 40 80 
P [%] P [%] P [%] 



Figure 6.3: KS Test for Gaussianity with Trispectrum II 

The same as figure 3.2 but for the group (b) (Za = Z3, h ^ h, and L 7^ 0). There are 
4,059 modes. The KS statistic gives the probability of the distribution being consistent 
with Gaussianity as 38%, 2.5%, and 5.2% for the three Galactic cuts, respectively. 
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Figure 6.4: KS Test for Gaussianity with Trispectrum III 

The same as figure 6.2 but for the group (c) (li = h = h = h and L / 0). There 
are 209 modes. The KS statistic gives the probability of the distribution being consis- 
tent with Gaussianity as 41%, 73%, and 63% for the three Galactic cuts, respectively. 
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Table 6.1: Gaussianity Test with Normalized Trispectrum 

Probability of the measured trispectrum being consistent with Gaussianity (the rightmost 
column) for the three different Galactic cuts. The probability is derived from the KS test for 
the P distribution in comparison with the uniform distribution. The group (a) comprises 
the modes of I2 ^ ^3, and L ^ (16,554 modes), the group (b) of I2 = h, h h, and L ^ 
(4,059 modes), and the group (c) of l\ = I2 = h = U, and L ^ (209 modes). 

group # of modes Galactic cut probability [%] 



(a) 


16,554 


20° 


5.4 






25° 


12 






30° 


48 


(b) 


4,059 


20° 


38 






25° 


2.5 






30° 


5.2 


(c) 


209 


20° 


41 






25° 


71 






30° 


63 



6.4 Discussion and Conclusions 

In this chapter, we have presented the first measurement of the CMB angular trispectrum on the 
COBE DMR sky maps. We have measured all the trispectrum terms, 21,012 terms, down to the 
DMR beam size. Since 190 L = modes have no statistical power of testing Gaussianity, we have 
used 20,822 L / modes to test Gaussianity of the DMR data, and found that the DMR map is 
comfortably consistent with Gaussianity. 

Our results do not directly constrain the connected trispectrum for L = 0, T//'/(0), which con- 
tributes to the power spectrum covariance through equation fl6.9p . We can thus conclude nothing as 
to whether the covariance matrix is diagonal on the DMR angular scales. Moreover, T}}(0) increases 
the power spectrum variance; we have no idea how much the contribution is. We need to use other 
statistics than the angular trispectrum to investigate the power spectrum covariance. Otherwise, 
we have to have a model for the connected trispectrum, and constrain Th,,(0) by measuring the 
other trispectrum configurations. 



One example for a trispectrum model is the one produced in a closed hyperbolic universe, tnoue 



(2001b ) suggests that the closed hyperbolic geometry produces non-zero connected trispectrum. In 
appendix |D|, we have derived an analytic prediction for the connected trispectrum produced in a 
closed hyperbolic universe (Eq.(p.9j)). For L = 0, we reduce the prediction to 



where 

Fhh =E p K%ThH^ 3 H (l£/ imi HI 2 ) (l W^)l 2 ) , (6-12) 
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and v = y/k 2 — 1 is a discrete wavenumber, is the primordial curvature perturbation, g^i{v) is 
the radiation transfer function, and ^i m (u) describes eigenmodes in a closed hyperbolic geometry. 
Hence, a single function, Fu<, determines the connected trispectrum completely. 

Fn> contributes to the power spectrum covariance directly. Using the prediction, we find the 
power spectrum covariance matrix in a closed hyperbolic universe (Eq. flD.lOj )), 

(CiC v ) - (Q) (C v ) = ^-y ((Ci) 2 + 2F n ) S w + 2F W . (6.13) 
The power spectrum is given by 

(Q)=E^(^M(l6mH| 2 ). (6.14) 
v 

Note that (Cj) (Cy) > Fu>. It follows from this equation that to constrain the contribution of the 
connected 2j/j/(0) to the power spectrum covariance, we need to constrain Fu>. 

In addition to L = modes, fortunately, a closed hyperbolic universe also produces the con- 
nected trispectrum for the group (c), l\ = I2 = I3 = 1 4 = I and L / 0. The prediction is 
TIHL)\ = 4(2L + 1)^;. We may constrain this term with our measurement of the group (c). 

/ c 

The group (c) is, however, the most noisy group, and the constraint from the DMR measurement 
is too weak to be useful yet. 

The most promising way to investigate the angular trispectrum in near future is to use the MAP 
CMB sky maps. The big advantage of MAP over COBE is the much higher angular resolution. 
The high-resolution measurement is important even on large angular scales, as we can minimize 
the effect of the Galactic cut. Since the rather big COBE Galactic cut has been the main cause of 
numerous complications in the data analysis, we expect that the M^4Pdata will measure the angular 
trispectrum with much better sensitivity, and with much smaller systematic errors than COBE. 
The MAP trispectrum will put strong constraints on several challenging non-Gaussian models: the 
non-Gaussianity induced from topology of the universe, the weak lensing effect, and so on, which 
are not probed by the angular bispectrum. 
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Appendix A 

Slow-roll Approximation 



In this appendix, we describe the slow-roll approximation, which has played a central role in the 
analysis of inflation dynamics. We then apply the approximation to the effective mass of scalar- field 



fluctuations, md(r), to show explicitly the approximation that we have used in equation ( 2.22 ), 

A slowly-rolling scalar field on a potential, V ((/)), is a key ingredient of a successful inflation 
model, for to achieve the accelerated expansion of the universe the kinetic energy of <j> needs to be 
smaller than the potential energy (see Eq. (|2.2|) ), 

£)' < V«). (A,) 

It then follows from this condition that the second-order time derivative, d 2 (j)/dt 2 , needs to be 
smaller than the potential slope, dV/dcj), 



(A.2) 



These conditions provide us an useful scheme of approximation, the slow-roll approximation. The 
approximation demands that the l.h.s's of the conditions be much smaller than the r.h.s's. 

The trace-part Einstein equation and the Friedmann equation give the exact relation between 
the time derivative of H and d(p/dt, 

(A.3) 



d 2 (p 




dV 


dt 2 


< 


dcf> 



dt \dt J 



Hence, we obtain another form of the slow-roll condition, \dH/dt\ < 4irGV((j)). In the slow-roll ap- 
proximation, the dimensionless variables such as V~ 1 (dcj)/dt) 2 , (Hdcp/dt) -1 (d 2 (p/dt 2 ), H~ 2 (dH/dt), 
and so on, are small order parameters which control the approximation. 

By the physical requirement, the accelerated expansion of the universe, the slow-roll conditions 
are defined by the time derivative with respect to the physical time, t. If we use the conformal 
time, r (dr = a~ l dt), then we have (f> = a(d(/)/dt), and 

4> = aH<j) + a 2 ^, (A A) 
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where the dots denote the conformal time derivative: x = d<p/dr. In contrast to d 2 6/dt 2 , 6 is not so 
small compared with dV/d<j) because of the first term. From this, we obtain a rule for the slow-roll 
analysis: evaluate the second- or the higher-order time derivative with respect to the physical time, 
not the conformal time. Using this rule, we can use the slow-roll approximation consistently in the 
conformal time coordinate. 

Using the slow-roll approximation, we derive the approximation to the effective mass of scalar- 
field fluctuations, mt, which we have done in chapter [2] (Eq. ([2.22D ). To derive 



2 _ H d 2 (6/H) a (*V dH\ 2 _2_ 



we begin with 



where 



J^ 2 - = H-6 +2 H\6-Hh (A ' 6) 



H . 66 a.6 ( d 2 6\ , . 

I = 2a 2 H 2 -2aH-a 2 ^, (A.8) 

6 dd 2 

\ ~ ^ 

Here, to calculate 0, we have used the Klein-Gordon equation for a homogeneous scalar field, 
6 + 2aH6 + a 2 (dV/d(j)) = 0. These equations are exact. 

We then neglect the higher-order slow-roll terms such as H 2 and (d 2 6/dt 2 )H, and obtain 

6 dr 2 \dd 2 dt J 

Next step is to relate H to the conformal time, r, through 

dt f da 1 f H 

^ = J^H = -7Ji-J da Wf- (A - n) 

By rewriting the conformal-time derivative in the second term with the physical-time derivative, 
/ da(dH I dt) / (a 2 H 3 ) , and neglecting the higher-order slow-roll terms, we obtain the first-order 
slow-roll correction to the conformal time, 

1 dH/dt , h , 

aH an 6 

and hence H 2 ~ (ar)" 1 — 2{dH/dt). Finally, substituting this for H 2 in equation ( A.1C| ), we obtain 

_ H d 2 {6/H) ^(fV + 9 dH\ 2 _2_ 

6 dr 2 [d6 2+ dt a r2- [AA6) 



Appendix B 

Wigner 3-j Symbol 



In this appendix, we summarize basic properties of the Wigner 3-j symbol, following Rotenberg et 



aL (19591) . The Wigner 3-j symbol characterizes geometric properties of the angular bispectrum. 



B.l Triangle conditions 

The Wigner 3-j symbol, 



h h h 
m\ m<i 777,3 



(B.l) 





is related to the Clebsh-Gordan coefficients which describe coupling of two angular momenta in 
the quantum mechanics. In the quantum mechanics, I is the eigenvalue of the angular momentum 
operator, L = r x p: L 2 Y/, m = 1(1 + l)Y/ m . m is the eigenvalue of the z-direction component of the 
angular momentum, L z Yi m = mY\ m . 
The symbol such as 

/ 7 1 i \ 

(B.2) 

describes coupling of two angular-momentum states, Li and L2, forming a coupled state, L3 = 
Li + li2- It follows from Li + L2 — L3 = that mi + 7772 — 7773 = 0; thus, the Wigner 3-j 
symbol ( p.l| ) describes three angular momenta forming a triangle, Li + L2 + L3 = 0, and satisfies 
777i + m 2 + m 3 = 0. 

Since Li, L2, and L3 form a triangle, they have to satisfy the triangle conditions, \Li — Lj\ < 
Lk < Li + Lj, where L, = |Lj|. Hence, 1%, I2, and I3 also satisfy the triangle conditions, 

\k-lj\<lk<k + lf, (B.3) 

otherwise, the Wigner 3-j symbol vanishes. The triangle conditions also include mi + 7772 + 7773 = 0. 
These properties may regard (/, m) as vectors, 1, which satisfy li + 12 + 13 = 0. Note that, however, 
L/l. 
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For h = I2 and 7713 = 0, the Wigner 3-j symbol reduces to 

("l) m ( 1 1 l ')= -fc=U'o- (B.4) 
V m -m J y/2l + 1 

In chapter ^ we have used this relation to reduce the covariance matrix of the angular bispectrum 
and trispectrum. We have also used this relation to reduce the angular trispectrum for L = (see 
Eq.flOTD). 



B.2 Symmetry 

The Wigner 3-j symbol is invariant under even permutations, 

( h h h \ = ( h h h \ = ( h h h 
\ mi tti-2 777,3 J \ 7773 mi 777.2 J \ m2 mi 

while it changes the phase for odd permutations if h + h + h = odd, 



(B.5) 



{ _ 1)h+ l 2 +lJ h h h \ (B6) 



mi ?7i2 mz 

h h h \ = ( h h h \ = ( h h h 
m2 mi 7773 I \ mi 7773 777.2 I \ ms tni mi 



(B.7) 



The phase also changes under the transformation of mi + tti-2 + ^3 — * —( m i + m 2 + m 3)> if 
h + h + h = odd, 

( h h k \ ( _ 1)h+h+ i 3 ( h h h \ (B 8) 

V mi T7T-2 TTT-3 / V — TTT-i — 777,2 — TTI-3 / 

If there is no z-direction component of the angular momenta in the system, i.e., rrii = 0, then the 
Wigner 3-j symbol of the system, 

I h h h \ 
/ ' 



(B.9) 



is non-zero only if li + I2 + ^3 = even. This symbol is invariant under any permutations of l{. 
In chapter ||, we have frequently used the Gaunt integral, G™^ 2 ™ 13 , defined by 



OT™" = / d 2 nY hmi (n)Y l2m2 (n)Y l3 rn 3 (n) 



' (2h + l) (2/ 2 + l) (2/ 3 + l) ( h h k 
4vr ^000 

( h h h \ 

\ mi TTT-2 TT13 / ' 



(B.10) 



to calculate the angular bispectrum. By definition, the Gaunt integral is invariant under both the 
odd and the even permutations, and non-zero only if l\ + I2 + I3 = even, ttt-i + TTI2 + 7773 = 0, 
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APPENDIX B. WIGNER 3-J SYMBOL 



and \li — lj\ < Ik < h + ly In other words, the Gaunt integral describes fundamental geometric 
properties of the angular bispectrum such as the triangle conditions. 

The Gaunt integral for rrii = gives the identity for the Legendre polynomials, 

2 



£y P h {x)Pi 2 {x)P l3 {x) = i h Q h Q k ) . (B.ll 



In chapter ||, we have used this identity to derive the bias for the angular bispectrum on the 
incomplete sky (Eq. Q3.4i| )). Here, we have used 



Y l0 (n) = J^-^P l (cos9). (B.12) 



B.3 Orthogonality 

The Wigner 3-j symbol has the following orthogonality properties: 



2,(2/3 + 1) , , }=Sm 1 m>$m 2 m', (B.13) 

^ — ' \ m-t mo mi I \ m, mL rrtn I ~ ~ 



and 



or 




(B.14) 





2/ 3 + l 

1. (B.15) 



h 

The orthogonality properties are essential for any basic calculations involving the Wigner 3-j sym- 
bols. Note that these orthogonality properties are consistent with orthonormality of the angular- 
momentum eigenstate vectors, and unitality of the Crebsh-Gordan coefficients, by definition. 



B.4 Rotation matrix 

A finite rotation operator for the Euler angles a, (3, and 7, D(a, (3, 7), comprises angular momentum 
operators, 

%A 7 ) = e -i«L. e -it>Lv e -*ri<*, (B.16) 

Since the Wigner 3-j symbol describes coupling of two angular momenta, it also describes coupling 
of two rotation operators. Using the rotation matrix element, D®, = (I, m! \D\ /, m), we have 

^l^Sl 2 = E(^ + i)E^t 3 f h h h \{ l \ h , \\ (B.i7) 

m 1 mi m 2 m 2 t-j z—^ m 3 m 3 y mi m2 771-3 / \ m l m 2 m 3 J 

In chapter ||, we have used this relation to evaluate rotationally invariant harmonic spectra. Note 
that the rotation matrix is orthonormal, 

E^Sm^S'm = <W" (B.18) 

m 



B.5. WIGNER 6-J SYMBOL 
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B.5 Wigner 6-j symbol 

The Wigner 6-j symbol, 

describes coupling of three angular momenta. We often encounter the Wigner 6-j symbol, when we 
calculate the angular bispectrum which has more complicated geometric structures flGoldberg and 



Spcrgel, 1999| ), The angular trispectrum also often includes the Wigner 6-j symbol ( |Hu, 2001|) . 




The Wigner 6-j symbol is related to the Wigner 3-j symbols through 

^1 h h 

h ^2 ^3 

—nil m 2 m 3 / I m l ~ m 2 






(B.20) 



In appendix |C], we use this relation to derive the angular bispectrum from isocurvature fluctuations 
in inflation (Eq.( |C.7| )). By using equation ( |B.15| ), we also obtain 



(_l)Ji +l'2+l'J L ) '2 '3 I = (_l}rn> 1+ m' 2 +m' 3 





(B.21) 



Appendix C 

Angular Bispectrum from 
Isocurvature Fluctuations 



In this appendix, we derive the angular bispectrum from isocurvature fluctuations generated in 



inflation. The mechanism of generating isocurvature fluctuations we consider here is that of Linde 
and Mukhanov (1997| ): a massive-free scalar field, a, oscillating about a = during inflation. 



Quantum fluctuations of a produce Gaussian fluctuations, 5a. Since there is no mean cr-field 
because of the oscillation about a = in the model, the model predicts density fluctuations which 
are quadratic in 5a: 

5p a ~ m 2 a5a + m 2 (5af = m 2 {5af. (C.l) 

5p a is thus non-Gaussian. Moreover, since the energy density of a does not dominate the universe 
during inflation, 5p a does not perturb the spatial curvature, being isocurvature density fluctuations. 

After inflation, 5p a may produce the spatial curvature perturbations in the Newtonian gauge, 
$, through the evolution. If 5p a becomes dominant in the universe at some point, then the linear 
perturbation theory gives = ^(5p a /p a )(a/a eq ) in the radiation era, and <I> = l(5p a /p a ) in the 
matter era ( Kodama and Sasaki, 1984j ). <I> then produces CMB fluctuations through the Sachs- 



Wolfe effect, AT/T = —2$. Since 5p a is non-Gaussian, $ is also non-Gaussian, so is AT/T. 

Our goal in this appendix is to calculate the CMB angular bispectrum from the <I>-field bispec- 
trum. We start with writing <I> as a Gaussian- variable-squared in real space, $(x) = r/ 2 (x) — (r/ 2 (x)) , 
where rj is Gaussian. We have subtracted the mean from ensuring (<I>(x)) = 0. Transforming 
<3?(x) into Fourier space, we obtain 

Hk)= /(^^( k + P)^(p)-(2^ (3) (k)(r? 2 (x)). (G2) 
Using the rj power spectrum, Pri(k), we write the $ power spectrum, P$(k), as 

Mk) = 2 lwf Pr > {p)Pr > (|k + p|) • (c - 3) 
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If we use a conventional power-law spectrum, P$(k) oc A:™ -4 , then we find P^{k) °c fc(™~ 7 )/ 2 . On 
the other hand, quantum fluctuations of a massive-free scalar field give P a (k) oc fc- 3 + 2m2 /( 3 - ff2 ) ) 
where is the Hubble parameter during inflation (chapter ||). It then follows from $ oc 5p a that 
P v (k) oc P CT (fc) oc fc- 3 + 2m2 /(3H 2 ) ; and n = 1 + 4m 2 /(3if 2 ); thus, the model predicts a tilted "blue" 
spectrum ( [Linde and Mukhanov, 1997| ). 
The <3?-field bispectrum is 

($(k 1 )$(k 2 )$(k 3 )) = (2^) 3 5( 3 )(k 1 + k 2 + k3)^|^P,(p) 

x [P v (|ki + p|) P v (|k 2 - p|) 

+P 7? (|k 2 + p|)P r ,(|k3-p|) 

+P„(|k 3 + p|)P„(|k 1 -p|)]. (C.4) 

The delta function has appeared in the formula to satisfy the triangle condition, ki + k 2 + k3 = 0. 
Since we will work in harmonic space eventually, it may be useful to expand P^dkipl) into 
harmonic space, 

P T ,(\k±p\) = ^P$\k,p)Y LM (l)Y£ M (p). (C.5) 

LM 

The next step is to project the 3-dimensional $-field bispectrum onto the 2-dimensional CMB 
angular bispectrum. We write the harmonic coefficients of the CMB anisotropy, a; m , with $ and 
the radiation transfer function, g^i^k), as 

, f r/ 3 k 

a lm = 4vr(-iy J -—^g^Y^lc). (C.6) 

On large angular scales, the Sachs- Wolfe effect gives gTi(k) = [k(ro — r dec)]! where To and r^ec 
are the present-day conformal time and the decoupling-epoch conformal time, respectively. On 
small angular scales, we need the full radiation transfer function for isocurvature fluctuations. 

We then calculate {ai imi ai 2m2 ai 3m3 ) . A key point on calculations is to expand the delta function 
in equation ( |C.4j ) into harmonic space with the Rayleigh's formula. After lengthy calculations, we 
obtain the CMB angular bispectrum from isocurvature fluctuations in inflation 



(a hmi a hm2 a l3m3 ) = S^ 2 ™ 3 - j Q r 2 dr ^ p 2 dp P v (p) 



OO ("OO 



Ejrhh 
l 2 l[ 



hhfo 
L 



y x v, 2 L 



x- r~kfdh P^\k u p)g Th (k 1 )j ll (k 1 r)(-i) h -^ 
vr Jo 1 

x- f k\dk 2 P^fe^W^tMH)^ 2 
vr Jo ' 2 

2 f°° 

x- / k 2 dk 3 g T i 3 (k 3 )ji 3 (k 3 r) 
7T Jo 

+ (i 2, 2 3) + (1 3, 2 1)] , (C.7) 
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where G™lJ! } 2rn3 and F]}]?^ represent geometric structures of the bispectrum. They are written 
with the Wigner 3-j and 6-j symbols, 



Qmim Wh + 1) (2Z 2 + 1) (2k + 1) h h h) h h h \ fcgl 
/l/2 ' 3 V 4tt I j \ mi m 2 m 3 J ' 1 ' ' 



■ h l 2 l 3 = (2l[ + l)(2l' 2 + l)(2L + l) ( h l 2 l z \ 1 f h h h \ 
l 'A L ~ 4vr ^ j \ 1' 2 l[ L J 

x f l'l 1*2 h\(h l[ L\fl 2 l> 2 L\ y> l+l , 2+L , } 



Appendix D 



Angular Trispectrum in Closed 
Hyperbolic Universe 



In this appendix, we derive the angular connected trispectrum in a closed hyperbolic universe. If 
topology of our universe is infinite fiat, then we can expand the harmonic coefficients of the CMB 
anisotropy into Fourier series, 

a lm = M-iyj <Z>(k)g Tl (k)Y lm (k), (D.l) 

where <3?(k) is the primordial curvature perturbation, and gTi(k) is the radiation transfer function. 
Hence, statistical properties of ai m take over statistical properties of $ directly; if $ is Gaussian, 
then ai m is also Gaussian. 

If topology of our universe is closed hyperbolic, then the relation between a/ m and <3? becomes 
Qlnoue, 2001b|) 

aim = $ (")STf MOmC"). (D-2) 



where v = \Jk 2 — 1 is a discrete wavenumber, and new expansion coefficients, £/ m (z^), describe 
eigenmodes in a closed hyperbolic geometry. The complex conjugate of £; m is given by £* m (^) = 
(— l) m 6-m( l/ )- Since a closed hyperbolic universe is globally inhomogeneous, £z m (f) is also a 
function of observer's positions and orientations in the universe. Inoue (2001b| ) have shown that 



this property makes ^ m (^) behave as if it were a Gaussian random number with the covariance 
matrix diagonal, 

{ilm{v)£l>m>{v')) = (l&mHI 2 ) (~ 1)™' fa' 5 m -m' <W , (D-3) 

and ^|C/m( zy )| 2 ^) u ~ 2 ■, where the bracket denotes the ensemble average over observer's positions 
and orientations in the universe. Hence, in addition to statistical properties of £j m also affect 
statistical properties of ai m . If <I> is Gaussian, then it follows from the relation among a; m , <E>, and 
£i m (Eq.(p.2j)) that a; m comprises two independent Gaussian random numbers. 
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APPENDIX D. ANGULAR TRISPECTRUM IN CLOSED HYPERBOLIC UNIVERSE 



The angular power spectrum, Cj, is given by 



(Ci, 



aw, 



(D.4) 



where P${y)5 uu i = (<J>(z/)3>*(i/)). The ensemble average for $ is taken over initial conditions. We 
should keep in mind that the ensemble average for is taken over all possible observer's positions 
and orientations in the universe. Hence, there exist those observers who measure C\ different 
from (Ci) substantially. In other words, there exists theoretical uncertainty in (Q) arising from 
uncertainty as to where we are. 

We calculate the angular trispectrum from equation (D.2) as follows. We begin with 

(^Zimi Oi2"i2 ^hms 0,141x14 

all v 

X ($(l/i)$(l^)*(l*)*(l/4)) 



If $ is Gaussian, then this equation yields 9 terms. First, we reduce the second line to 



(D.5) 



($(i/ 1 )$(i/ 2 )$(zy 3 )$(i/ 4 )) 



P$ (f l)-P$ 3)^1/1 V 2 



+P$(^l)P$(^ 2 )(^ 1 ^3 "^2^4 

We then evaluate the first three terms in the trispectrum: 

E 9T« 1 (l / l)5T/ 2 (^l)5T«3( l/ 3)9T«4( zy 3)Pl'(l / l)Pl'(l / 3) 

X [(^imiC^O^mz^l)) (6 3 m 3 (^3)64^14(^3)) 
+ (6 1 m 1 (^l)6 3 m 3 (^3)) (£/ 2 m 2 (f 1)^4(^3)) 
+ (6imi(^l)6 4 (^3)) (62 

= (C h )(C h )(-ir^S hh S hh S mi 



+ 5 UlU4 5 U2U3 ) . 



(D.6) 



-m2 w m3— m,4 



~^~Fi 1 i 2 { 1) 3 4 {Si 1 i 3 Si 2 i 4 5 mi — m3 5 m2 — mi ^5i 1 i 4 5i 2 i 3 5 rni — rn4 5 m2 — m3 ), 



where 



,hmi (y) 



(D.7) 



(D. 



The first term in the r.h.s. of equation ( D.7| ) is the unconnected term, while the second and third 
terms are the connected ones. By collecting all the connected terms, we obtain 

i a h mi Q>l2rri2 ^i 3 m 3 0^47714 / c 
= Fihi -0 {^Iil3^l2l4^mi— m->,^m2— m 4 ^h^^hh^mi — 014^2- m 3 ) 

+(2 <->3) + (2 «->4). 

Actually, the functional form of this equation is exactly the same as the unconnected terms, if 
we replace 2F/./. with (CQ (C/ . V Hence, by substituting 2i*y. for (Cjl) (C; / in the formula of 
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the unconnected trispectrum (Eq.( |3.25| )), we find the angular connected trispectrum in a closed 
hyperbolic universe 



(T l ^(L)) c = (-l) l ^yj(2h + l)(2/ 3 + 1) (2F hl3 ) S hh S l3k 5 L0 

+(2L + 1) (2F hl2 ) \(-l) l ^ +L S lll3 5 hh + S hh S hh ] . (D.9) 



Unfortunately, the unconnected terms always contaminate the connected terms produced in a 
closed hyperbolic universe: there is no non-Gaussian signal in the modes for which the unconnected 
terms vanish. Since what we measure is the sum of the unconnected terms and the connected terms, 
(Cjj) ) + %Flih> we nave to subtract (QJ (Cljj from measurement to obtain 2Fij j . Although 
Fiii- is difficult to measure because of the contamination, if we detect non-Gaussian signals that 
exist only in these limited modes (L = 0, or h = fa = fa = fa), then it may suggest that the universe 
is closed hyperbolic. 

The connected trispectrum for L = contributes to the power spectrum covariance matrix 
(Eq.( |3.12"|) ). Using our prediction for the connected trispectrum (Eq.( p.9j )) for L = 0, we obtain 
the power spectrum covariance matrix in a closed hyperbolic universe, 



The variance is amplified, and the off-diagonal terms appear in the covariance matrix; these prop- 
erties may be used to search for a signature of closed hyperbolic geometry in the CMB angular 
power spectrum covariance. 




(D.10) 



